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1.0 ABSTRACT 

This report describes research performed at the 
University of Missouri-Columbia to devise and design a lower 
cost inverter-induction motor drive for electrical pro- 
pulsion of passenger vehicles. A two-phase inverter-motor 
system is recommended. It is predicted to provide com- 
parable vehicle performance, improved reliability and nearly 
a 10% cost advantage for a high production vehicle because 
of the reduction in total parts count, decreased total 
rating of the power semiconductor switches and somewhat 
simpler control hardware compared to the conventional three- 
phase bridge inverter-motor drive system. The major disad- 
vantages of the two-phase inverter-motor drive are that the 
two-phase motor is larger and more expensive than a three- 
phase machine, the design of snubbers for the power switches 
is difficult because motor lead and bifilar winding leakage 
inductances produce higher transient voltages, and the 
torque pulsations are relatively large because of the 
necessity to limit the inverter switching frequency to 
achieve high efficiency. 

The conclusions of this research are based on analyses 
and simulations. An actual model of the two-phase system 
must be constructed and evaluated to prove the advantages 
claimed for the proposed system. The most challenging 
engineering design task will be to design the inverter, 
motor and snubber circuits to minimize transient voltages 
with high system efficiency. 
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2.0 EXECUTIVE SUMMARY 

This report describes research performed at the 

University of Missouri-Columbia to devise a lower cost in- 
duction motor drive for use in electrical propulsion of 
passenger vehicles. The work was done during the period 
from January through August of 1984 in two phases. 

Phase I - Concept and Feasibility Study 
( January-May ) 

Phase II - Detailed Design and Performance Analysis 
(June-August ) 

As a result of Phase I, three preferred candidate 
systems were identified. 

A. Half Wave Three-Phase Inverter-Motor (three 
power switches, three pairs of bifilar wound 
stator windings) 

B. Four Switch Three-Phase Inverter-Motor (four 
power switches - two legs of conventional three- 
phase inverter, conventional three-phase 
machine, motor line returned to center- tap of 
dc source) 

C. Two-Phase Inverter-Motor (four power switches, 
two pairs of bifilar wound stator windings) 

After a detailed analysis of the candidate systems at 
the beginning of Phase II, the two phase inverter-motor was 
selected as the lowest cost approach to meet the electrical 
passenger vehicle system requirements. 
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Detailed analyses and digital computer simulations were 

/ t t * \ < \ “f i * 

carried out to design the inverter circuit and the two- 

; <_• : 

phase induction motor required. The two-phas'e system has 

the following advantages over the conventional three-phase 

- -1 ’ ’ 

inverter-motor drive: 

1. The number of power semiconductor switching 
devices is reduced from six to four. 

2. The current rating of each power switch is 65% 

i ,! 

of that required for each switch in the 
equivalent three-phase system. 

3. The total number of components is estimated to 

, -f -*-r * ^ » 

be about 70% of the number required for the 

- ' i. 1 ' 

three-phase circuit. This also implies higher 

- i ^ J f * 

reliability and lower manufacturing asembly 
cost . 

4. The motor transient temperature rise is reduced 
due to the increased volume of copper. 

5. Somewhat simpler control hardware is possible 
because of the common emitter circuit so that 
isolation is not required in the base drive 
circuits . 

6. The total system cost is estimated to be about 
90% of that for the three-phase system. 

The overall efficiency is almost the same for the 


two-phase system and the three-phase system - approximately 
91% at the 12kw, 80% speed operating condition. 



5 


The major disadvantages of the two-phase approach are 
increased motor size and cost, more difficult snubber design 
for transient voltage suppression, and greater torque 
pulsations . 

For a high production electric passenger vehicle, the 
two-phase inverter-motor system developed on this research 
should provide significant cost advantages for the electric 
vehicle application when compared to the conventional three- 
phase bridge inverter-induction motor drive. It is 
important to note that the conclusions of this research are 
based on analysis and simulations. An actual model of the 
two-phase system must be constructed and evaluated to prove 
the advantages claimed for the proposed inverter-motor 
system. The most challenging engineering task will be to 
design the inverter, motor and snubber circuits to minimize 
transient voltages while retaining high system efficiency. 



3.0 TWO-PHASE BIFILAR INDUCTION MOTOR/INVERTER SYSTEM 


3 . 1 Specifications 

The inverter-motor drive system requirements are 
defined by the contractual motor specifications in Appendix 
A. These system requirements are as follows: 

Inverter input: 

240V battery 
Motor output: 

a. Shaft power output 

1) 30kW for 1 hour with 20% duty 
cycle at maximum speed 

2) 60kW for 3 minutes with 30/% duty 
cycle at 80% of maximum speed 

3) 16kW for 3 minutes with 30% duty 
cycle at 10% of maximum speed 

b. Maximum motor speed between 8,000 and 
16,000 rpm 

c. Two or four pole motor, totally enclosed, 
convection cooled with 40° C maximum ambient 
air and 4000 hours operating life at 80% of 
maximum speed with 12kW output 

Fig. 1 is a graphical representation of the motor 
requirements. It illustrates one method of satisfying the 
specifications. However, later in this section, an 
alternate motor design is described which results in a 
reduction in motor size and weight below that implied by 
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Fig. 1. For Fig. 1, it is assumed that a constant volts/Hz 
is applied to the motor from zero speed to 37.5% speed. In 
this range the air-gap flux is essentially constant, im- 
plying a constant torque capability and a linearly increas- 
ing output power with speed. At 37.5% speed, the maximum 

voltage available is required, and so the motor voltage is 
assumed constant from 37.5% speed to 80% speed. It should 
be noted that whenever the equivalent battery voltage is 
less than 240V, it is assumed that PWM control of converter 
output voltage is used to provide this voltage reduction. 
Since the output power requirement is only 30kW at 100% 
speed, it is possible to reduce the motor applied voltage 
from 80% to 100% speed. The capability implied by Fig. 1 
would meet the specification requirements since the three 
specified operating points lie on the power vs. speed curve, 
which defines the power capability at each speed. 

The inverter-motor system also is to be designed to 
minimize the function 

F = Cost + 10 X Weight - 120 X Efficiency (1) 

where each of the parameters on the right side of this 
relationship is for the total inverter and motor system ex- 
pressed as follows: 

Cost = retail cost in dollars 
Weight « total weight in kilograms 

Efficiency * % efficiency at 80% of maximum motor 


speed with 12kw output 
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As a result of this research, a two-phase inverter/ 
motor drive is proposed. The remainder of this report de- 
scribes the system design to tt)e„ specifications. In 

addition, the report explains why .the system proposed is con- 
sidered the preferred induction motor approach for this 
application. , 

3.2 System Description ... , 

3.2.1 Circuit Configuration and. Operation ... , 

During the first phase. of._,,;the research, several 
alternative inverter-induction motor, combinations were con- 

sidered. The following preferred candidate systems were 
selected for more detailed, analyses : v 

A. Half-wave inverter driven three-phase bifilar 
wound induction moto.r (Appendix B, Fig. B-l). 

B. Four switch inverter . driven three-phase 
induction motor (Appendix B,. Fig. B-2). 

C. Two-phase inverter, driven, induction motor with 

two pairs of bifilar wound .stator windings as 
shown in Fig . 2. • ...... 

As a result of an evaluation of each of these systems, in- 
cluding simulations to determine.; the peak, current in the 
power switches at the 60kW - 80% speed operating point, the 
two-phase approach shown in Fig. 2 was selected. It has two 
significant advantages compared to • the . alternative candidate 
systems: (1) the lowest total -power -switch rating (the 

power switch peak current times.. its peak voltage times the 
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number of switches required); (2) it is one of the two 
alternatives (Fig. B-l-b is the other) which may be arranged 
to have a common emitter connection with transistor switches 
or common cathodes if thyristors are used. This second ad- 
vantage makes it possible to simplify the control hardware 
since gate or base drive isolation is not required. The 
primary disadvantages of the scheme in Fig. 2 are that the 
motor size and weight is increased by 10%-20% compared with 
the three-phase bridge system (Fig. B-3), the power switch 
snubber design is more difficult since line and bifilar 
winding leakage inductances produce additional transient 
voltages, and there are relatively large torque pulsations 
because of the necessity to limit the inverter switching 
frequency to achieve high efficiency. Appendix C contains 
an idealized comparative analysis of the two-phase and 
three-phase inverter arrangements. 

The two-phase bifilar induction motor has only two 
stator phases, each of which includes a pair of tightly 
magnetically coupled coils. Fig. 2 shows the circuit con- 
figuration. Fig. 3 is an alternate circuit. However, Fig. 
2 is the preferred approach since it simplifies the control 
hardware as mentioned previously. With transistor switches, 
all emitters are common, and all cathodes are common with 
thyristor switches. Thus, isolation is not required for the 
base drive or gate drive circuits. 

Four switching devices and four free-wheeling diodes 
are required in the circuit of Fig. 2. The center-taps of 
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the stator windings are connected to the positive side of 
the battery. The operation of the inverter is as follows. 
Suppose switch SI is closed, and current is flowing from the 
battery through stator winding Al. When SI is opened and S2 
is closed, the current through stator winding Al is 
transferred to the bifilar winding A2. As a result, 
free-wheeling diode D2 conducts, and the stored energy in 
the winding returns to the battery. After this diode 
current reduces to zero, the current starts to flow through 
switch S2. The current in stator windings Bl and B2 is 
controlled in a similar manner. However, the conduction 
intervals for switches S3 and S4 are displaced 90 electrical 
degrees from those for Si and S2. Controlling the two-phase 
currents in this fashion, a rotating magnetic flux is 
produced in the machine. 

3.2.2 PWM Control of Inverter 

PWM control is employed in the inverter system so that 
the current through the stator windings can be shaped to the 
optimum waveform. In the light load case, the current is 
shaped to approach a sinusoidal waveform as shown in Fig. 
4-a. The switching pattern is controlled such that the 
stator current always stays within a given band from the 
sinusoidal function. With heavy load, a current waveform 
approaching a square wave is delivered to the motor to 
provide maximum power with minimum peak motor current as 
shown in Fig. 4-b. 
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With a microprocessor controller, it is also possible 
to control the PWM voltage waveform and the applied 
frequency to maximize efficiency at the 80% speed, 12 kW 
operating point. This might be accomplished by storing the 
required control program in memory, or possibly by an 
on-line efficiency optimization technique. 

3 . 3 Motor Design 

3.3.1 Operating Points 

As indicated in section 3.1, three specific operating 
points are required. 

16kW output at 10% speed (3 minutes, 30% duty) 
60kW output at 80% speed (3 minutes, 30% duty) 

30 kW output at 100% speed (1 hour, 20% duty) 

The two-phase bifilar induction motor proposed must satisfy 
the operating conditions at each of these points. 

3.3.2 Nature of Induction Motor 

The polyphase induction motor maximum power output as 
a function of inverter frequency is given by 


p »ax * C 1 ‘r> f 


( 2 ) 


where 


constant = 


3P 


8ttL 


lr 


lr - 


P 
L 

u 
V A 

f 4 


number of poles 
rotor leakage inductance 
rms value of stator voltage 
source frequency 
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When an induction motor is connected to an inverter, a 
similar maximum output is achieved, neglecting the effects 
of harmonics. Since the output voltage and frequency of the 
inverter may be adjusted, the maximum output power given by 
(2) can be changed, depending on the speed control strategy. 
If V/f is maintained constant over the low frequency range, 
the maximum power is proportional to the frequency (or rotor 
speed). After maximum stator voltage V is reached, the max- 
imum power decreases proportional to 1/f since V is constant 
and f increases. In Fig. 5, generalized curves for the max- 
imum power P , rated power P_, stator current I_, 

luclX K t* 

stator voltage V, and slip frequency f are plotted as a 
function of rotor speed. 

The rated output power as a function of slip frequency 
f SR , which may change at different source frequencies, is 


p s c (-)^f f 
R '■'2 ' f r SR 


(3) 


where 


C 


2 


A 3P 

= constant = 


A 

R 2 = rotor resistance 
A 

f SR = rated slip frequency at different source 
frequencies 

3.3.3 Motor Design Procedure 

For these discussions, the following variables are 


defined: 
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Peak outnut power; Corresoondinq 
Deak stator current; Corresoondinq 
Deak stator voltaoe and slip 
frequency f sn ,„ 

Rated outnut nower; Corresoondinq 
rated current; Corresoondinq 
rated stator voltage and rated 
si in frequency f $R 





Fiq. 5 Generalized motor characteristics 
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Base frequency f fi A 


Base output P g A 

Base current I B A 


Base voltage V g A 
Peak output power P p A 
Peak current I p A 


the design frequency for the 
motor (the frequency 
when V reaches a maximum) 
the rated output at the base 
frequency 

the current which produces 
rated output at the base 
frequency 

the motor voltage at base 
frequency 

the maximum output power at 
f B 

the current which achieves 


There is a simple relation between these variables: 


P L. f 

P _P smax 

P B I B f SR 


at f B 


(4) 


where f gmax is the maximum slip frequency which produces 
P P at f B 

Fig. 6 shows three possible motor designs. One 
possible design to satisfy the specifications is shown by 
the B curves in Fig. 6, where the maximum power output is 
16kW at 10% speed and 60kW at 80% speed. However, this is 
not the optimum motor design, because the peak current is 






Fig. 6 Three- 



20 


excessive at low speed and the motor weight is not a min- 
imum. It is very important to have minimum peak current 
since the power switch device ratings are determined by the 
maximum current through the devices. The peak current can 
be reduced for a given output power by lowering the base 
frequency, but this increases the size of the machine. 
Therefore, the base frequency should be increased, resulting 
in a smaller machine. Then, the voltage should be increased 
at low frequency to saturate or overexcite the machine. 
This permits higher power output with reduced peak current 
and minimum motor size and weight. The optimum machine 
design should have nearly equal peak currents at 10% speed. 

Constant V/f operation means nearly constant air-gap 
flux. Thus, by increasing the V/f ratio, oversaturated air- 
gap flux can be achieved. Since output power is propor- 
tional to the product of air-gap flux and rotor current, 
larger air-gap flux helps to decrease the rotor current, and 
as a result, lower stator current can deliver the same out- 
put power. This mechanism is shown in Fig. 7. Case I on 

this figure corresponds to case B on Fig. 6. In the digital 
simulation program, this saturation effect can be easily 

checked by sensing the magnetizing current. Thus, by a 
trial and error method, oversaturation of air-gap flux is 
simulated, and stator current can be reduced until the 

allowable saturation flux is reached. 

Also, when the motor is operated at the overloaded 
output, the temperature rise due to the greater loss is a 





22 


very important factor for machine design. In summary, the 
induction motor design is optimized when the following 
conditions are satisfied: 

1. 60kW output is designed to be maximum output 

power at 80% speed. 

2. 16kW output at 10% speed is achieved by over- 

saturation of air-gap flux so that the peak 
current is almost the same as that for 60kW at 
80% speed. 

3. 30kW output is equal to the rated output at 

100% speed. 

4. The temperature rise at the overload condition 

is within an allowable range. 

5. The weight of the motor is minimized. 

3.3.4 Computer-Aided Design 

The design of the two-phase induction motor involves 
four steps. A flow chart is given in Fig. 8. 

Step 1 ; Choose the Base Operating Point 

Since a machine can be almost uniquely designed at a 
base operating point (base frequency, base voltage and base 
output power), the base point should be selected carefully 
so that the designed machine is the optimum one satisfying 
all required operating conditions. Usually one must come 
back to Step 1 after checking the given conditions in Step 
3. 

Maximum frequency and pole number 

For a given power rating, the size and the weight of 
an ac machine decrease as the rated frequency increases. On 


\ 



Requirements (Output Power - Speed) 



Fig. 8 Design flow chart 
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the contrary, the iron loss increases at high frequency, and 
the mechanical bearings become very expensive at high rotor 
speeds. A rotor maximum speed of 12,000 rpm is chosen be- 
cause this is believed to provide a reasonable practical 
trade-off between losses, weight and bearing cost. For a 
given power rating and maximum speed, as the pole number in- 
creases, the maximum ac source frequency increases, which 
causes greater iron loss and switching loss. Thus, a two- 
pole machine is chosen, which sets the maximum ac source 
frequency at 200Hz. 

Base voltage, base frequency and base output power 

The base operating condition to establish the motor 

design is chosen as the 30kW output, 40% speed point shown 

in Fig. 9. The base voltage is selected as follows. With 

the PWM technique, the peak value of the fundamental 

component of the output voltage can change from zero to a 

maximum of (4/tt )(V dc ) ( where V DC is the battery voltage 

of 240V. The motor is designed such that the applied 

voltage as a function of speed changes as shown in Fig. 10. 

The V/Hz (Volts/ Hertz) is constant from zero to 40% speed. 

From 40% to 80% speed, the V/Hz is reduced with speed as the 

voltage is linearly increased to a maximum of 306V peak at 

80% speed. Then the voltage is fixed at the maximum value 

from 80% speed to 100% speea. A base voltage of 241V peak 

at 40% speed is chosen by an iterative technique. Since, 

roughly speaking, the induction motor output power is 

proportional to the square of the applied voltage and 

2 

inversely proportional to the applied frequency ( p 00 V"/f), 



25 


(kw) 

60 


40 


20 


0 


Outout Power 


- - - - - -- — 

© 

60 kw ( 3 minutes 


0 

proposed base noint 

30 kw ( 1 

• 


16 kw ( -3 minutes ) - — 

• 

12 kw for efficiency 

: 


0 


20 40 60 80 

Fiq. 9 Output power-speed requirement 


100 ( % ) 
speed 


Note: • are specification operating points. 


hour ) 





27 


the proposed voltage variation in Fig. 10 implies that the 
rated output power at different speeds is 
7.5kW at 10% speed 
30kW at 40% speed 
38kW at 80% speed 
30kW at 100% speed 

The 16kW, 10% speed point requires about twice the 
rated output power that would be delivered with the low 
speed range constant V/Hz of 241V/80HZ. (This implies 
approximately 60V at 10% speed, or 20Hz.) Thus, the motor 
voltage must be increased above 60V at 10% speed to provide 
16kW. When the voltage is increased at a given frequency, 
the air-gap flux increases, which drives the machine further 
into saturation. This may result in excessive motor current 
due to the reduction of the magnetizing reactance. Thus, 
the motor voltage required for maximum power at 10% speed 
must be high enough to drive the machine into saturation 
without excessive current. However, if the V/Hz is too high 
at 10% speed, it will be too low at 80% speed because the 
maximum voltage available with a 240V source will be reached 
at a lower speed. A low V/Hz at 80% speed will mean that 
the 60kW output cannot be achieved. 

Step 2 : Design the Machine at the Base Operating 

Point 

The detailed design of the induction motor, both in 
mechanical and electrical characteristics, is included in 
Appendix D. When the base operating point is given, the 
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machine is essentially uniquely defined. The final design 

is determined in Appendix D after several iterative trials. 

The resulting machine parameters are as follows: 

External diameter D 11.5 inches 

o 

Internal diameter D 5.75 inches 

Rotor length L 6.3 inches 

Weight W * 190 lb <= 86kg) 

Stator resistance R s 0.0274 ft 

Rotor resistance R r 0.0213ft 

Magnetizing inductance L m 13.86mH 

Stator leakage inductance L^ s 0.25mH 

Rotor leakage inductance L^ r 0.27mH 

Step 3 : Check the Output Power by Simulations, 

Considering Saturation Effects 
Since the machine is designed for the base operating 
point, the maximum output power should be checked by digital 
simulations, including the saturation effect of the air-gap 
flux. Appendix G describes the detailed simulations. 

Fig. 11 shows the simulation results for the output 
power-speed relations at 10%, 40%, 80% and 100% speed of the 
finally designed motor, when the stator current is limited 
to 245A peak. All of the specified operating points are 
satisfied. At the maximum power output of each curve, the 
current is almost a square wave since this current waveform 
results in more power output than with sinusoidal current 
for the same peak current. At the 16kW, 10% speed point. 


the greatest saturation of the magnetic core occurs. The 
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saturation effect of the core is implemented in the 
simulation study by changing the magnetizing and leakage 
inductances of the machine. If the magnetizing current 
exceeds the allowable limit, then the magnetizing and 
leakage inductances are replaced with adjusted values, which 
are calculated from the B-H curve. This modification is 
easily implemented in Step 2. If the magnetizing current 
still exceeds the allowable limit, which occurs for very low 
slip frequency, then this is not considered a permissible 
operating condition. 

At the 60kW output, 80% speed operating point, under- 
saturation of the core occurs because the voltage is not 
high enough to produce rated flux. In this case, when the 
magnetizing current in the simulation is below the allowable 
limit, the magnetizing and leakage inductances are recalcu- 
lated in Step 2 using the B-H curve. 

Step 4 : Check the Temperature Rise 

Most of the heat generated in the machine comes from 
the copper loss, which is proportional to stator resistance 
times the square of the current. According to the simula- 
tion results, the rms currents at the given operating points 
are as follows: 

1) 16kW at 10% speed (3 mins., 30% duty) I i rms “ 200A 

2) 60kW at 80% speed (3 mins., 30% duty) I i rms = ^ 75A 

3) 30kW at 100% speed (1 hour) I lrms = * 15A 

Assuming a conventional rotor fan type construction, 

the air flow is reduced at lower rotor speeds so that the 
worst-case cooling condition is the 10% speed point. If the 
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iron and mechanical losses are considered negligible, then 
the highest transient temperature rise occurs for the 16kW, 
10% speed operating point. Since this condition continues 
for three minutes with a 30% duty cycle, the load pattern is 
assumed to be 16kW for three minutes, followed by no output 
for seven minutes, repeated several times. The temperature 
rise is estimated by solving the heat balance differential 
equations (Appendix E). The estimated temperature rise is 
50 °c The bifilar windings help to reduce the temperature 
rise because of the increased volume of copper. The 
temperature rise is higher if the same heat loss is genera- 
ted in a smaller volume of copper. 

In the calculation of the temperature rise, the most 
pessimistic conditions are considered. First, static air is 
assumed for the heat dissipation. Next, the heat transfer 
to the core and frame is neglected. Thus, the actual 
temperature rise should be somewhat less than 50° C. 

3.3.5 Efficiency 

In accordance with the JPL specifications, the most 
important operating point from tne standpoint of efficiency 
is the 12kW, 80% speed point. The theoretical approach to 
optimize the efficiency of an induction machine implies that 
the magnetizing current should be equal to the rotor 
current.* But, since the V/hz is reduced at the 80% speed 


* T. W. Jian, N. L. Schmitz, D. W. Novotny, "Characteristic 
Induction Motor Slip Values for Variable Voltage Part Load 
Performance Optimization," IEEE PAS, Jan. 1983, pp. 38-46. 
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point, the magnetizing current does not reach the optimized 
value because of the undersaturation of the air-gap flux. 
Thus, the best that is possible is to make the magnetizing 
current its maximum value and to control the current to get 
12kW output at 80% speed. Simulation results, performed by 
an iterative method changing the slip frequency and the 
amplitude of the sinusoidal motor current, give 98.6% effi- 
cieny at a slip frequency of 1.16Hz (0.725% slip). The cur- 
rent and output power are shown in Appendix G. The pulsating 
power is +65%. These simulations neglect the iron loss, 
stray loss, and mechanical losses. If these losses are 
included, the efficiency is expected to be about 95%. 

3 . 4 Inverter Design 

3.4.1 Switching Device Requirements 

In Step 3 of the design procedure in the previous 
section, three given output power-speed operating points 
were checked by simulations. The maximum peak current 
through the stator winding is 245A when either 16kW at 10% 
speed or 60kW at 80% speed is required. Since the forward 
blocking voltage of the device is two times 240V, the total 
VI rating of the devices for the two-phase inverter will be 
(240 X 2 X 245) x 4 = 470kVA. 

From the viewpoint of the power dissipation in the 
switching devices, 16kW output at 10% speed produces the 
maximum switching loss. From the simulations, the switching 
frequency is 3.4kHz at 16kW, 10% speed, and 160Hz at 60kW, 
80% speed. 
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In summary, the switching device requirements are 
Peak current = 245 A 
RMS current = 200. 5A 
Maximum switching frequency = 3.4kHz 
Peak voltage {not including transients) = 480V 
3.4.2 Power Circuit 

As a result of a study of the available transistor and 
thyristor devices to meet the switching requirements for 
this application, MG300M1FK1 Toshiba power transistors are 
recommended. These transistors have 1000V blocking capabil- 
ity, 300A dc rating and 1.6kW maximum power dissipation ca- 
pability. Thus, this Toshiba transistor provides a nominal 
dc voltage safety factor of 1000/480 = 2.08 and an rms cur- 
rent safety factor of 300 X /F = 2.12. (Note: This does 
not include the transient voltages across the transistors. 
The transient voltages are discussed in the next section.) 

The two-phase inverter requires four power switch 
modules, each of which includes a Darlington power transis- 
tor, a freewheeling diode, a snubber circuit and a base 
drive. However, only one heat sink is required since the 
MG300M1FK1 transistor is an isolated package. Also, no base 
drive isolation is necessary for the circuit of Fig. 2 with 
common emitters. 

The two-phase and three-phase inverter circuits both 
can deliver power to the motor load or take power from the 
source. This is accomplished by the proper control of the 
inverter switching instants relative to the ac voltage 
generated by the load machine, which produces a negative 
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slip. Appendix G.5 gives simulation results for regenera- 
tive operation. 

The detailed inverter design is described in Appendix 
F. 

3.4.3 Transient Voltage 

The two-phase system does present a more difficult 
snubber design problem because of the leakage inductances 
that will be present in the bifilar windings and because of 
the inductances in the motor leads from the inverter to the 
motor. It will be necessary to use special techniques to 
minimize the motor lead inductance, such as using parallel 
closely spaced flat conductors between the inverter and the 
motor or packaging the inverter essentially on the motor 
housing. 

The amplitude of the voltage spike primarily depends 
on the equivalent leakage inductance of the bifilar winding 
and the lead inductance between the inverter and motor. The 
energy stored in these inductances must be dissipated in the 
snubbers. Based on estimated practical values of total 
leakage and lead inductance, the voltage spike phenomena can 
be simulated by the digital computer, as shown in Appendix 
F. However, a precise design of the snubbers required can 
be accomplished only through experiment with an actual model 
of the inverter. 

As shown in Appendix F, the peak transistor voltage is 
about 900V with 4 pH total lead and leakage inductance. 
However, if the inductance is reduced to lyH, the peak tran- 
sistor voltage is approximately 600V. This would reduce the 
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transistor voltage safety factor from 2.08 (as given in the 
previous section) to about 1.6, which is still a quite 
reasonable voltage safety factor. 

3.4.4 Control Circuit 

A microprocessor controller will be used, which con- 
sists of an MPU, ROM, RAM, timer, I/O port and current 
sensors which can provide a variety of speed and torque 
controls. In either case using PWM, the speed or torque 
controller generates the frequency and amplitude commands 
for the stator current. This current is regulated so as to 
stay within a hysteresis band of the reference command by 
PWM control of the inverter. Also, the switching frequency 
is monitored and limited to less than 3.4kHz. The maximum 
fundamental inverter frequency command is 200Hz. 

The use of the microprocessor permits the implementa- 
tion of very sophisticated controls to provide optimal 
efficiency at any operating point and protection of the 
drive system. 


4.0 FINAL DESIGN AND EXPECTED PERFORMANCE OF A IWO-PHASE 
INVERTER/BIFILAR INDUCTION MOTOR SYSTEM 


4.1 Machine Data 

External diameter 
Internal diameter 
Rotor length 
Weight 

Output rating 
Number of phases 
Base frequency 


D =11.5 inches 
o 

D = 5.75 inches 
L = 6.3 inches 
W = 190 lb (= 86kg) 
30kW 
2 


80Hz 
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Maximum speed 
Base voltage (rms) 

Base current (rms) 

Stator resistance 
Rotor resistance 
Magnetizing inductance 
Stator leakage inductance 
Rotor leakage inductance 
4.2 Inverter Data 

Power switching devices 


12,000 RPM 

170V 

110A 

R s = 0.0274 & 

R r » 0.0213 n 

L = 13 . 86mH 
m 

L ls = 0 . 25mH 
L lr = °‘ 27mH 

(4) Toshiba 
MG300M 1FK1 
Transistors 


Input voltage 
Output voltage 
Output current 
kVA rating 
kW rating 

Maximum output frequency 
Maximum switching frequency 
Estimated efficiency 
at 12kW, 80% speed operating 


240V 

240V (peak or rms) 

245A peak; 200A rms 

96kVA 

60kW 

200Hz 

3.4kHz 


point 95.4% 

4 . 3 Inverter/Motor Performance 

The digital computer simulations described in Appendix 
G provide the following performance data at the specified 
operating points: 



16.4kW output at 10% speed 

Slip frequency = 1.7Hz (slip * 8.50%) 
Switching frequency = 3.4kHz 
Stator peak current* 245A* 

Stator rms current = 200A 
Pulsating torque = +20% 

60.5kW output at 80% speed 

Slip frequency = 3.2Hz (slip = 2%) 
Switching frequency = 160Hz 
Stator peak current = 245A** 

Stator rms current = 175A 
Pulsating torque = +30% 

13.1kW output at 80% speed 

Slip frequency =1.16 Hz (slip = 0.725%) 
Switching frequency = 1.6 kHz 
Stator peak current = 103A*** 

Stator rms current = 53. 2A 
Pulsating torque = +65% 


*The current is a PWM controlled wave form with a square 
wave envelope 205+40A. 

**The inverter voltage is a square wave producing an 
unmodulated current waveform. 

***The current is a PWM controlled sinusoidal waveform with 
an envelope of + 38A. 



Motor efficiency (only with copper loss) = 98.6% 
Motor efficiency (with 3.4% iron and stray 
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mechanical loss) = 95.2% 

Inverter efficiency (with snubber, on-state and 
switching loss) = 95.4% 

Total efficiency = 91% 

4.4 System Schematic Diagram 

The complete schematic diagram of the inverter/motor 
drive system is shown in Fig. 12. 

As mentioned previously, transient voltages are 
produced by lead inductances and leakage inductances due to 
imperfect coupling between the bifilar motor windings. With 
the input capacitor connected across the inverter dc bus 
very close to the inverter transistors and motor windings, 
the dc lead inductance is minimized. Snubber circuits are 
included to suppress voltage transients due to the imperfect 
coupling between bifilar windings and the inductance of the 
leads between the inverter transistors and motor windings. 
This also implies that the inverter transistors must be 
mounted very close to the motor. 

The design of the snubbers is discussed in Appendix F. 
However, the precise snubber configuration can be determined 
only through experimentation with an actual model of the in- 
verter. A final optimum arrangement may include varistors 
and additional capacitors. 
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The base drive circuits are critical for reliable power 
transistor operation with minimum switching loss. Current 
sensors are required to detect each motor winding current 
for the PWM controlled current operations as well as for 
overcurrent protection. Sophisticated control strategies 
can be implemented in the microprocessor software. 

5.0 COMPARISON OF TWO-PHASE BIFILAR AND THREE-PHASE 
BRIDGE INVERTER MOTOR SYSTEMS 
5 . 1 Advantages and Disadvantages of Two-Phase System 

Advantages 

A. The number of power switching devices is reduced 
from six to four. 

B. The current rating of each power switch is 65% of 
the value required for the three-phase bridge 
inverter (Fig. B-3). This has been verified by 
simulation as shown in Table II of the next 
section. 

C. The total number of components is estimated to be 
about 70% of the number required for the three- 
phase system. Since four transistor switches are 
used instead of six as in the thre-phase bridge 
inverter, this would imply about two-thirds as 
many components. However, some circuits are not 
repeated for each switch and the snubbers are 
more complex for the two-phase inverter. The 
component count reduction also implies higher 
reliability and lower manufacturing assembly 


cost. 
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D. The motor transient temperature rise is reduced 
due to the increased volume of copper. 

E. The control hardware is somewhat simpler due to 
the common emitter arrangements so that isolation 
is not required in the base drive circuits. 

F. The total component cost for the two-phase system 

is estimated to be about 90% of that for the 
three-phase system. This cost reduction cannot 
be proven until prototype equipments are manufac- 
tured. The cost will be reduced because of the 

reduced total number of components, the elimina- 
tion of the need for base drive circuit isola- 

tion, and the reduced peak current in the tran- 
sistor switches. However, the cost will be 

greater because of the higher peak transistor 
voltages required, the need for very effective 
snubbers to minimize transient voltages, the re- 
quirement for a somewhat larger motor than for 
three-phase, and the necessity to minimize motor 
bifilar winding leakage inductances and lead in- 
ductances in the inverter current switching 
paths. It is presently estimated that these 

tradeoffs will result in reducing the two-phase 
system cost to about 90% of that for a three- 

phase inverter/motor drive. 

Disadvantages 

A. The size, weight and cost of the two-phase motor 
is greater than for the three-phase machine. It 
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is estimated that an equivalent three-phase motor 
would have 85% of the weight and 90% of the cost 
of the two-phase bifilar machine. 

B. The design of the snubbers for the power switches 
is more difficult since motor lead inductances 
and bifilar winding leakage inductances produce 
additional transient voltages. 

C. The torque pulsations are relatively large 

because of the necessity to limit the inverter 
switching frequency to achieve high efficiency. 

5 . 2 Simulation Results 

Digital simulations to compare the two-phase and three- 
phase systems are presented in Appendix G. Table I lists 
the parameters used for the simulations. The results of 
these simulations are summarized in Tables I and II. 

Fig. 13-a shows the stator voltage and current per 
phase, and Fig. 13-b shows the instantaneous output power 
for the two-phase machine. Fig. 13-c and Fig. 13-d show the 
corresponding variables for the three-phase machine. 

The output power or torque pulsations are greatest for 
the two-phase system at the 60kW operating point since a 
square wave voltage is supplied from the two-phase inverter. 
However, at the 12kW operating point, with sinusoidal 
current control, the torque pulsations of the two-phase 
system are nearly as large because the inverter switching 
frequency must be limited to achieve high efficiency. 
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Table I Machine Parameters 


Bifilar Two-Phase Induction Motor 
(at 160Hz ) 

R $ =0.0274 

R '=0.0213 
r 

L 1s =2.94x10" 4 H 

L lr =3.23xl0 _4 H 

L =26.8x10" 3 H 
m 

Weight=86kg 


Three-Phase Induction Motor 
(at 320Hz) 

R s =0.00578 

R '=0.0038 
r 

L 1s =5.95x10' 5 H 

L ] r =7.0xl0 _5 H 

L=3.32xl0“ 3 H 

m 

Weight=73kg 


Note: A four pole three-phase motor was used for this comparison because 
more commercial data was available to us on this type of machine. However, 
we believe a two pole three-phase machine is somewhat preferable. 
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Table II Comparison of Two-Phase and Three-Phase Inverter-Motor System 


80% speed, 60kW output (for peak current comparison) 


Case I : 

Two-Phase Induction Motor 

Peak current I p =245A 

rms current I =175.8A 
rms 

slip = 2.0% 

P out = 60 ' 5kW 
frequency = 160.0Hz 

Case II: 


Three-Phase Induction Motor 

Peak current I p =375A 

rms current * rms =24 3A 

slip = 1.25% 

P Qut ■ 60.2kW 

frequency = 329.8Hz 


80% speed, 12kW output (for efficiency comparison) 


Two-Phase Induction Motor 

Peak current I = 103A 
P 

rms current I rms = 53. 2A 

Slip = 0.725% 

’’out ■ 13.1KW 
frequency = 160Hz 

efficiency (only copper loss)=98.6% 


Three-Phase Induction Motor 
Peak current I p =105A 

rms current I rms = 66A 

slip = 0.8% 

P out - 12.5kW 

frequency = 329.8Hz 

efficiency (only copper loss)=98.9% 
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VOLTAGE & CURRENT (80% S.60KW, MOTOR OPERATION) 



Fiq. 13-a Two-ohase 






Fig. 13-c Three-phase 


anrou 
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OUTPUT POWER (80% SPEED MOTOR OPE? 
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Fig. 13-d Three-phase 
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5 . 3 Performance Function Comparison 

Machine: The following data are used 


1) 

Cost of two-phase machine 

$360* 


Cost of three-phase machine 

$325* 

2) 

Weight of two-phase machine 

86kg 


Weight of three-phase machine 

73kg 

3) 

Efficiency at 12kW, 80% speed 



Two-phase machine 

95% 


Three-phase machine 

95% 

4) 

Performance function F (F=C + 

10W - 120E) 


Two-phase machine 

-10. IK 


Three-phase machine 

-10. 3K 

Inverter: The following data are 

used 

1) 

Estimated manufacturing cost 

of two-phase 


inverter 

$1000 


Estimated manufacturing cost 

of three- 


phase inverter 

$1200 

2) 

The weight and efficiency of 

the two-phase 


and three-phase inverters are 
be essentially the same. 

> estimated to 

3) 

Performance function (neglecting the weight and 
efficiency terms) 


Two-phase inverter 

$1000 


Three-phase inverter 

$1200 


♦These are best judgment estimates assuming at least 10,000 
units/year. 
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Inverter-Motor System : (Including only the para- 

meters in the Performance Function which are differ- 
ent. ) 

Two-phase system $360+$1000+$10X86=$2220 
Three-phase system $325+$1200+$10X73=$2255 

Thus, the two-phase system has only a very slight 
advantage over the three-phase system based on the 
Performance Function. 

6.0 CONCLUSIONS 

The proposed two-phase inverter-motor drive has the 
following advantages over the traditional three-phase bridge 
inverter-motor system: 

1. The number of power semiconductor switching 
devices is reduced from six to four. 

2. The current rating of each power switch is 65% of 
that required for each switch in the equivalent 
three-phase system. 

3. The total number of components is estimated to be 
about 70% of the number required for the three-phase 
system. This also implies higher reliability and 
lower manufacturing assembly cost. 

4 . The motor transient temperature rise is reduced 
due to the increased volume of copper. 

5. Somewhat simpler control hardware is possible 
because of the common emitter circuit so that 
isolation is not required in the base drive 
circuits . 



51 


6. The total component cost is estimated to be about 
90% of that for the three-phase system. 

The overall efficiency of the two systems is about the 
same and the Performance Function in the JPL Specifications 
is also nearly equal for the two-phase and three-phase 

systems . 

The major disadvantages of the two-phase system are: 

1. The weight and cost of the two-phase machine are 
greater. It is estimated that the three-phase 
machine would have 85% of the weight and 90% of 
the cost of the two-phase bifilar machine. 

2. The snubber requirements for the inverter power 

transistors are much more severe for the 

two-phase inverter since motor lead inductances 
and bifilar winding leakage inductance produce 

additional transient voltages. 

3. The torque pulsations are relatively large because 
of the necessity to limit the inverter switching 
frequency to achieve high efficiency. 

For a high production electric passenger vehicle, the two- 
phase inverter-motor system developed on this research should 
provide significant cost advantages for the electric vehicle 
application when compared to the conventional three-phase bridge 
inverter-induction motor drive. It is important to note that 
the conclusions of this research are based on analysis and simu- 
lations. An actual model of the two-phase system must be con- 
structed and evaluated to prove the advantages claimed for the 
proposed inverter-motor system. The most challenging engineering 
task will be to design the inverter, motor and snubber circuits 
to minimize transient voltages while retaining high system effi- 
ciency . 


APPENDIX A 


MOTOR SPECIFICATIONS PROVIDED BY JPL 



EXHIBIT I 


Induction Motor Specification 

The following is a specification for an ac induction 
motor to be used in conjunction with a Pulse Width Modulated 
(PWM) Inverter and a 240 volt battery for propulsion of 
electric vehicles using a fixed ratio gear reduction. 

I. BACKGROUND 

Numerous traction drives have been developed using in- 
verter driven induction motors. In all cases the inverter 
converts dc battery power into three-phase power of con- 
trolled frequency and voltage. While control algorithms 
vary, with all schemes, the electrical frequency is main- 
tained close to the mechanical frequency and the voltage to 
frequency ratio is maintained below a critical maximum (lest 
magnetic saturation occur). 

Compared with dc systems, such ac drives are 
advantageous in that inherent four quadrant performance is 
achieved, motor cost is greatly reduced, motor weight and 
losses are approximately halved while motor maintenance is 
virtually eliminated. Until recently, however, these 
advantages were offset by cost, weight, and inefficiency 
problems of the inverter. 

II. SYSTEM DESIGN PHILOSOPHY 

With traction drives, where PWM inverters provide the 
ac power, the usual constraints of voltage and frequency are 
nonexistent. Accordingly, the requirement for producing 
torque under high slip conditions is deleted, and thus the 
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I 

i 

need for compromise between "starting characteristics” 
versus "running efficiency" is eliminated. Accordingly, 
traction type machines are designed for maximum efficiency 
with minimal rotor resistance. 

For each torque-speed point of operation, a range of 
excitation voltage-slip frequency values may be used; 
however, only one value of voltage and slip frequency 
produces optimal efficiency. With modern inverter 

technology it is possible to control both voltage and 
frequency such that optimal system efficiencies prevail over 
the entire torque speed plane. From experience it has been 
found that the optimal solution corresponds to maintenance 
of a constant slip frequency with the excitation voltage per 
Hertz varying roughly as the square root of the torque. 

Since selection of the inverter control parameters is 
flexible, it follows that the appropriate optimization 
sequence is : 

1. Design optimized motor. 

2. Characterize this motor. 

3. Design appropriate inverter control algorithm. 

The above philosophy is embodied in the attached 

specification . 

With traction systems, cost, weight, and efficiency 
parameters are of key importance. In order to meaningfully 
optimize these parameters, relationships must be established 
between these parameters. From JPL studies it has been 
found that for passenger vehicles the "cost of weight" is 
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approximately $10/kg and the value of efficiency is 

approximately $120/%. Accordingly, for both the motor and 
the inverter, the total cost, C^, to be minimized is: 

C 1 - C + 10W - 120E (1) 

where C is the retail dollar cost 
W is the weight in kg 
E is the efficiency in % 

Since efficiency is primarily important in terms of range 
and since maximum range is under freeway conditions, E is 
defined under the speed and power conditions associated with 
freeway driving (55 mph , 12 kW shaft output). 

A totally enclosed convective cooled design is 
recommended since: 

1. Expected heat loss is low as motor efficiency 
will be in the range of 92% to 95%. 

2. Duration of operation is limited by the battery 
storage capability. (With 150 km of driving, 
heat produced by motor will be less than 1.2 
kWh. ) 

3. With open design, dirt particles may cause 
stator winding damage. 

III. DISCUSSION OF SPECIFICATIONS 
1. Voltage 

The nominal battery voltage is 240 V dc. This 
corresponds to a maximum inverter output voltage of 165 V 
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rms (sine modulation) and 185 V rms (six-step mode). The 
inverter will control to lower voltages in accordance with 
the control algorithm. 

2. Power Output 

Since the bidder shall determine actual speeds and 
torques, specifications are presented in terms of shaft 
output power and fraction of maximum speed. (The gear 
reduction ratio and inverter control constants will be 
determined after establishing the motor.) 

2a. Power at Maximum Speed 

This specification establishes one hour thermal 
and torque breakdown capability of handling the 
rated power at full speed. Since efficiency and 
breakdown torque both improve with decreasing 
speed, it follows that this power level can 
be maintained for all speeds above base speed. 

2b. Power at 10% Maximum Speed 

This specification establishes the torque 
breakdown capability below base speed and 
corresponds to vehicle acceleration. 

2c. Power at 80% Maximum Speed 

This specification establishes the torque 
breakdown capability at 0.8 N max (55 mph) and 
corresponds to vehicle acceleration at this speed. 
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3. Maximum Speed 

As maximum speed is increased, cost and weight 
decrease. Accordingly, moderately high speeds are desirable 
consistent with mechanical and windage limits. With 
previous efforts, maximum speeds between 9,000 and 20,000 
rpm have been used. The actual maximum speed selected will 
be on the basis of Equation 1. 

4. Number of Poles 

As the number of poles is increased, the conductor 
losses generally drop while the magnetic losses increase. 
(The inverter switching losses are proportional to frequency 
and, therefore, proportional to pole number; this impact is, 
however, small and will be neglected.) The optimum pole 
number will be determined on the basis of Equation 1. 

5. Weighting Function (See Equation 1) 

This function relates cost and efficiency to the retail 

price of the motor. The motor is to be designed such that 

the value of this function is minimized. The efficiency is 

defined at 0.8 N „ and 12 kw shaft output. Excitation 

max 

for this point of operation will be the optimal efficiency 
voltage (probably less than 185 V rms ) . The corresponding 
slip frequency will be noted. 

6. Efficiency and Power Factor 

Efficiency and power factor are defined under "optimal 
conditions" where the excitation voltage is such that the 



resulting slip frequency matches that associated with 
optimal efficiency at 0.8 N jnax , 12 kw. Since excitation 
voltage is limited to 185 V rms, slip frequencies higher 
than optimal will prevail for some combinations of high 
torque and high speed. Efficiency and power factor are to 
be specified for at least five torque and five rpm values. 



IV. EXPECTED RANGES OP PARAMETER VALUES 


Based on JPL studies, the following ranges of 
parameters are expected for an optimized machine: 

1. Weight 65 to 95 kg 

2. Frame Size 10" to 12" 

3. Pole Number 2 or 4 

4. Maximum Speed (N„) 8,000 to 16,000 rpm 

RlaX 

5. Base Speed 0.35 to 0.45 N 

max 

6. Optimal Slip Frequency 1 to2 Hz 

7. Efficiency at 0.8 N 

max 

12 kW 92% to 95% 

8. Power Factor for 7. 85% to 90% 
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1 . 


2 . 


3. 


4. 


5. 


Induction Motor Specification 
Input line voltage, maximum rms: 185 V 

(will be varied as required for speed and torque 
control) 

Power output: 

1 hr. (20% duty) at maximum speed (N ) 30 kW 

max 

3 min. (30% duty) at 0.8 N 60 kW 

max 

3 min. (30% duty) at 0.1 N 16 kW 

max 

Speed, N „ . between 8,000-16,000 rpm 
max 

(to be selected by contractor) 

Number of poles: two or four 
(to be selected by contractor) 

Weighting function F to be minimum: 

F = cost + 10 (weight) - 120 (efficiency ) 
where 


cost is retail cost in $ 
weight is in kg 

efficiency is in % measured at 0.8 N , 12 kW 
shaft power with optimal voltage excitation 

6. Construction: totally enclosed. 

7. Cooling: convection air, maximum ambient 40 C. 

8. Life: 4000 hours at full load. 


In presenting the quotation the bidder shall: 
Describe the selected design. 

Specify the size and weight. 


Specify the torque breakdown profile over the whole 
range of speed (either graphically or in tabular 
form) . 

Specify efficiency and power factor versus speed 
and torque for optimal excitation. 


T. W. Macie 
9 / 23/83 



APPENDIX B 


ALTERNATIVE INVERTER/ INDUCT ION MOTOR CONFIGURATIONS 








APPENDIX C 


IDEALIZED COMPARISONS OF TWO- PHASE 
AND THREE-PHASE INVERTERS 
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The two-phase inverter circuit shown in Fig. 2 and the three- 
phase arrangement in Fig. B-3 are used for these comparisons. 

Assuming Equal DC Power to a Resistive Load 
Two-phase 

For this situation/ both the load current and load voltage 
are assumed to be square waves as shown in Fig. Cl. This would 
be the case if the motor center- tapped windings in Fig. 2 were 
the primary windings of two ideal transformers with a purely re- 
sistive load on each secondary. 

Let 

4 rms volts/phase on each half of center-tapped 
winding 

- V DC 

and 

I 2^ 4 rms current/phase through the resistive load 

Assuming that the ideal transformers each have unity turns ratio 
from half of the primary winding to the secondary, 

I 2 4 > J DC^ 2 

The current I DC is the dc source current. Since the two pairs of 
switches in Fig. 2 are operated with 90° phase displacement, the 
total dc source current is always equal to twice the current in 
one phase. Basically, I DC divides equally at all instants, with 
I dc /2 flowing through one half of each transformer primary. 
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Thus 


V T,P 4 switch peak voltage 


2V, 


DC 


and 


I T ,p A switch peak current 


= W 2 


and let 


TR~ A total switch rating V_ -I™ for two-phase 

fc(p 1 / XT X. f XT 

(4) (2V dc ) (I dc /2) 

4v dc x dc 

Three-phase 

In this case, thecircuit of Fig. B-3 is considered, except 
with the motor winding replaced with a wye connected balanced set 
of load resistors. With a current I DC flowing from the dc source 
in Fig. B-3 the line-to-neutral load voltage and current are as 
shown in Fig. C-l. In this case 


v = v 
T,P DC 


X T,P “ I DC 


and let 


tr 3< £ 4 total switch rating V T , p I T/P for three-phase 


- (6) (V DC ) (I DC ) 
= 6V DC I DC 
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Therefore, the total switch rating for two-phase is 2/3 of that 
for three-phase. 

Note: This can be checked in another very simple manner 

as follows. If the power delivered to the resistive 
loads is the same for the two-phase and the three-phase 
circuits, then the power taken from the dc supply must 
be the same in each case. 


P = V I 
DC DC DC 

Also 

I T, P, 24> = I DC /2 

and 

X T, P, 3<}> = I DC 
Therefore 


TR 2sJ * 4 ^ * 2V DC* * T DC/2* 4V DC 1 


DC 


TR 34> ^ ^ V DC^ ^DC^ 6V DC I DC 


and 


TR 2(j) 2 ^ 3 TR 34> 

Assuming Equal Fundamental Power to a Resistive Load 
Two-phase 
Let 

V 1 A amplitude of fundamental voltage/phase 
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•%* I 


v 2( j ) sin cot d (cot) 


sin cot d (cot) 


4V, 


DC 


and similarly for the current 


41 


l,2<j> 


T,P,24> 


Then 


P, _ . A total fundamental power to load 
1,29 *= 


= (2) ,!!£=, 

■n/2 tt/2 


= v 1 

TT 2 DC T,P,2* 

Three-phase 

Let 

V, ,, A amplitude of fundamental voltage/phase 
r 2 tt 


_ 1 

“ ?r J v 34> 


sin cot d (cot) 


2V 


DC 


ir/3 2 tt/3 

j ^ sin cot d (cot) + j j sin cot c 
L 0 tt/3 


+ | j sin cot d (cot) 
2V3 


2V 

- DC 


(cot) 
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and for the current 


2tt 


I, ^ 1 
i, 3 * - ? 


21 


■3<{> 


sm 

cot d 

V 3 


f 

sin 


' + 


31 


PT, 34> 


Jjsin cot d (cot) 
2tt/3 

T,P.3<t> 


2tt/3 


tt/3 


Then 


>, A total fundamental load power 
1 1 3$ = 


=. (3) (^S£) t^3t) 
tt/2 tt/2 


9V DC 3 T,P, 34> 
2 


For equal load power 
P l,2c p = P l,3<p 

16 = 9V dAp ,34> 

*2 V DC 1 t,P,2<D ^2 

or 


[ = — t 

T,P,2<{> 16 ' s r,P,34> 
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Finally, 


the total switch ratings TR 


2<J> 


and TR 


34 


TR 2<t> = (4) (2V dc ) (I T< p i2t ) 


tr 34 = (6) <V DC> (I T,P,34> 


Then 


TR 2<f> (4) (2V DC J ( 16 I T,P,34> > 


= i V 1 

2 DC T,P, 3 


and 

TR 3<t> = 6V DC 


Thus 


TR 


2 * 



are 
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Two-phase load 
voltaqe and 
current 


Three-nhase line-to-neutral 
load voltaqe and current 




V 2 



Fig. C-l Load voltages and currents for Fig. 2 and Fig. B-3 respectively 
with motors replaced by resistive load 
Note; Two ideal transformers are used in Fig. 2, with unity 
turns ratio from each half of their primaries to their 
secondaries connected to resistive loads. A wye- 
connected set of balanced load resistors replaces the 
motor in Fig. B-3, 
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APPENDIX D 


DETAILED DESIGN OF TWO-PHASE BIFILAR INDUCTION MACHINE 
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1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 

11 . 

12 . 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20 . 
21 . 
22 . 
23. 


MAIN DATA OF UMC TWO-PHASE BIFILAR INDUCTION MOTOR 


Output 

Number of phases 
Frequency 

Equivalent voltage of battery 
Phase current 
Efficiency 
Power factor 

Number of poles (p a pole pairs) 
External diameter of stator 
Internal diameter of stator 
Length of stack 
Length of air gap 
Number of slots of stator 
Number of slots of rotor 
Turns per coil 
Number of coils 
Series turns per phase 
Number of branches 
Resistance of stator 
Resistance of rotor 
Leakage inductance of stator 
Leakage inductance of rotor 
Magnetizing inductance 


P m = 30 kW 
m 

m a 2 
f = 80 Hz 
V D = 189 V 

! ph = 110 * 2 A 

n = 0.91 


cos $ = 0.88 

2p = 2 

D Q = 29.21 cm (11.5 in) 
D = 14.6 cm (5.75 in) 

L = 16 cm (6.3 in) 


l g = 0.03 cm 
Zj = 24 
Z 2 = 20 


n c I 2 

T ph = 33 


V 2 

R s = 0.0274ft 


(0.0118 in) 


R^ = 0.0213ft 

L )s = 2.5 X 10' 4 H 

L, _ = 2.7 X 10‘ 4 H 
1 r 

L = 13.86 X 10’ 3 H 
m 





75 


UMC TWO- PHASE BIFILAR MOTOR 

1. Given data 

P m 3 30 kW 

m 3 2 (number of phases) 
f 3 80 Hz (basic frequency) 

P 3 1 (pair of poles) 

V D = 189 V (equivalent base battery voltage) 



E lm * £ V D * ? X 189 V = 241 V 
E l(™s) “ X '89 V = 170 V 

2. External diameter of stator D Q 

choose 0 Q 3 29.21 cm = 11.5 in 

3. Internal diameter of stator D 

let 0 3 14.6 cm 3 5.75 in 

2P<|> 

4. Specific magnetic loading S’ 3 

_ 2 
assumed B 3 0.45 wb/m 

_ 2 
B 3 0.3 —0.6 wb/m for common motor 
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5. Slots of stator Z 1 = 24 

6. Slots of rotor Z 2 

Z 2 = z i + 4P = 20 

7. Assumed Efficiency n = 0.91 

Power factor cos <j> = 0.88 

8. The kVA rating S 

p 

S = EL_ = 30 = 37>46 kVA 

cos 4> 0.91 x 0.88 

9. The motor constant G 


G = -J- = 11 kw B X ac X 10" 3 = 11 X 0.9029 X 0.45 X 30,000 X 10‘ 3 = 134 
D^Ln 

where L is length of stack 
kw is winding factor 
ac is specific electric loading 
n is synchronous speed (rps) 
kw = kp x kd kp = 1 (single layer winding) 



ac = 2m T ph I ph / D = 5000 45,000 Amp-cond/m 

Let ac = 30,000 
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in n^i - _i - 37.46 _ - v ln -3 3 

10 * D L " Gn " 134 x 80 “ 3,4944 x 10 m 


. _ ^r- _ 3.4944 X 10' 3 

L - —R 

U 0.146* 


= 0.164 m 


Let L = 0.16m = 16cm = 6.3in 

11. The flux per pole *m 


s B x || X L = 0.45 X T JL. 0 . J 46 x 0.16 = 0.0165 wb 

12. The turns per phase (in series) T ph 

t - E 1 - 170 „ , 

'ph “ 4.44 kw f Z " 4.44 x 0.9029 x 80 x 0.0165 " 1 

r m 


choose T ph = 33 


13. The real flux per pole $ 


m 


170 


^ni = 4.44 kw f T“ = 4.44 x 0.9029 x 80 x 33 = 0,016 wb 


ph 


14. The number of coils n. 


n £ * ^ s ~ 12 (single layer) 
The number of branches a p 

15. The coils per phase in series n p ^ 

n c 12 , 

n ph “ ap x m ” 2 x 2 “ 3 
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16. The turns of coi 1 Z sl 


7 S , ph. = H S 11 

51 "ph 3 


17. B 


4» 


m 


0.016 


irDL 

W 


xO. 146x0. 16 
2 


= 0.436 wb/m 


18. Phase current I 


ph 


t = S x 10 3 
x ph ~ m x Ej 


37.46 x 10 3 
2 x 170 


= 110.2 A 


ac = 2m Tp^ x Ip^/irD = 2 x 2 x 33 x 110.2 /tt x 0.146 = 31714 Amp-cond/m 

19. The density of current of stator winding 5^ 

2 

choose = 5 A/mm 

20. The section of coil S c 


S 


c 



110.2 
2 x 5 


11.02 mm 2 


choose 3 x #13 wire and 2 x #15 wire 

#13 wire d = 1.83 mm A^ = 2.63 mm 2 

insulated d' = 1.83 + 0.18 = 2.01 mm 

2 

#15 wire d = 1.45 mm A^ = 1.65 mm 
d' = 1.63 run 

S c = 3 x 2.63 + 2 x 1.65 = 11.19 irm 2 

21. The number of wires in parallel np 
n p = 3 + 2 = 5 



79 


22. The total number of wires in a slot n g 
n $ = 2 x Op x Z s j = 2 x 5 x 11 = 110 

The first number 2 on the right side represents the bi filar winding. 

23. The real density of current 


6 1 = 


_1e!l = _ 

as I 
P c 


TTTT9 ' 4 -« A/- 2 


24. Length of mean conductor 
L mc = L + 1.1Y + 3 (cm) 

= 16 + 1.1Y + 3 = 16 + 1.1 X 22.93 + 3 = 44.2 cm 


where Y = = - * 2 14,6 = 22.93 cm 

We use concentric coils: 1-8 

1-10 coil spans 

1-12 

25. Length per phase L ph 

L ph = T ph x 2 x L mc = 33 x 2 x 0.442 = 29.17 m 

26. Resistance per phase R $ (75°C) 

R s = P ?5 x L ph /S c x a p = 0.021 x 29.17/11.19 x 2 = 0.0274fi 

I 

27. Reduced rotor current I 2 

v z a I ph cos <j> = 110.2 x 0.88 = 97 A 

28. The tooth pitch of rotor t $ 2 


. " D 2 . it x 14.54 . „ 

T s2 ” 25 m 

where diameter of rotor D, = D - 2x1 

2 9 

lg = 0.3 mm (length of air gap) 


14.6 - 2 x 0.03 = 14.54 cm 
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29. Number of slots per pole S 2 (rotor) 


S =-1 = 20 = 1Q 
b 2 2p 2 IU 


30. The current of the bar I b 

r _ 2m k w T ph T ' _ 2 x 2 x 0.9029 x 33 .. n , _ 

i. - 7 b ~ x y/ - 


31. The current of the end ring I, 


I K S, 

I = = 578 x 10/ir = 1840 A 

C IT 


32. Casting aluminum alloy for rotor cage 


PA1 (75°C) = 0.04Q mnr/m 

(The value can be in negotiation with manufacturer) 
Let 6 2b = 3.2 A/mm 2 


34. The section of each bar a. 


a _ b _ 578 _ ,o n c _,2 
” f. “on** I OU • o m 
2 <$ 2b 3.2 


35. The resistance per bar A b 


A b = P A1 x L/a 2 = 0.04 x 0.16/180.6 


3.54 x 10" 5 n 


36. The section of end ring a 2c 


Let 6 2c =2.8 A/mnr 

a _ c _ 1840 _ ,, 7 2 

a 2c " ^Zc " 2T8 - ' 657 "" 



37. The resistance of ring r Q 


r c = p A1 x it x 0 r /a 2( . = 0.04 x n x 0.126/657 = 2.41 x 10"^ ft 

where D r = 0.126 m (the mean diameter of ring) 

38. The reduced resistance of rotor r 2 

, . Z 2 (r b + 2r c S 2 /7T2z 2^b = 20 (3.54 + 2 x 2.41 x 10 2 /tt 2 

2 m I£ 2 2 x 97 2 

= 0. 02124ft = 0.0213ft 

39. Carter's factor k c 


k c = k cl • k c2 


T S1 


cl b~ 

- b n + 21 n x 1 [1 + — ] 

S1 0 9 n 2 xl g 


19.1 


19.1 - 3.5 ♦ 2 x 0.3 x ! n [l + 


1.14 


where t, 


SI 


ttD tt x 146 
Z, " 24 


= 19.1 mm 


(tooth pitch of the stator) 
b Q = 3.5 mm (width of slot opening - refer to Fig. 1) 


S2 


c2 


- b o2 2 ] g x ] n [1 + 2T 3 


2.284 X 10 

22.84 - 1.5 + 2 x 0.3 x l p [1 + 2 3^3] 


1 .034 


k c = k cl . k c2 = 1.14 x 1.034 = 1.18 
Effective length of air gap l g 




I O il 19:2x3] . * _ t ^ s ( 14 6 : |xO,3-6K . 9>2 = 12 .7 mm 


17 mm 

(D-21g-2x3-2xd 2 )iT ^ (146-2x0.3-6-2x17) tt 

Z- t Z2 


7.356 mm 


Fig. D-2 
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l g = k c l g = 1.18 x 0.3 = 0.354 mm 

40. The area of air gap A g 

A g L = 2 ffi- . 146 . x 0.16 = 0.0367 m 2 

41. Modified density of flux of air gap B 

B g = 1.36 x * m /A g = 1.36 x 0.016/0/0367 = 0.593 wb/m 2 

where 1.36 is for considering saturation in both teeth of stator and rotor 

42. Area of teeth per pole of stator A^ 

A/i = t 2 j x L x 0.92 x Sj = 0.78 x 16 x 0.92 x 12 

= 138 cm 2 = 0.0138 m 2 
2 1 24 

where = ^- = -g = 12 

43. Modified density of flux in teeth 

B Z1 = 1.36 x<J> m /A zl = 1.36 x 0.016/0.0138 = 1.577 wb/m 2 

44. Length of tooth 1 Z1 

^ Z1 = + kg = 31 + 2 = 33 mm = 0.033 m 

45. Height of core h cl (stator} 


D 0 - [D + 2 ( h i + h 2 + h 3 + h 4 )j 

h cl = ? 


= 29.21 - ri4.6 + 2 (3.1 + 0.2 + 0.2 + 0.1)1 B 


3.705 cm 


= 0.03705 m 
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46. Area of core A cl (stator) 

A cl = h cl x L x °* 92 = °- 03705 x °* 16 x °- 92 = 0.005454m 2 

47. Maximum density of flux through the core 

O _ ^m _ 0.016 _ , .,.0 u/K /m 2 

B cl ' ' 2 x 0.005454 ' 1,4668 wb/m 

48. Length of core per half pole l c -| 

, c] = = < 29 - 21 j ■ 3 - - - 705 )" = 20.03 on 

I 

49. Effective length of core 1 -j 


1 1 = d cl = 0.38 x 20.03 = 7.6 cm = 0.076 m 

50. 5 = 0.38 considering the densities of flux along the length of 

core are not equal 

51. The fill coefficient of slotk^ 


n^d 


_ 2 x 11 x (3 x 2 . 01 2 + 2 x 1.63 2 ) . „ 

r w s +w b x ( h 1 + h 2 ) x {31 + 2) 

~T~ 

kf = 0.702 < 0.75 (This means that bifilar winding can be inserted 

in slots without difficulty.) 

52. Area of rotor teeth per pole A^ 


A Z2 = b Z2 x b x 8,92 X S 2 = 0.92 x 16 x 0.92 x 10 = 135.4 cm = 0.01354 m 
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53. 


54. 


Modified density of flux in rotor teeth B 



1 ' 36 = 1.36 x 0.016 
A z2 0.01354 


Length of rotor tooth 1^ 


= 1 .607 wb/m 


Z2 

2 


1 Z2 - d 2 3 17 rmr. = 0.017 m 


55. The height of rotor core h c2 

2 x 2 - d 


c2 


D - 2 1 - 2 d 9 

g 2 




146 - 2 x 0.3 - 2 x 17 - 4 


— + \ x 45 


= 38.7 mm 

where d Q = 45 mm (d Q is diameter of the shaft) 


56. The maximum density of flux through the rotor core B c2 


57. 


m 


0.016 


c2 2A 


c2 


2 x 0.0057 


= 1 .404 wb/nf 


Area of rotor core A 


c2 


A c2 ■ h c2 x L x 0.92 = 0.387 x 0.16 x 0.92 = 0.0057m 2 

According to the B-H curve on Fig.D-3(Sec. 4-264, Fig. 4-31, Standard Hand - 
book for Electrical Engineers , Tenth Ed., Donald G. Fink and John M. Carroll, 
McGraw-Hill), we get 

B Z1 = 1.577 wb/m 2 = 15770 gauss 

H Z1 = 32 oe * 32 x 79.577 = 2546 A/m 

B , = 1.4668 wb/m 2 = 14668 gauss 
cl 

H = 17 oe = 17 x 79.577 = 1352 A/m 
cl 

B z2 = 1.607 wb/m 2 = 16070 gauss 
H z2 = 40 oe = 40 x 79.577 = 3183 A/m 
B c2 = 1.404 wb/m 2 = 14040 gauss 
H o = 13 oe = 13 x 79.577 = 1034 A/m 

u C2 
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MAGNETIC MATERIALS 

H 



Magnetising force^H, oersteds 


Fig. D_3 Normal-induction curve 
and hysteresis loop of electrical sheet, 
No. 22 gage. 
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58. Magnetizing force AT 

AT = AT g + AT Z1 + AT cl + AT Z2 + AT c2 

59. Magnetizing force of air gap ATg 

ATg =8x10^x1gXB g =8x 10 5 x 0.000354 x 0.593 = 168 A 

60. Magnetizing force of stator tooth AT^ 


AT Z1 = H Z1 x ] Z1 = 2546 x 0,033 = 84 A 


61. Magnetizing force of stator core AT c ^ 


AT d = H cl x 1^ = 1352 x 0.076 = 102.7 A 


62. Magnetizing force of rotor tooth AT Z2 


AT Z2 = H Z2 x ] Z2 = 3183 x 0,017 = 54 A 


63. Magnetizing force of rotor core AT c2 


AT c2 = H c2 ] c2 = 1034 x 0.0152 = 16 A 


where 1^ 2 = ? l c2 = 0.38 x 4 = 1.52 cm = 0.0152 m 


0 - 2 1 - 2 d, - 2 h', 

1 „ = 9 2 c2 

c2 4p 


- 4 cm 


AT = 168 + 84 + 102.7 + 54 + 16 = 424.7 A 


64. Magnetizing current I 


or 


I or - P X 2AT/ 1.17 x kw x T „ - 2 X 424.7/, x 


.9029 x 33 


24.4 A 



65. Magnetizing reactance X 


X- ~ ■ Jo - = 6.967ft 


m 


or 


66. Magnetizing inductance L 


m 


I = — E = - § - 967 = 13.86 x TO" 3 H 

2rrf 2 xttx80 


67. The specific slot permeance of stator X $ g 


'si ■ upt r + 


W, h U. h. h.tl 

£_ In f-il + — + -i 

W ' ' W -b ln [ bj b n W. 
s s o o o b 


B 2 . 6_ _ lng 


(1-8)(1-B 2 ) 2 


— In [ ■ + £ In ( 1 ) + 

.1 M2.1 ’ 12.1-3.5 1 3.5 


12.62-12 


12. lx J_ 
3.5 


32 

21 


-2 CL6 
4 


0.6 - ^ - In 0.6 - £ 


( 1 - 0 . 6 ) ( 1 - 0 . 6 2 ) 2 


= 0.162 + 0.288 + 0.286 + 0.822 = 1 


where 


u _ [146 + 2 (h. + h. + h ? )) ir 

w 2 ' ^ - - 7.8 = 12. 


62 mm 


W $ = 12.1 mm 
h 2 = 2 mm 
hg = 2 mm 
h 4 = 1 mm 
b Q = 3.5 ram 
W 

8 = ^ = 0.6 (refer to Fig. 1) 


•Hu>, 
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68. The specific slot permeance of rotor X g2 


h' W c9 , d, 

X s2 = W s2 -f.5 ln ( 1T5 ) + TT5 + 


b 2 -£- me - f 

(l-0)(l-6 2 ) 2 


12.7-1.5 


# 12. 7 > . J_ + _ 

In ( 1.5 } 1.5 7. 


17 


356 


1.726 2 - 1^- - In 1.726 - | 
(l-1.726)(l-1.726 2 ) 2 


where W g2 = 12.7 mm 
W^ 2 = 7.356 mm 


d 2 = 17 mm 
W s2 

e - ^ = 1-726 
w b2 


69. The reduced value A 1 


s2 


v .*h x - 0.9029 2 x 24 x , .. . , 
* s2 Zo s2 20 A '* 4y 


70. The slot leakage reactance of stator x sl 
r 2 

_ ° r m L r 
c sl 


Qf I T 4 * O 

x c1 = ■ m - 7 - £h x = 8 x 80 - x 2 x 6 * — - x 33 -. x 1.558 = 0.057 Q 

** in' *7 SI in' w 


10' Z 


1 


10' x 24 


where L is in inches 

71. The slot leakage reactance of rotor x s2 

2 


8fm L T, 


x s2 


— 7 — 
10 7 Z. 




x' 


s2 


_ 8 x 80 x 2 x 6.3 x 33‘ 

3 


10' x 24 


X 


= 1.49 


x 1.46 = 0.535 B 
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72. The end leakage reactance x. 


X* = 


4fm D T £h_ = 4 x 80 x_2 x_5.75 x 33 2 a Q 1 


e 10 7 (2p) 2 


10 7 ( 2) 2 


where D = 14.6 cm = 5.75 in 


73. Differential leakage reactance (sum of the zigzag and belt leakage 
reactance) 

,5.. , 6K C1 - V 5g2 . Z^ 1 , . | x 6M1 „ ( 6 X . LH - l^f l ajy 


x d ~ (-T x m ^ 5 S 1 2 


+ HV#- = °’ 0516 


a = 0.5 is the angle of rotor slot skew, as a fraction of one stator slot 
pitch 

74. The leakage reactance of stator xj 


Xi * x sl ♦ (x e ♦ X d > * 0.057 + 55-1^5 (0.1 +0.0516) = 0.126 


75. The reduced leakage reactance of rotor x * 2 


x 2 = x s2 + Z~ TT~ ^ x e + x d^ = °- 0535 + (0.1 + 0.0516) = 0.136 n 


76. L 


. X 1 . 0.126 


Is 2rrf 2 x ir x 80 


2.5 x 10' H H 


, _ X 2 _ 0.136 _ o 7 „ ln -4 .. 

hr “ srr 1 xTTso ■ 2 - 7 x 10 H 



f = 20 Hz V D = 60 V 

f = 80 Hz V Q = 189 V 

f = 160 Hz V D = 240 V 

R $ = 0.0274 

R $ = 0.0274 

R s = 0.0274 

RJ. = 0.0213 

R' r = 0.0213 

R^ = 0.0213 

L ]r = 2.41 x 10~ 4 H 

L ]r = 2.7 x 10‘ 4 H 

L ]r = 3.23 x 10‘ 4 H 

L 1$ = 2.256 x 10" 4 H 

L ]s = 2.5 x 10' 4 H 

L ]s = 2.94 x 10' 4 H 

L m =6.51 x 10‘ 3 H 

L = 13.86 x 10~ 3 H 
m 

L = 26.8 x 10" 3 H 
m 

I n „ = 66 A 
or 

(rms) 

I or = 24.4 A 
(rms) 

I 0 , = 8 A 
or 

(rms) 

I_„ = /? x 66 
or 

(peak) 

• I = n x 24 
or 

(peak) 

I _ = JZ x 8 
or 

(peak) 

= 93 A 

= 34 A 

- 11 A 






Modified Parameters of the Motor Operating at 


Modified Parameters of the Motor Operating at 


f = 20 Hz V D = 60 V 


f = 160 Hz V D = 240 V 


(1) ^ V n = — - — X 60 = 54 V 

1(ras) /Jx, 0 /!!,. 


'1 (rms) 


/? X TT 


x 240 = 216 V 


mi * b l 54 

" 4.44 kw f TT " 4.44 x 0.9029 x 20 x 33 
ph 

= 0.0204 

Neglect the change of resistance, as the 
operation frequency is changed. 


= 4.44 x 0.9029 x 160 x 33 = °' 0102 


(26) R $ = 0.0274 0 


R $ = 0.0274 n 


(38) R' r = 0.0213 0 


R' = 0.0213 a 
r 


tAi\ n _ , o/r , /* _ 1.36 x 0.0204 
(41) Bg = 1.36 - — 0357 


B g * 1,3 |.0 3g 7 --' =0 - 378 


0.756 


(43) B 71 = 1.36 A/A 


m' "Z1 


1.36 x 0.0204 
0.01 38 


R = * 0.0102 _ i QQg 

B Z1 0.0138 '- uua 


= 2.01 


0.0204 


(47) B cl - 4> m /2A cl = 2 x q!o05454 = 1,87 


R = 0.01 02 _ q goc 

B cl 2 x 0.005454 



(53) 



1,36t> m _ 1.36 x 0.0204 
A Z2 0.01354 


2.049 


(56) B c2 


_ 0.0204 

2 x A c2 2 x 0.0057 


1.79 


(57) B Z1 = 2.01 

B Z1 = 0.9 x B Z1 = 1.809 


(When teeth are over saturated, part of 
flux passes through slots.) 

H Z1 = 130 x 79.577 = 10345 

B cl = 1.87 

B^ = 0.9 x 1.87 = 1.683 

(Part of flux spreads to housing.) 

H cl = 55 x 79.577 = 4376 
B Z2 = 2.049 


1.36 x 0.0102 
Z2 ~ 0.01354 


1 .0245 



0.0102 
2 x 0.0057 


0.895 


B Z1 = 1.005 


H Z1 = 3 x 79.577 = 239 
B cl = 0.935 


H cl = 2.6 x 79.577 = 207 
B Z2 = 1.0245 


Bi, 2 = 0.9 x 13 Z2 = 1.84 


to 

CO 



H Z2 = 160 x 79.577 = 12732 



H z2 = 3.2 x 79.577 = 255 
B c2 =0.895 

H c2 = 2.5 x 79.577 = 199 

ATg = 8 x 10 5 x 0.000354 x 0.378 
= 107 

AT n = 239 x 0.033 = 8 

AT cl = 207 x 0.076 = 16 

AT Z2 = 255 x 0.017 = 5 

AT c2 = 199 x 0.0152 = 3 

AT = 107 + 8 + 16 + 5 + 3 = 139 

t - 2 x 139 

or ' 1.17 x 0.9029 x 33 0 


VO 

4 * 



tn 


(66) L =^ = 


m _ 0.818 


m 2nf 2 x it x 20 


6.51 x 10 


-3 


27 


'm 2 x x 160 


= 26.8 x 10' 


The overload current is around two times larger than rating, so slot leakage inductance still can be 
considered as a constant. 

(70) X sl * 8fm L _ T ph X S1 

in'- 


10'Z, 


v _ 8 x 160 x 2 x 6.3 x 33 . C m 

A - - 7 X 1.030 

Si lo' x 24 


= 8 - x 2 P x - 2 2L 6 j. 3 x_33 , x 1.553 = 0.01425 
10' x 24 


= 0.114 


8fm L T‘ 


(71 ) X s2 = 0,1,1 L ’ph = X 
10 7 x Z, 


s2 


y _ 8 x 160 x 2 x 6.3 x 33 x 1 46 
s2 10 7 x 24 


8 x 20 x 2 x 6.3 x 33* 
1 0 7 x 24 


X 1.46 = 0.1336 


= 0.107 


do\ v - 4fm 0 T 2 . 4 x 20 x 2 x 5.75 x 33 2 

(72) X e - Pn = 


10 7 (2p) 2 


10 7 ( 2) 2 


X e 


4 x 160 x 2 x 5.75 x 33 c 

10 7 ( 2) 2 


- 0.025 


= 0.2 


6k- - 1 - Sa c ok- -1 

03 ) 


d 6m' 5S ^2 


5S 2 2 


V . 5 „ .. ,6 x 1.14 - 1 + 5x 0.5* 1 t 0.5x1 .034-1 \ 

X d “ 6 X L x ' 5 5 — ' 


5 x 12 


5 x 10‘ 


= I- x 0.818 x + Orlx b OM-l) 

b 5x1 2^" 5 x 10^ 


= 0.2 


to 

cn 


0.00605 



(74) X x = X sl + j-£- (X e + X d ) = 0.01425 

+ (0.025 + 0.00605) = 0.02835 


(75) X£ - X s2 + (X e + X d ) 

= 0.01336 + (°- 025 + 0.00605) 

= 0.0303 n 

( 7fi \ I - A - 0-02835 - 9 9ZR x in'' 

t76) L 1s " 2trf " 2 x it x 20 " 2,256 x 10 


= _^ = 0.0303 = 2 41 1(J -4 

L lr 2irf 2 x it x 20 x lu 


R s = 0.0274 n 
RJ. = 0.0213 n 
= 6.51 x 10' 3 


X 1 = 0,114 + ^ 0T?4 (°‘ 2 + 0 ‘ 2 > 
= 0.296 

X- = 0.107 + 2^ (0.2 + 0.2) 
= 0.325 



R $ = 0.0274 Q 

Rj. = 0.0213 n 

L = 26.8 x 10' 3 
ro 


UD 

o> 
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MOTOR TRANSIENT TEMPERATURE RISE 
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Peak Momentary Motor Temperature Rise 

The specifications require three motor operating points. 

1. 100 per cent speed, 30 kW output for 1 hour 

2. 10 per cent speed, 16 kW output for 3 min. 

3. 80 per cent speed, 60 kW output for 3 min. 

The first one is obviously no problem. 

'ph ’ m * 110 ' 2 * 86 ' 7 fl 

$1 = 2 x 8 f j jg = 3,87 A/ 1 ™ 2 (See Appendix D) 

The copper losses are small due to the low current, and the motor operates 
under saturation, so the core losses are also small. Although the mechanical 
losses are greatest at top speed, high speed creates very good heat dissipation 
conditions. One hour operation is expected to produce nearly steady state 
temperature rise. 

The second one is worse than the third case for temperature rise because 
at this operating point there are maximum losses and worst case heat dissipation 
conditions. 

The heat balance differential (1) can be applied to evaluate the 
transient temperature rise. 


d8 . SA _ p 

dt 0 GC " GC 
P P 

where 0 is temperature rise, °C 

2 

S is cooling surface of the hot body, m 
G is weight of the hot body, kg 
C p is specific heat of the hot body, J/kg-°C 
X is specific heat dissipation or emissivity, 

J per second perm of surface per °C difference 
between surface and ambient cooling medium 
P is heat developed, J per second or W 


( 1 ) 
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During a short time we neglect the exchange between stator and rotor (the 
heat- resistance of air gap is very large), and even neglect the heat 
dissipation. This means that all the heat developed causes temperature 
rise, and of course, this is the most serious situation. 

Then (1) becomes 

de p 
it 

9 = f GC“ dt 
P 

If we assume P and C p are constant, then 


where t * 3 minutes = 180 s 

_3 

G = m x 2 x Lpij x S c x a x 10 x g c 

= 2 x 2 x 29.17 x 11.19 x 2 x 10 -3 x 8.9 = 23.24 
2 

g c = 8.9 kg/dm (specific gravity of copper) 

The total weight G of the two-phase bi filar winding is greater because each 
single winding conducts for only a half cycle. Thus, the bifilar windings have 
an advantage for short time transient temperature rise - not significant for 
steady state. 

The specific heat C p of copper is 

C p = 0.092 Btu/lbF 
1055 

• 0.092 x Q>454 x 5 - 384.8 J/kg°C 
P is stator winding losses 
P = m x r $ x I^ s = 2 x 0.0274 x 200. 4 2 
= 2201 

(As shown in Appendix G, 1^ is 200.4 A) 


( 2 ) 


(3) 
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From (3), we get 

o - P f. - 2201 » 1 fin 

GC p 1 ' 23.24 x 384.8 x l0U 

= 44. 3°C 


We suggest the insulation class F. Specifications for the F class insulation 
are as follows:* 

The hot-spot temperature is 155°C. 

The maximum permitting temperature rise for totally enclosed 
fan-cooled motors is 105°C, based on ambient temperature 
of 40°C,3330 ft. altitude, and determined by the resistance 
method. 

If we directly apply (1) for 30% duty 


then 


de , S\ a 

dt or 6 = 


p 

P/GC r 


GC 

P 

P/GC 


SX 


SX/GC p SX/GCp 

P P - — ■ 
= SX " SX e GC p 


p Q GC 
n e p 


where X = 1/c = 1/0.04 = 25 J/sec m °C 
(assume no cooling air) 


S = Z. x L x L cn = 24 x 0.442 x 0.091 
1 me sp 

= 1.05 m 2 

L $p is the perimeter of stator slot 


(4) 


* Standard Handbook for Electrical Engineers , 11th Edition, McGraw-Hill, pp. 20-45. 
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L sp = 2(h l +h 2 ) + W b + W s 

= 2 (31 + 2) + 21 + 12.1 = 99.1 mm 
= 0.0991 m 


For Zj, L mc , hj, h 2 , W fa , and W g refer to the design sheet of the motor in 
Appendix D. 


From (4) 

1.05x25 

a _ 2201 2201 _ ' 23.24x384.8 

" 1.05 x 25 ’ 1.05 x 25 e 


= 83.85 - 83.85 e 


■0.00294t 


30% duty 3 minutes operation 
7 minutes off 


For the first operation 


t = 3 x 60 = 180 sec 


= 83.85 - 83.85 e 


-0.00294 x 180 


= 34.5°C 


During the 7 minute off interval 

0' _ 0 „ -0.00294 x t 

1 - i e 

= 34.5 e ■°- 00294 x 420 . 10"C 
The second duty cycle 

0 2 = 0 : + 10 = 44.5°C 

0' = 44.5 e -°- 00294 x 420 = 13°C 

2 



The third duty cycle 


0 3 = 9 1 + 13 = 47 - 5 ° c 

6' = 47.5 x e - 0 ' 00294 * 420 

The fourth duty cycle 

e 4 = + 13.8 = 48. 3°C 

04 = 48.3 e -°' 00294 x 420 = 
The fifth duty cycle 

0 5 = 0j + 14.05 = 48.55°C 
0- = 48.55 e -O' 00294 * 420 

The sixth duty cycle 

0 g = Q 1 + 14.12 = 48.62°C 

9' = 48.62 e ■°- 00294 x 420 

D 


= 13.8°C 


14.05°C 


14.1 2°C 


14.14°C 
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Temperature rise 



Fiq. E-l Transient Temperature Rise 


APPENDIX F 


DETAILED INVERTER DESIGN 
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F.l INTRODUCTION 

The complete system schematic diagram is shown in Fig. F-l. This is 
identical to Fig. 12 of section 4.4 in this report, but it is included again 
here for convenience. This appendix describes the detailed design of the 
two-phase inverter-bifilar induction motor system. 

F.2 POWER CIRCUIT DESIGN 

Fig. F-2 is the power circuit (same as Fig. 2). The most important 
aspect of the design of the power circuit is the determination of the rating 
requirements for the power switch devices. Analysis and simulation of the 
two-phase bifilar induction machine define the ideal switch voltage and current 
ratings to supply the power required at the specification operating points. 

The power switch ratings (section 3.4.1 ) .assuming ideal circuit conditions, 
are 

Peak Voltage V p = 480V 

Peak Current I p = 245A 

RMS Current I RMS = 200A 

Maximum 

Switching f<- = 3.4kHz at 245A peak 

Frequency 


However, it is necessary to analyze the specific circuit configuration to 
determine the precise transient and steady state voltage and current requirements 
of the switching devices. One of the primary disadvantages of the circuit in 
Fig. F-2 is the existence of greater transient voltages. These are the result 
of leakage inductance due to imperfect coupling between the bifilar windings and 
motor lead inductance which is now in the commutation path where current switching 
occurs. The lead inductance does not produce voltage transients in the three- 
phase bridge system. 
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F.3 TRANSIENT VOLTAGE SPIKE ANALYSIS 

Consider Fig. F-2 at the moment when switch SI is opened to transfer its 
current to the S2-D2 branch. If the A1 and A2 bi filar windings are perfectly 
coupled, and if there is no lead inductance in the circuit including the dc 
source, motor leads and switch-diode branches, then the current commutation 
can occur instantly with no voltage transient. Unfortunately, the bifilar 
winding coupling is not perfect and lead inductances do exist. Therefore, 
voltage transients are generated proportional to the rate of change of current 
and the equivalent leakage inductances of the bifilar windings and the battery, 
motor and switch-diode wiring inductances. Thus, it is essential to minimize 
these inductances and determine accurate estimates of their maximum values. 

Snubber Design 

Fig. F-3 shows a simplified model of a single snubber circuit including a 
lumped inductance to represent bifilar winding leakage inductances and lead 
inductances. 

The order of magnitude of the spike voltage is determined from assuming that 
all of the initial energy stored in the inductance is transferred to the capacitor, 

that is, 

^L I p 2 =lc V p 2 (1) 

where 

Ip A peak current at instant switch opens 
Vp A peak voltage on capacitor 
L A equivalent total leakage and lead inductance 
C A snubber capacitance 
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Fig. F-3 Simplified Switching Circuit 



no 


Thus 



First, it is necessary to determine the maximum value of L. This inductance 
is primarily t! ■» leakage inductance due to imperfect coupling in the bifilar 
windings. Normally, the lead inductances are small, particularly if the 
inverter is mounted close to the motor and if a high quality capacitor Is 
connected direct.y from the motor winding center-tap to the negative dc bus, 
which minimizes the effect of battery lead inductance and the inductance in 
the winding center-tap lead. 

The normal total leakage reactance for a 10-10Qhp induction motor is 
0.23 per unit (pu). Of this total, the stator leakage is about 0.1 pu. This 
stator leakage can be reduced to about 0.05 pu for hiah torque motors. Finally, 
the equivalent leakage inductance between two tightly coupled windings in one 
motor phase should be less than .005 pu. This results in an estimated maximum 
value of L equal to 15pH. However, it is estimated that the leakage inductance 
between two bifilar motor windings could be reduced an order of magnitude below 
15pH with special motor winding techniques. For the following calculations 
an L of 4pH is used. With I = 24 5A, from (2) 

V p 245 = (3) 

K L /C 

If C is chosen equal to 0.33pF, 

v p * ow * 853V < 4 > 

This is a reasonable limit for Vp. 

The snubber R is necessary to limit the peak snubber capacitor discharge 
current at the moment of turning on the power switch. If R is chosen equal 
to 10 ohms, then 



Ill 


Ip(O) - 



48A 


(5) 


This is a reasonable I p (0) to allow for the device rating. 

Voltage Spike Simulation 

The snubber design carried out in the previous section is based on an 
approximate transient voltage estimation. The precise switching behavior 
can be determined by digital simulation including all non-linear' effects, 
or by hardware experiments. Since no hardware construction was done, 
digital simulations were performed using SPICE to provide more accurate 
results than could be obtained from aoproximate analysis. 

The switching transient is short relative to changes which occur in 
the rotor circuit of the induction motor. Thus, the steady state equivalent 
circuit model is used for the motor load, i.e., ^/s where r« is rotor re- 
sistance and s is the slip. The equivalent stator circuit is derived using 
the three winding transformer input-output relations. Fig. F-4 shows the 
flux linkages of the transformer, which is equivalent to the bifilar wound 
induction motor with the assumptions that the switching transient is very 
fast and that all windings have equal turns. The #1 and #2 windings are 
the stator bifilar windings, and the #3 winding is the rotor winding, v^ , 

V£, v 3 and i^, i^, i ' 3 are the terminal voltages and currents for each winding. 
d>-| • 4>2 and $ 3 are defined as the total flux of fl , £ 2 and *3 windings, re- 
spectively. , <t>£ 2 > ^3 are the corresponding leakage flux for the * 1 , 

#2 and #3 windings, respectively. is the mutual flux between either 
stator winding and the rotor winding, and 4 >g is the mutual flux between the 
two stator windings. The following relations are obtained. 

*1 = *41 + *B + 

*2 = hz + *B + 

*3 = *£3 + *M 


( 6 ) 

(7) 

( 8 ) 
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where v. 


<ty 


1 


d$. 


'1 dt * v 2 = dt * and v 3 = dt 
Since the #1 and #2 windings are bifilar windings, they are very tightly 
coupled. Thus, all of the flux links the #3 winding as well as both the 
#1 and #2 windings. 

Defining the mutual (or magnetizing) inductance I_ M between the #1, P2 
and f?3 windings, the mutual flux is 


d(J>- 


1 1 + ( i l + i o ) 


(9) 


Defining the mutual inductance between the rl and *2 windings, the 


mutual flux <£□ is 


**B = L MB n 1 + L MB ^2 


( 10 ) 


Defining the leakage inductances , L„ 2 , and L ?3 of the *1 , *2 and 
*3 windings respectively, then ^ and ^3 are 


ai = 4i n i 

(ID 

12 = l l2 1 2 

( 12 ) 

'13 L £3 ^ 

(13) 


Combining ( 6 ) - (13) 

*1 = L fi.l n l + W’l + V + L M ^1 + ’2 + ^ 

*2 = L £2 *2 + L MB^1 + V + L M (l l + n 2 + V 

^3 = l l3 '3 + L M (l l + ^ + V 

These three equations provide the equivalent circuit shown in Fig. F-5, 

where the rotor load and stator resistance r-| are added. v-| and v 2 are 
controlled to be E, -E or 0. This equivalent circuit becomes the conven- 
tional equivalent circuit if one of the bifilar winding branches is removed. 



IY> 
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SPICE is a general purpose circuit simulation program. Circuits may 
contain resistors, capacitors, inductors, mutual inductors, independent 
voltage and current sources, four types of dependent sources, transmission 
lines, and the most common semiconductor devices: diodes, BJTs, JFETs and 

MOSFETs. 

Fig. F-6 shows the equivalent circuits used for simulation, including 
snubber circuits and base drive circuits. When one device is turned off, 
the base drive voltage source is changed linearly from 10 volts to zero for 
a 5 microsecond interval. After a 15ys delay, when the other is turned on, 
the base voltage is controlled linearly from zero to 10 volts for a 5ps 
interval. When one device is switched off, the current first flows through 
the feedback diode of the other device. Thus, the delay time in turn-on of 
the next device has no effect. The interval during which the base voltage 
source is linearly controlled is changed from 5ys to 2ys, but the transient 
voltage changes no more than 0.4? in its peak value. The lumped leakage 
inductances LLS1 and LLS2, shown in Fig. F-6, are chosen to be 4yH, while 
the sum of LLS1 and LMB is equal to the stator leakage inductance. LM and 
LLR are equal to the magnetizing and rotor leakage inductances respectively. 
RSI and RS2 are the stator winding resistances, and RR2 is calculated from 
r 2 /s where r 2 Is the rotor resistance and s is the slip at the steady state 
operating point. The series resistances Rll and R22 help to suppress the 
transient voltage. By trial and error, 3.9 ohms are selected. A peak tran- 
sient voltage of 943 volts across the transistors is obtained for the follow- 
ing operating conditions and system parameters: 

10% speed, 16 kW output 
Percent slip a 8.5% 

Stator current = 245A 


/ 


/ 
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Magnetizing current (through LM) = 93A 

Rotor current = 142A 

LM = 6.51mH; LLR = 241yH 

LLS1 = LLS2 = 4yH; LMB = 221. 6yH 

RSI = RS2 = 0.0274ft; RR2 = 0.0213/0.085 = 0.2510 

R1 = R2 = 10ft; Rll = R22 = 3.90 

Cl = C2 = 0.33yF 

The transient waveforms of the collector-emitter voltage and bifilar 
winding current (i-j in Fig. F-6) are plotted in Fig. F-7, and the corres- 
ponding simulation source program is listed in Fig. F-8. Without Rll and R22, 
the transient voltage is about 1160 volts. 

Next, the equivalent leakage inductances LLS1 and LLS2 are changed from 
zero to lOyH and the peak transient voltage is plotted in Fig. F-9. The 
operating condition of 80?» speed, 60kW output is also simulated, and the 
transient voltage variation is shown in Fiq. F-9 using the following parameters: 
Stator current = 245A 
Magnetizing current = 20A 
Rotor current = 225A 
LM = 26.8mH; LLS1 + LMB = 294yH 
RR2 = 0.0213/. 02 = 1 .07(2(2% slip) 

Other parameters same as for 16kW, 10^ speed condition 
Fig. F-9 indicates that the spike voltaqe is somewhat less than 1000V if the 
equivalent leakage inductance is 4yH. 


/ 



Fig. F-7 Voltage spike simulation : 16kw, 10%speed, Ip=245A, 

C1=C2=0.33 pF, LLS1=LLS2=4 pH, Rll=R22=3.9n , Rl=R2=10fl in Fig.F-6 


o 



0.0000-05 

8. 1000- 05 
O.200D-05 

6.3000- 05 
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0.5000-05 

0.6000-05 

0.7000-05 

0.0000-05 

0.9000-05 

9.0000-05 

9.1000- 05 
9.200D-05 

9.3000- 05 
9.6000-05 
9.5000-05 
9.6000-05 
9. 7000-05 

9.0000- 05 
9.9000-05 

1.0000- 06 
1.0100-06 
1.0200-06 
1.0300-06 
1.0*00-06 
1.0500-06 
1.0600-06 
1. 0700-06 
1.030D-06 
1.0900-06 

1. 1000- 06 
1.1100-06 
1 . 12C0-06 
1.1300-06 
1.1600-06 
1.1500-06 
1.1600-06 
1.1700-06 

1.1000-06 
1.1900-06 
1.2000-06 
1.2100-06 
1.2200-06 
1.2300-06 
1.2600-06 
1.2500-06 
1 . 260D-06 
1.2700-06 
1.2000-06 
1.2900-06 

1.3000- 06 


9.6260-02 
9.6330-02 
9.63°0-02 
9 6660-02 
9.6520-02 
9.6590-02 
9.6660-02 
9.6720-02 
9.6790-02 
9.6050-02 
9.6920-02 
9.6920-02 
9.5050-02 
9.5110-02 
9.5120-02 
9.5260-02 
9.5310-02 
9.5370-02 
9.5660-02 
9.5500-02 
9.5570-02 
1.1990-01 
3.6160*02 
9.6250*02 
8.5790*02 
5.2650*02 
6.6660*02 
6.6660*02 
6.5550*02 
6.6220*02 
6.6600*02 
6.6370*02 
6.7060*02 
6.7160*02 
6.7220*02 
6.73.0*02 
6.73 r I)i02 
6.7620*02 
6.7660*02 
6.7650*02 
6 . 766D*02 
6.7670*02 
6.7670*02 
6.7620*02 
6 . 7600*02 
6.7600*02 
6.7600*02 
6.7600*02 
6.7620*02 
6.7670*02 
6.7670*02 
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***** 1/28/85********* SPICE 2E.2 (26SEP7S) 


15445436***** 


JPL PROJECT VOLTAGE SPIKE 12/18/84 


INPUT LISTING 


TEMPERATURE » 27.000 DEG C 


«w 


*** BIFILAR MINDING 3 TRANSFORMER MODEL 
*** MAIN CIPCUIT 

* CURRENT SENSE FOR COLLECTOR CURRENT VCOLLE 
V0C1 1 0 240VOLT 

Q1 9 2 3 QNPN 
VCOLLE 1 9 1UV 
****** AOJUST Rll ******* 

Rll 1 8 3.9000HM 

01 8 4 OMOO 

R1 1 4 10.00HM 

« CHANGE C FOR SN'JEBER ////// 

Cl 4 3 0.33UF 

02 3 1 OMOO 
RSI 3 5 0.0274 

* CHANGE LLS1************************** 

LLS1 5 7 4.00UH IC=198A 
» CURRENT SENSE FOR BIFILAR MINDING VBIFI 
VBIFI 7 6 1UV 

* 

VOC2 0 11 240V 
Q2 13 12 11 QNPN 
****** AOJUST R22 ****** 

R22 13 16 3.9000HM 

03 16 14 OMOO 

R2 13 14 10.00HM 

* CHANGE C FOR SN'JEBER /////// 

C2 14 11 0.33UF IC=450V 

04 11 13 OMOO 
RS2 13 15 0.0274 

* CHANGE LLS2*************************** 

LLS2 15 6 4.00UF 

* CHANGE LM3*********LMB+LLSl(LLS2)=225.6UH AT 20 HZ, 294UH AT 160HZ 

* *******************t.M=6.51MH, LLR=241UH AT 20 HZ 

» *******************lm=26.8MH, LLR=323UH AT 160 HZ 

LMB 6 21 221. 6UH IC=198A 
LM 21 0 6.51MH IC=85A 
LLR 21 22 241UH IC=113A 
RR2 22 0 0.251 
*** GATE DRIVE 

RG1 2 7 10HM 
RG2 12 17 10HM 

VG1 7 3 PML (0 10V 100US 10V I05US -10 420US -10 425US 10V) 

VG2 17 11 PHLtO -10 120US -10 125US 10V 400US 10V 405US -10) 

*«* DEVISE MODEL 

•MODEL QNPN NPN( VA=2000V,VB=2000V) 

.MODEL OMOO 0(BV=2000V) 

»** OUTPUT 

.PRINT TRAN V(9,3),I(VC0LLE),V(4,3),I(VBIFI),V(6) 

.PLOT TRAN V( 9,3 ) ,I( VBIFI ) 

•TRAN l.OUS 130US UIC 
1 .OPTIONS LIMPTS=201 LIST NOOE 
• END 


Fig. F-8 Source listing of SPICE simulation program 


SPI00060 


SPI00070 

SPI00080 


SPI00170 

SPI00190 


SPI00280 


SPI00330 

SPI00380 

SPI00410 
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F.4 LOSS CALCULATION* 

Switching Loss 

There are four possible switching operations for which different loss 
calculations are required. Fig. F-10 shows the switching devices and the 
snubber circuit components only for the A phase. I A1 and I A2 stand for the 
stator currents for bi filar windings A1 and A2 respectively, and V C1 and 
V C 2 are the voltages across capacitors C 1 and C 2 respectively. The follow- 
ing four cases are considered to calculate the switching loss only for 
swi tch SI . 

(1 ) Turn-On of SI When > 0: Initially, while switching device SI is 

open, diode 04 conducts and capacitor voltage is charged up to 2 Vq£. 

The stored energy in this capacitor is 


"c ■ 1 c < 2 V oc > 2 = 2 CV DC 2 < 17 > 

Some of this energy will be dissipated in SI during its turn-on and some 
will be dissipated in the snubber R1 . Assuming an exponential drop of 
collector-emitter voltage v $ ^ as SI turns on. 


v 


SI 



* "Analysis and Design of Optimized Snubber Circuits for dv/dt Protection 
in Power Thyristor Applications," S. J. Wu, GE 660.24, 2/71. 

"Optimum Snubbers for Power Semiconductors," William McMurray, 1971 IAS 
Conference Record, pp. 885-893. 

"Snubbers for Pulse-Width Modulated Bridge Converters with Power Transis 
tors or GTOs," Tore M. Undeland, IPEC Tokyo '83 Conference Record, pp. 
313-323. 
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Fin. F- 10 Switching devices and snubber circuits onlv for "• nhase 



1 


C J 


where t is the exponential time constant. Then the aischarging current i^ is 


- V C " V S1 _ f^DC / -t/x s 
S' R R ' E 



(19) 


where t<- = RC. 


The energy dissipated in SI is 



00 


■ 

o 


V S1 ’si dt 



T T, 


T S + T 


( 20 ) 


Since all of the initial stored energy in the snubber capacitor must be 
dissipated, the total turn-on switching energy loss is given by 
(2) Turn-Off of SI when I A1 > 0: While switch SI is on, I A1 flows through SI. 

The energy initially stored in the equivalent leakage L is 

W L =i-Li L 2 (0) (21) 

where i L (0 ) is the current I A1 at the instant SI turns off. 

The final energy stored in C after the switch is open is 


“c ■ 7 c < 2 v oc > 2 ■ 2Cv dc 2 < 22 > 

An amount of energy equal to that in (22) is dissipated in the equivalent charg- 
ing resistance for capacitor Cl. Some of the initial stored energy in L is 
returned to the dc source, but the worst case turn-off switching energy loss is 

“l - “c = 7 L 'l<°> * 2cv dc 2 


(23) 
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(3) Turh-On of SI When I A2 > O: While switch S2 is on, I A2 flows through S2 
and capacitor Cl is charged up to 21/^. Even though S2 is turned off and 
the base current is supplied to switch SI, SI is still open. If \ l^\ is less 
than 48A, first capacitor Cl discharges its stored eneray, and after the 
capacitor voltage V C1 minus the voltage drop across R1 goes to zero, diode D2 
conducts. If |I A1 | Is larger than 48A, then diode D2 conducts very quickly 
after switch S2 is turned off. 

The worst case energy loss is estimated to be equal to the stored enemy 
in the capacitor. Thus, the loss energy is 

Vi 1 < 2 v dc> 2 ■ 2cv dc 2 < 24 > 

(4) Turn-Off of SI when I^ < 0: While S2 is open, diode D2 is on, and capaci- 
tor voltage is zero. When S2 turns on, 1^ extinguishes and is charged 

up to 2V 0C . 

The loss calculation is very similar to case (2), and the total loss 
energy is 

W L + W c = y L i L 2 (0) + 2CV qc 2 (25) 


On-State Loss 

The total on-state loss is the integral of the instantaneous power 
dissipated in the four switch branches, which is 


P 0N = T 


U 0N,1 


{ V S1 hi dt + 


t 0N,l +t 0FF,l 


V D1 i Dl dt 


t 0N,l 


r^ON, l +t 0N,2 +t 0FF,l 


V S1 n Sl dt + 


1 '0N,l +t 0FF,l 


X 4 


= W(V S1 I S1 ) + (1-d) V M I 01 J X 4 


(26) 
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where d A duty cycle of switching 

V_. A on-state voltage drop across SI (assumed constant) 

0 1 

I S1 A average current through SI 

V Qi A diode drop (assumed constant) 

Ipi A average current through diode 

and (26) assumes that the average currents over the cycle through the 
transistor and the diode are equal. 

Summary 

The total switching energy loss w,^ for one switching is given by addina 
the turn-on energy loss, (17), and the turn-off energy loss, (23) or (25). 


w -II 
w SW 2 L 


2 

( 0 ) 


+ 4CV 


DC 


and the total switching loss over a period j is 



(27) 


1 I t ^ f + ary ^ ~ 
1 L a, RMS S HLV DC T S 


X 4 


(28) 


The total on-state loss is given by (26). Thus, the total inverter loss 
is the sum of (26) and (28) 


P T = 


7 L *L?RHS f S + 4CV 0C 2 f S + dV Sl 'SI ' °- d) V Dl ] Dl] X 4 
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Calculated Losses 

13.1 kW output a t S0% speed 

I n j 1 $ = 53. 2A; I p = 103A; f $ = 1.6 kHz 
d = 0.72/2 = 0.36; V $] = 1.8V; V m =0.8V 

! ave = 48 - 8/2 = 24.4A/device 
On-State Loss from (26) 

P Qn = [0.36 x 1.2 x 24.4 + 0.64 x 0.8 x 24.4] x 4 = 113.2 
Switching loss from (28) 

P et( = [0.5 x 4 x 10‘ 6 x 53. 2 2 + 4 x 0.33 x 10" 6 x 240 2 ] x 1 .6 x 10 3 
sw 

x 4 = 522.8 
Total: 636.0 W 
Efficiency: 95.4% 

16.4 kW output at 10°/ speed 

1 rms = 200A; I p = 245A; f $ = 3.4kHz 

d = 0.6/2 = 0.3; V $1 = 1.8V; V Q1 = 0.8V 
= 197.4/2 = 98. 7A/ device 

aVc 

On-State Loss from (26) 

P on = l0 ' 3 x ] - 8 x 98 * 7 + °- 7 x °- 8 x 98 - 7 l x 4 = 434 - 3 
Switching Loss from (28) 

P sw = [0.5 x 4 + 10‘ 6 x 200 2 + 4 x 0.33 + 10' 6 x 240 2 ] x 3.4 x 10 3 
x 4 = 2,122.0 
Total: 2,556.3 W 
Efficiency: 86.5% 
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60.5 kW Output at 30% Speed 

! rms = 175A 5 ! p = 245A ; f $ = 160Hz 
d = 0.83/2 = 0.42; V $1 = 18. V; V Q1 = 0.8V 
T ave = 173/2 = S6.5A/device 
On-State Loss from (26) 

p on = 10,42 x 1,8 x 86,5 + 0,58 x 0,8 x 86>5] x 4 = 422J 
Switching Loss from (28) 

P sw = [0.5 x 4 x 10" 6 x 175 2 + 4 x 0.33 x 10' 6 x 240 2 ] x 160 x 4 = 
87.9 

Total : 510.0 W 

Efficiency: 99.2% 


F.5 CONTROLLER 

The microprocessor based controller system must generate the required 
driving signals for the power transistors to implement current controlled PWM. 
The hardware necessary for this controller is shown in Fig. f- 11 which consists 
of the main CPU, ROM, RAM, I/O ports and other electronic components. The 
controller sets the output frequency of the inverter drive, in reference to 
the driving command, from 0-200Hz. In addition, the controller sets the ampli- 
tude and the envelope of the controlled current at any given load and speed to 
provide maximum efficiency with acceptable torque pulsations. The closed loop 
control diagram is shown in Fig. F-12. 



Ad re 



Fig. F-ll System hardware 
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G.l DERIVATION OF EQUATIONS FOR TWO-PHASE MACHINE 


The circuit diagram of the two-phase inverter/bifilar induction motor 
is shown in Fig. 2, section 3.0 of this report. Fig. 6-1 shows the expected 
phase voltage and current waveforms without PWM current control. The phase 
relationships between the stator and rotor windings are shown in Fig. G-2. 

The derivation of the simulation model for the two-phase induction motor 
proceeds as follows: 

Let L 1SA = L 1SB = L 1S’ L lrA = L lrB = L lr 


L r “ hr 


L m is the mutual inductance of A-phase or B-phase 
i^S» igS Have the frequency of supply (w) 
i Ar » ^Br have the fre( l uenc y the slip (w-^ r ) 


V AS = VaS + L 1S TF + it X AS 

V BS = ^BS + L 1S “dF + it X BS 

V Ar = R r i Ar + L lr IF + it X Ar 
V Br = R r 1 Br + L lr “iF + it X Br 


X AS " VaS + VAr cos 9 " L m 1 Br sin 0 
X BS = VbS + Vat 5in 6 + Vfir cos 9 
X Ar “ ^Ar + Vas cos 9 + VbS sin 0 
X Br = VBr " VaS sin 9 + VbS cos 9 


( 2 ) 







f 
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Equations (1) and (2) can be expressed in matrix form as follows: 




L m 

0 

L 2 = 

L COS0 
m 

-L sin0 
m 

0 

L m 


L sin© 

L cos9 


m 


m 

m 



L cos 0 
m 

L sin0 
m 


L m 

0 

4 4 

L 4 4 



-L,sin9 

Lcos9 


0 

L 


m 

m 



m 
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If i ^ and i gr are transformed to d-q coordinates which are fixed on the stator, 
that is, d-q axes rotating in the reverse direction relative to the rotor, then 

I I 

i Ar and i gr with angular velocity (co-co r ) become i Ar and ig r with angular velocity 

(«). 

( 5 ) 



r at 1 ! 

• 

= Z t ' + Z T„ 

r = 

r 

r 

r r 


cos(-0) 

-sin(-e) 

cose sine 

T - 


= 



sin(-e) 

cos(-9) 

-sine cose 


cose 

-sine 

• 

- cosine 

-w r cose 

sine 

cose 

Z = 

o r cos6 

-ursine 
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4s + 4 ' L 2 


i s 

T'kj c - 1 (T lr+ L 4 ) 


• • 

_ *>» _ 

C i + C i 
r • r 


.© 


0 + I 

■s 4 


C -1 L 


© 


4 C + L 2 C 


© 


© 

c _ 1 (R r + 4 )c 

+ c _1 (r lr +L 4 )c 

© 


4 s + L 1 L 2 C I 

'V ®t ".r~ "® H 

c 4 C < L lr +L 4> C 


■f 

lr 


1 


-t 

■\r 


( 6 ) 


O 


• 

R r 

0 

• 

0 

0 

+ 4 = 

b 

0 

R s 

Therefore, = 

0 

0 


(7) l 2 xt + L 2 r = 


-cj 1 sine 
r m 


oj l cose 
r m 


-co L cose 
r m 


-co L s i n6 
r m 


cose 

sin9 

-sine 

cose 


Lcos6 

m 

-L sin6 
m 


-co r sine 

-co cose 
r 

Lysine 

lxos6 


co cose 

-corsine 

m 

m 


r 

r 


* |0| 


© t’v 


cose 

sine 


■sine 

cose 


'“r L m sine “r L m c ° s9 
-“r L n, cos9 -“r L m s1ne 


0 co L 
r m 


Vm 0 
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• • 

© c -1 (R r + l 4 ) c + c “ 1 (E |r + L 4 ) c 



R r 

0 


COS0 

-sine 


L r 

0 


-a) sin6 a) cose 
r r 


0 


+ 

sine 

COS0 


0 

L 


-co cose -a) sine 



r 




r 


r r 



R 

0 


0 CO 1 


R„ 

co 


r 



r r 


r 

r r 

= 



+ 


s 




0 

R 


-co L 0 


-co L 

R„ 



r 


r r 


r r 

r 


© 4s + 4 ■ 

4s 0 

+ 

r~ 

3 

o 

i 

L s 0 


0 4s 


0 L 

| m 


0 L S 


© L,C = 

L m cos9 "Sn 5 ^ 6 


cose sine 

= 

L 0 

m 


Lsine L cose 


-sine cos6 


0 L 


m m 




m 



COS0 

sin9 


-sin0 

cose 


L„cos8 

m 

-L m sin8 

m 


L_sin8 

m 

L cos0 
m 




c _ 1 (r lr + l 4 )c = 


o 


0 


L 


r 


So equation (?) becomes 


Z _1 1 


0 

Re 


Cl) L 

r m 


-coJ-_ 0 
r m 


(O..L,.. 
r r 


-co L R 
r r r 
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L S 0 I L m 0 





• 


• i 

‘AS 

;• 

i = 

V 

• 

r 


‘bs 


• i 

i Br 


(7) 


di 


*AS £ ~dT = pi AS 




If the machine has a squirrel -cage rotor, then 



Therefore 


V AS 


R s +L s p 

0 

L p 
m r 

0 


‘as 

V BS 


0 

R S +L S P 

0 

L m p 


‘bs 

0 


L m p 

“r L n, 

R +L p 
r r r 

r r 


‘Ar 

0 


'Vm 


-w_L 
r r 

R +L p 
r r 

■ 

‘Br 


( 8 ) 


The commonly used torque equation is 

T = pL m ^BS *Ar ’ ‘AS ‘b^ ^ 

Equations (8) and (9) are the mathematical model for simulation of the two-phase 
machine. Equation (8) can be changed to the form 


V AS 


R S 

0 

0 

0 


j AS 

V BS 


0 

R s 

0 

0 


^S 

0 


0 

0 

R r 

0 


‘Ar 

0 


0 

0 

0 

R r 


‘Br 


L K + L m 

1 5 m 


m 

0 


m 


L, +L 0 
1 s m 


L lr +L m 


0 


‘as 

L m 


‘bs 

0 


‘Ar 

lr +L m 


S_ 
- CO 
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0 0 0 0 


Us 

0 0 0 0 


Us 

0 “rU 0 VhrV 


Ur 

"“rU 0 -VVLJ 0 


Ur 


( 10 ) 


where 


L s 

0 

L 

m 

0 

-1 

L r 

0 -L m 0 
m 

0 

L S 

0 

Sn 


0 

0 -L 
r m 

Sn 

0 

L r 

0 

A 1 

LrL -L 2 

~ L m 

0 L s O 

0 

L m 

0 

L r 

S r m 

0 

-U° L S 


L S - L 1S + L m 


L r ” Hr + L m 


Then, (10) may be Written 

' Us 


BS 


Ur 

Ur 


L S L r - hn 


L S L r " L m 


L r 

0 

- L m 

0 


V AS 

0 

L r 

0 

" L m 


V BS 

-u 

0 

L S 

0 


V Ar<°> 

0 

- L m 

0 

L S 


V Br<°> 

| L r R s 


, 2 

" L m “r 

~L R r 
m r 

_L mh 


e 

T 

3 ~ 

1 

L r R s 

L L “ 
m r r 

-L_R„ 
m r 

"Vs 

L sU“r 

*-S^r 

L S L r“ 

- l sW 

"UU 

” L S L r“r 

4 R r 


'AS 


BS j 
I i 1 

hr I 


i i 

Url 


(11) 
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L C L - 

L 2 

Sr 

m 

1 


L C L - 

L 2 

S r 

m 

1 


L C L - 

1 2 

S r 

m 

1 



? (L J/ - L-RcI’ac + L m 6J r 1 R«; + + Ry) 


" ^Rclnc " L L tl) 1 . + L R i D „ ) 


L S L r - L m 


A 22 = 

L S /A 50 

A 41 ^ A 31 ’bS + A 26 x “r 

A 26 " 

2rrL r 

A 42 = ^ A 31 ’AS + A 26 ’a^ x “r 

A 31 = 

27rL m 

1 

A 47 = A 42 ‘ R rV 

A 36 = 

L r /A 50 

A, n = 1<-L - L 2 
bO 5 r m 

1! 

00 

co 

<c 

L m /A 50 

A 53 = V BS ‘ R s’bS 

A 39 ~ 

P x A 31 / 2 tt 

A 61 = ' A 41 ' R r’ Ar 


A 40 “ V AS ‘ VAS 


^AS = A 40 X A 36 " A 61 x A 38 = x ^ 

’BS = A 53 x A 36 " A 47 x A 38 = x 

’Ar = A 61 x A 22 ‘ A 40 x A 38 = x ^ 

’Br = A 47 x A 22 ' A 53 x A 38 = x ^ 


T = PL m ^BS ’ 


BS ’Ar _ ’AS W = A 39 ^BS 1 Ar “ ’AS ’Br^ = 
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The simulation diagram for (12) is shown in Fig. 6-3. 

The following sections of this appendix include simulation programs and 
results for different operating points, using the mathematical model of (12). 



Fig. G-3 Simulation diagram of two-phase 


induction motor 



















G. 2 TWO-PHASE SIMULATION RESULTS FOR 10% SPEED, 16 KW 


Source Program 
Plots 

Stator Voltage 
Stator Current 
Output Power 
Current Spectrum 
Power Spectrum 


Partial Listing of Plot Data 
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FILE: LS MATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK PA6E 001 


//IS JOB ( XXXXLH ) • 'GOOOARZI * , 

// PASSWORD =XXXX 
/* JOBPARM Ts5,L=6 

//"ROUTE PRINT UMVMA.R12 
// EXEC WATFIV 
//GO.SYSIN DO * 

4 JOB , NOLIST, T=C5,0J 

C***»*»*»* **»»*# *»*<•*««««« «<t««*K*«*«xx»*««*«*»**«**««*inf **•*»«»»«»*««» 
C « 

C 10X SPEED, 16KW JPL PROJECT * 

C * 

C*******************************"***************"*****"**"************ 
DIMENSION ARK 800 ) ,AR2( 800 ) ,AR3( 800 ) , Y( 20 ) ,X( 20 ) ,GA( 20 ) ,PE( 20 ) , 

1 XXI 4, 20 ),H(9)>ZZ1(20 ) >ZZ2( 20). ZZ 3(20), IF 1(7 ) , AR4< 800 ), 

2 AR5t 800 ), AR6( 800 ),AR7( 800), AR81 800) 

LOGICAL* 1 LR1 ( 800 ) , LR2( 800 ) , LR3( 800 ) , LR4( 600 ) , LR5( 800 ) , LR6( 800 ) , 

1 LR7( 800 ) > LR8( 800 ) 

REAL LOS , LOR . LM , IAS . IBS. ICS . IAR , IBR . ICR , IDS , IQS . IDR , IQR , LMDQ , INS , 

1 POUT .IAREF . IBREF , LI . L2 .M.IMAX. IMAGD , IMAGQ . IBU , 

2 MAGI ( 512 ) . MAG 2 ( 512 ) 

COMPLEX Fl( 512 ) »F2< 512 ) 

C WRITEI 6 .50 ) 

C 50 FORMAT!/, 3X, 'TIME' .8X, 'VDS' .8X, '105' ,8X, 'VQS' ,8X, 'IQS' ,8X, 

C 4 'IDR' ,9X, 'POUT' ) 

IMAX=205. 

SLIP-. 085 
E=240. 

DI=40 

NR UN- 1 

ISTOP=0 

KEEP=1 

IFL(1)=-1 

NINT=7 

IFNT1=0 

IPNT2=0 

IFNT3=0 

IPNT4=0 

IPNT5=0 

IPNT6=0 

IPNT7=0 

C ***************** O************************************************* 
C PARAMETERS FOR TWOPHASE INDUCTION MOTOR 
C ***»**»*»**»» 

L0S=.2256E-03 

LOR=.2410E-03 

LM-6.5100E-03 

RS=.0274 

RRs.0213 

PAIs3. 141593 

P=l. 

FREQ=20. 

OMEGA=2»PAI*FREQ 
OMEGR=OMEGA*( 1 . -SLIP ) 

L1=L0S»LM 

L2=L0R*LH 
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FILE: L$ NATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK PAGE 002 


M=LM 

FACT=1./(L1«L2-M*M) 

FACT1=L2«FACT 

FACT2=L2*RS»FACT 

FACT22=L1*RR»FACT 

FACT3=M*M*OMEGR*FACT 

FACT4=M»L2»OMEGR«FACT 

F ACT4<*=N«L 1 *OMEGR«FACT 

FACTS-M^RRoFACT 

FACT55=M«RS*FACT 

FACT6=M»FACT 

FACT7=L1»L2«0MEGR«FACT 

1 IFLC2)=-1 

2 CALL RKF0UR(NINT,Y,X,GA,PE,XX,H,ZZ1,ZZ2,ZZ3,TIME,PNTS,KEEP,IFL) 
IF(IFU4> J3.4.5 

3 INIT=1 
TIME1=-1 
NC=1 
NCC=1 
TIHE=0. 

FINTIM=2.0500 

PTMAX=800. 

Y( 1 )=0. 

Yt 2 1=0. 

YC3)=0. 

Y(4)=0. 

Y(5)=0.0 
Y(61=0.0 
Y(7)=0.0 
AREAIN=0.0 
ARAOUT=0.0 
H( I )=0.2000E-04 
C H<7)=-1 
VDS=0.0 
VQS=0.0 
IN=0 
SUM-0. 

sum=o. 

NUMBER=0. 

SH=0.0 

SUITCH=0.0 

PSW=0.0 

TOTSW=0.0 

IAE=-1 

IBE=-1 

DO 33 1=1.512 
Fill 1 = ( O.O.O.O ) 

F2( I )=( 0.0.0.0 1 
33 CONTINUE 
C 

4 IFtTIME-TIMEl 17,8.7 
7 TIMEl=TIME 

C 

OEG=FREQ*TIHE 

IDEG=DEG 
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FILE: LS WATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK PAGE 003 


DEG=(DEG-I0EG)»360. 

C 

IAREF=IMAX*SIN< OMEGA*TIME > 

IBREF*IMAX*SIN( 0HEGA»TIME-PAI/2. ) 

C 

XF( IAREF.GT.O. 0 ) IAREF=IMAX 
IF ( IAREF . LT.0.00 ) IAREF=-IMAX 
IFdBREF. GT. 0.0) IBREF=IMAX 
IFdBREF. LT.0.00) IBREF=-IHAX 
C 

DIA=Y< 1 1-IAREF 
DIB=Y(2)-IBREF 
C 

VDS=E 

VQS=E 

IF( IAREF. LT.O ) VOS=-E 
IFdBREF. LT.O) VQS=-E 
C 

IF(TIME. LT.O. 35000) GOTO 8 
IFIOIA. GT.DI. AND. IAREF. GT. 0.0) GOTO 2 4 
IF( DIA.GT.DI ) IAE=1 
GOTO 27 
24 IAE=1 

IF (TIME . LT.2.0000 ) GOTO 27 
SN=SW*1 

PSM=ABS( ( 240. «Y( 1 )*5.5E-06 )/(6»0. 050000 ) )*PSW 
27 IFIOIA. LT.(-DD) IAE=0 

IF( VOS. GT. 0 . 0 .AMD . IAE . EQ. 1 ) GOTO 18 
IF(VOS.LT.O.O.AND.IAE.EQ.O) GOTO 18 
GOTO 22 
18 VDS=-VDS 

IFITIME.LT. 2. 0000) GOTO 22 
SWITCH=SWITCH+1. 

22 IF(DIB. GT.DI) IBE=1 

IF( DIB. LT. ( -01 ) ) IBE=0 

IF! VQS.GT.O .0. AND.IBE.EQ. 1 ) GOTO 19 
IFIVQS.LT. 0.0. AND. IBE.EQ.O) GOTO 19 
GOTO 23 
19 VQS=-VQS 

23 IFITIME. GT. 1.90000. AND. TIME. LT.l. 9001 ) AREAIN=Y(6) 

IF( TIME .GT . 1 . 90000 . AND . TIME . LT . 1 . 9001 ) AR AOUT=Y( 5 ) 

IFITIME. GT. 1.90000. AND. TIME. LT.l. 9001 ) ARMS=Y(7> 

C 

IFITIME. LT.2.0000) GOTO 8 
TOTSW=TOTSW*l. 

C 

8 CONTINUE 
C 

Xd)=FACTl«VDS-FACT2»Y< 1 UFACT3«Y( 2 )*FACT4*Y( 4 )*FACT5*Y( 3) 

XI 2 )=FACT1»VQS-FACT2*Y< 2 )-FACT3*Y( 1 )-FACT4*Y< 3 )*FACT5«YI 4 ) 
X(3)=-FACT6«V0S*FACT55*Y< 1 )-FACT44«Y( 2 )-FACT7*Y( 4 )-FACT22*YI 3) 
X(4)=-FACT6«VQS*FACT55*Y(2)+FACT44»Y(1 )*FACT7»Y(3)-FACT22»Y<4) 
X(5)=LM*(Y(2)»Y(3)-Y(1)»Y(4))«0MEGR/.15 
XI 6 )=( Y( 1 )»VDS*Y( 2 )«VQS )/ . 15 
XI 7)=Y( 1 )#Y( 1 )/. 15 



ooooo non oonn n n onno 


147 


FILE l LS 


HATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PAGE 004 


C 

GOTO 2 

5 IPR=0 
NCsNC-l 
IF(NC)9,9.10 

9 NC=NCC 
IPR=1 

10 IF ( TIME . GE . FINTIM ) ISTOP=l 
IF( ( IPR+ISTOP ) . EQ. 0 ) GOTO 11 

C 

IF( TIME. LT.2.0000 (GOTO 11 
I0S=Y(1> 

IQS=Yt2) 

I0R=Yt3> 

IQR=Y<4) 

IMAG0=I0S*IDR 
IMAGQ=IQS*IQR 

TQUM»< IQS*IOR-IDS*IQR I 
POUT-TQ«OMEGR 
PIN=IOS*VOS*IQS»VqS 
SUMI=SUM1*PIN 
SUM-SUM*POUT 
NUMBER-NUMBER* 1 . 

IF( TIME. LT.2.00 1 GOTO 38 
IN=IN*1 

Fit IN l^CMPLXt Yt 1 1 >0 .0 I 
F2( IN)=CMPLXt POUT >0.0 ) 

38 CONTINUE 

WRITEt 6,250 )DEG f IDS. VDS> POUT 
IFtTIME.LT. 2. 000351 POUT=0.0 

CALL FPLOTt 800 , IPNT1 , AR1 , LR 1 , ISTOP, 1 , 1 .TIME , IDS 1 
CALL FPLOTt 800 . IPNT2 , AR2 . LR2 . ISTOP .1.1 .TIME , VDS ) 
CALL FPLOTt 800 , IPNT4.AR4 , LR4 . ISTOP , 1 . 2 .TIME .IQS) 
CALL FPLOTt 800 , IPNT4 , AR4 , LR4 , ISTOP, 2 , 2 , TIME , VQS > 
CALL FPLOTteOO,IPNT3,AR3.LR3, ISTOP, 1,2, TIME, IMAGO) 
CALL FPLOTt 800, IPNT3.AR3.LR3, ISTOP, 2, 2, TIME, IDR) 
CALL FPLOTt 800, IPNT3.AR3, LR3, ISTOP, 1,1, TIME, POUT) 

11 INIT=0 
IFf ISTOP >1 , 1 ,6 

6 CONTINUE 

250 FORMATt ' MX.F6. 1 , 3t F9. 1 , 2X ) ) 

AVER =SUM/NUMBER 
AVER1=SUM1 /NUMBER 
EFF=AVER/AVER1 
PIN=Yl 6 J-AREAIN 
P0UT=Yt5)-ARA0UT 
RMS=t Yt7)-ARMS)*«0.5 


NOTE .7 IS THE VOLTAGE DROP ACROSS THE OIOOE 
1.85 IS THE CE VOLTAGE OF THE TRANSISTOR 

SWITCH IS THE • OF INTERVAL OF THE SWITCHING t DIODE CONDUCTION) 
TOTSW IS TOTOAL * OF INTEGRATION INTERVALS 
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FILE: LS HATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK PAGE 005 


C PDIOOE IS THE LOSS OF OIOOE DURING CONDUCTION 

C PTPANS IS THE ON TIME LOSS OF THE TRANSISTOR 

C PSH IS THE SWITCHING LOSS OF EACH TRANSISTOR (E*I*TSH/T) 

C TOTSL IS THE TOTAL LOSS OF THE SWITCH 

C ****************»<.(«»**»HHHt*.HHHH».HH»**l*«****.l******* 

POIOOE=! SWITCH/TOTSW 1*0 . 7*IMAX/PAI 
PTRANS=( I -SWITCH/TOTSW)*! .85*IMAX/PAI 
TOTSL=PDIODE*PTRANS*PSW 
T0T4SL=4»T0TSL 
DTRIO=SWITCH/TOTSW 
EFFIO=POUT/I T0T4SL+PIN I 
C 

WRITE! 6. 702 JPOUT.EFFIO, SLIP. IMAX.DTRIO, RMS 
WRITE! 6 .703 )SW. PSW , PTR ANS > PO 100 E » TOT4SL > IN 

702 FORHAT!///,' POWER OUT= * >F8. 1 >2X. ' EFFIO=‘ »F5.3.2X,'SLIP =, » 
1F7.5.2X. ' IMAX=' >F6 . 1 . 2X. *DTRIO=' . F7. 1 . 2X. ' I! RMS )=' »F6.2»///) 

703 FORMAT! /, • #SW/CYC= ’ . F5 . 1 , 8X , • PSW= ' , F6 . 1 , 3X . * PTRANS= * ,F5 . 1 , 

1 2X , • PDIODE = • , F A . 1 , 2X . ' T0T4SL= * ,F6 . 1 , 5X , * IN= *,I5 , /////// ) 

WRITE(6»705) 

705 FORMAT! 1 1 ' »3X, *Z* »9X» ' I HARMONIC * ,9X, *P HARMONIC* ,//) 

C706 FORMAT! 2X,F5.1,5XtF12.8,5X>F15.8) 

ISTOP=0.0 

CC CALL FFT! FI .9) 

CC CALL FFT! F2.9) 

CC DO 401 1=1.150 

CC MAGI! I )=CABSI Fl( I ) )/512. 

CC MAG2! I )=CABSI F2! I ) 1/512. 

CC Z=FLOAT! 1-1 ) 

CC IF 1I.EQ.1) MAGI! 2 )=270. 

CC AMP1=MAG1! I l/MAGl! 2 ) 

CC AMP2=MAG2!I)/MAG2! 1 ) 

C WRITE16.706) Z.AMP1.AMP2 

CC IF! I. EG. 150 ) ISTOP=l 

C 

C CALL FPLOTt 800.IPNT5.AR5.LR5.ISTOP. 1 . 1 .Z.AMP1 ) 

C CALL FPLOTISOO.IPNT6.AR6.LR6.ISTOP.1 . 1 .Z.AMP2 ) 

CC401 CONTINUE 
STOP 
END 

SUBROUTINE FFTIX.M) 

COMPLEX X(512 ) .U.W.T 

N=2**M 

NV2=N/2 

NM1=N-1 

J=1 

DO 7 1=1, NM1 
IFII.GT.J) GOTO 5 
T=XI J I 
X! J)=X(I) 

X!I)=T 

5 K=NV2 

6 IFIK.GE. J )GO TO 7 
J=J-K 

K=K/2 
GO TO 6 
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7 J=J*K 

PI=3. 16159265358979 

00 20 1=1. M 

LE=2**L 

LEl=LE/2 

U=( 1.0, 0.0) 

2=C0S( PI/FLOAT( LEI ) ) 
Y=-SIN(PI/FL0ATCLE1 »> 
M=CMPLX(Z,Y) 

00 20 J=1.LE1 
00 10 I=J,N,LE 
IP=I+LEl 
T=X(IP)*U 
X( IP )=X( I )-T 

io xm=xm+T 

20 U=U*W 

RETURN 
END 

* ENTRY 

1 I CURRENT OF PHASE A STATOR 

1 V VOLTAGE OF PHASE A 

1 P POWER VS TIME 

1 I CURRENT FREQUENCY SPECTRUM 

1 P POWER FREQUENCY SPECTRUM 

tSTOP 
/* 

// 











STATOR CURRENT (10% SPEED, 16 KW) 








(LDIUJC 
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OUTPUT POWER (10% SPEED, 16KW) 




moc-t-zciDi r>»- ;*ci -£ JJ « J. 
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CURRENT SPECTRUM (10~ SPEED. IS KW) 

1 T f 
3 I 

3 I 

3 ! 

H I 





niOCH-ZHD^ O— Z 3 1 j) 1) 
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POWER SPECTRUM (IOC 5PElD,16KW) 

1 C-‘. 

2 1 




o 



0 2 4 6 8 !C 12 14 15 18 
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//PLOTLS JOS (XXXXLH). 'GOOOARZI' > 

// PASSWOPD=XXXX 
/•JOBPARM R=900, B=99000 
// EXEC PLOTSAS 
//SYSIN CD * 

««««««*«***«»«««*«**«»«*«***««««***«*»**««*<HMH*«*«IH>««**««««*»««*»ft* 
* « 
» 10 X SPEED 16KM JPL PROJECT * 

« PLOT DATA * 

• DEGREE CURRENT VOLTAGE POWER * 


GOPTIONS DEVICE=VARIAN; 
DATA ONE; 

INPUT X T1-Y3I 

CARDS; 


0.1 

-184.3 

240.0 

9429.2 

0.4 

-164.8 

240.0 

10725.8 

0.6 

-145.3 

240.0 

12025.3 

0.9 

-126.0 

240.0 

12876.4 

1.2 

-106.7 

240.0 

11486.5 

1.5 

-87.6 

240.0 

10119.5 

1.8 

-68.5 

240.0 

9224.4 

2.0 

-49.6 

240.0 

10584.4 

2.3 

-30.8 

240.0 

11946.4 

2.6 

-12.1 

240.0 

13310.3 

3.0 

14.2 

240.0 

15241.9 

3.3 

32.7 

240.0 

14390.6 

3.6 

51.0 

240.0 

13116.9 

3.9 

69.3 

240.0 

12308.7 

4.1 

87.5 

240.0 

13723.0 

4.4 

105.5 

240.0 

15138.2 

4.7 

123.5 

240.0 

16553.9 

5.1 

148.7 

240.0 

18556.9 

5.4 

166.4 

240.0 

17790.5 

5.7 

184.0 

240.0 

16608.7 

5.9 

201.6 

240.0 

15884.9 

6.2 

219.0 

240.0 

17345.3 

6.5 

236.3 

240.0 

18805.1 

6.8 

246.6 

-240.0 

20059.5 

7.1 

221.9 

-240.0 

20278.0 

7.4 

197.4 

-240.0 

20476.7 

7.6 

172.9 

-240.0 

20227.3 

7.9 

183.4 

240.0 

18895.5 

8.2 

200.7 

240.0 

17797.6 

8.5 

217.9 

240.0 

17147.4 

8.8 

235.0 

240.0 

16629.2 

9.0 

245.1 

-240.0 

19888.2 

9.3 

220.3 

-240.0 

20026.4 

9.6 

195.5 

-240.0 

20148.9 

9.9 

170.9 

-240.0 

20249.9 

10.2 

181.1 

240.0 

19392.7 

10.5 

198.2 

240.0 

18367.0 

10.7 

215.3 

240.0 

17363.8 

11.0 

232.2 

240.0 

18449.1 

11.3 

242.1 

-240.0 

19709.7 
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11.6 

217.1 

-290.0 

19770.9 

11.9 

192.1 

-290.0 

19811.7 

12.1 

167.3 

-290.0 

19833.7 

12.9 

177.9 

290.0 

19028.3 

12.7 

199.3 

290.0 

18080.2 

13.0 

211.2 

290.0 

17159.6 

13.3 

227.9 

290.0 

18265.2 

13.6 

299.6 

290.0 

19777.7 

13. a 

226.3 

-290.0 

20012 9 

19.1 

201.2 

-290.0 

19975.6 

19.9 

176.1 

-290.0 

19919 0 

19.7 

185.9 

290.0 

20750.0 

15.0 

202.7 

290.0 

19869.6 

15.2 

219.3 

290.0 

19011.5 

15 5 

235.8 

290.0 

18567.9 

15. a 

295.3 

-290.0 

19823.9 

16.1 

219.9 

-290 0 

19729.2 

16.9 

199.6 

-290 0 

19605 5 

16.6 

169.9 

-290.0 

19967.9 

16.9 

179.0 

290.0 

18766 9 

17.2 

195.6 

290.0 

17975 6 

17.5 

212.1 

290.0 

17588.6 

17.8 

228 5 

290.0 

19119.3 

16.1 

299.8 

290.0 

20695.3 

18.3 

226.1 

-290.0 

20792.3 

18.6 

200.6 

-290.0 

20536.9 

18.9 

175.2 

-290.0 

20312.8 

19.2 

189.7 

290.0 

19658.1 

19.5 

201.0 

290.0 

18939.9 

19.7 

217.3 

290 0 

18292.9 

20.0 

233.5 

290.0 

19908.9 

20.3 

292.6 

-290.0 

20691.3 

20.6 

216.9 

-290.0 

20371.6 

20.9 

191.2 

-290.0 

20083.9 

21.2 

165.6 

-290.0 

19776.9 

21.9 

175.0 

290.0 

19173.0 

21.7 

191.2 

290.0 

18531.9 

22.0 

207.3 

290.0 

17911.2 

22.3 

223.9 

290.0 

19091.0 

22.6 

239.9 

290.0 

20622.7 

22.8 

220.9 

-290.0 

20582.2 

23.1 

199.6 

-290.0 

20211.5 

23.9 

168.8 

-290.0 

19822.9 

23.7 

178.0 

290.0 

19269 9 

29.0 

199.0 

290.0 

18705.7 

29.3 

210.0 

290.0 

18162.5 

29.5 

225.9 

290.0 

19353.9 

29.8 

291.7 

290.0 

20833.5 

25.1 

222.5 

-290.0 

20772.9 

25.9 

196.5 

-290.0 

20319.0 

25.7 

170.6 

-290.0 

19897.5 

25.9 

179.6 

290.0 

19395.0 

26.2 

195.5 

290.0 

18857.9 

26.5 

211.3 

290.0 

18390.6 

26.8 

227.0 

290.0 

19590.5 
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27.1 

242.6 

240.0 

21114.3 

27.3 

223.3 

-240.0 

20932.6 

27.4 

197.2 

-240.0 

20395.0 

27.9 

171.1 

-240.0 

19517.4 

2S.2 

179.9 

240.0 

18747.8 

28. S 

195.7 

240.0 

18342.1 

28.8 

211.3 

240.0 

18274.2 

29.0 

226.9 

240.0 

19796 9 

29.3 

242.3 

240.0 

21311.4 

29.4 

222.9 

-240.0 

21057.7 

29.9 

196.6 

-240.0 

20436.0 

30.2 

170.4 

-240.0 

19488.9 

30.4 

179.0 

240.0 

18782.5 

30.7 

194.6 

240.0 

18451.9 

31.0 

210.1 

240.0 

18141.3 

31.3 

225.6 

240.0 

19354.8 

31.4 

240.9 

240.0 

20862.5 

31.8 

221.3 

-240.0 

20541.9 

32.1 

194.9 

-240.0 

19841.8 

32.4 

168.5 

-240.0 

19125.3 

32.7 

177.1 

240.0 

18775.6 

33.0 

192.5 

240.0 

18519.2 

33.3 

207.9 

240.0 

18282.4 

33.5 

223.2 

240.0 

19495.2 

33.8 

238.4 

240.0 

20986.8 

34.1 

218.7 

-240.0 

20592.9 

34.4 

192.1 

-240.0 

19308.4 

34.7 

165.6 

-240.0 

18729.2 

34.9 

174.1 

240.0 

18154.0 

35.2 

189.4 

240.0 

17976.7 

35.5 

204.6 

240.0 

17818.6 

35.8 

219.8 

240.0 

19033.6 

34.1 

234.9 

240.0 

20511.6 

34.4 

243.0 

-240.0 

21592.9 

34.6 

216.1 

-240.0 

20728.4 

36.9 

189.4 

-240.0 

19585.0 

37.2 

169.8 

240.0 

17517.2 

37.5 

185.0 

240.0 

17418.1 

37.8 

200.1 

240.0 

17595.1 

38.0 

215.2 

240.0 

19061.1 

38.3 

230.2 

240.0 

20516.4 

38.6 

238.2 

-240.0 

21564.1 

38.9 

211.2 

-240.0 

20615.2 

39.2 

184.4 

-240.0 

18415.1 

39.5 

164.7 

240.0 

16380.2 

39.7 

179.8 

240.0 , 

16363.8 

40.0 

194.9 

240.0 

17572.0 

40.3 

209.8 

240.0 

19012.7 

40.6 

224.7 

240.0 

20442.1 

40.9 

239.6 

240.0 

21860.2 

41.1 

219.5 

-240.0 

20066.5 

41.4 

192.5 

-240.0 

17645.8 

41.7 

165.6 

-240.0 

15240.1 

42.0 

173.7 

240.0 

16036.4 

42.3 

188.7 

240.0 

17460.4 
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42.6 

203.6 

240.0 

18872.6 

42.8 

218.4 

240.0 

20273.1 

43.1 

233.1 

240.0 

20567.1 

43.4 

240.8 

- 240.0 

20240.8 

43.7 

213.7 

- 240.0 

18064.5 

44.0 

186.6 

- 240.0 

16961.5 

44.2 

166.6 

240.0 

16264.5 

44.5 

181.5 

240.0 

17652.0 

44 8 

196.3 

240.0 

19027.3 

45.1 

211.1 

240.0 

19368.4 

45.4 

225.8 

240.0 

19519.8 

45.6 

240.4 

240 0 

19688.1 

45.9 

220.1 

- 240.0 

18738.5 

46.2 

192.9 

- 240.0 

17551.2 

46.5 

165 8 

- 240.0 

16350.9 

46 8 

173.6 

240 0 

17277.2 

47.1 

188.4 

240.0 

17670.4 

47.3 

203.1 

240.0 

17888 2 

47.6 

217.7 

240.0 

18122 4 

47.9 

232.3 

240 0 

19289.0 

48.2 

239.8 

- 240.0 

20194.6 

48.5 

212.4 

- 240.0 

18921 6 

48.7 

185.2 

- 240.0 

17636.0 

49.0 

165.0 

240 0 

15900 9 

49.3 

179.7 

240 0 

16184 6 

49 6 

194.3 

240.0 

16484. 1 

49 9 

208.9 

240.0 

17640.4 

50.2 

223.4 

240.0 

18953.6 

50 4 

237 8 

240 0 

20253.1 

50.7 

217.4 

- 240.0 

19337.1 

51.0 

190.0 

-240 0 

17171.5 

51.3 

169.7 

240.0 

15304.2 

51.6 

164.4 

240.0 

15666.2 

51.8 

198.9 

240.0 

16196.5 

52.1 

213.4 

240.0 

17490.5 

52.4 

227.8 

240.0 

18770.3 

52.7 

242.1 

240.0 

20035.9 

53.0 

221.6 

- 240.0 

19054.5 

53.2 

194.1 

- 240.0 

17473.1 

53.5 

166.8 

- 240.0 

15187.1 

53.8 

174.4 

240.0 

15162.9 

54.1 

188.9 

240.0 

15602.2 

54.4 

203.4 

240.0 

16721.9 

54.7 

217.7 

240.0 

17962.5 

54.9 

232.1 

240.0 

19188.5 

55 2 

239.3 

- 240.0 

19946.9 

55.5 

211.8 

- 240.0 

17808.1 

55.8 

184.3 

- 240.0 

15561.1 

56.1 

163.9 

240.0 

13787.5 

56.3 

178.4 

240.0 

14420.9 

56.6 

192.8 

240.0 

15639.9 

56.9 

207.2 

240.0 

16843.9 

57.2 

221.5 

240.0 

18032.8 

57.5 

235.6 

240.0 

18663.3 

57.7 

242.9 

- 240.0 

18738.1 
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58.0 

215.3 

-240.0 

16643.1 

58.3 

187.8 

-240.0 

15071.4 

58.6 

167.3 

240.0 

13950.0 

58.9 

181.7 

240 0 

15114.6 

59.2 

196.1 

240.0 

16263.9 

59.4 

210.4 

240.0 

16933.3 

59.7 

224.6 

240.0 

17520.7 

60.0 

238.8 

240.0 

18121.2 

60.3 

218.1 

-240.0 

16509.8 

60.6 

190.5 

-240.0 

14869.4 

60.8 

170.0 

240.0 

13684 0 

61.1 

184.3 

240.0 

14804.2 

61.4 

198.7 

240.0 

15908.6 

61.7 

212.9 

240.0 

16623.4 

62.0 

227.1 

240.0 

17273.6 

62.3 

241.2 

240.0 

17936.3 

62.5 

220.4 

-240.0 

16351.2 

62.8 

192.8 

-240.0 

14643.3 

63.1 

165.3 

-240.0 

12927.2 

63.4 

172.7 

240.0 

13536.2 

63.7 

187.0 

240.0 

14596.3 

63.9 

201.3 

240.0 

15357.6 

64.2 

215.5 

240.0 

16071.8 

64.5 

229.6 

240.0 

16797.7 

64.8 

243.7 

240.0 

17584.9 

65.1 

222.9 

-240.0 

16272.8 

65.4 

195.2 

-240.0 

14483.9 

65.6 

167.6 

-240.0 

12687.6 

65.9 

175.0 

240.0 

13190 5 

66.2 

189.3 

240.0 

13960.5 

66 5 

203.5 

240.0 

14741.6 

66.8 

217.7 

240.0 

15533 6 

67.0 

231.8 

240.0 

16493.3 

67.3 

238.9 

-240.0 

16998.4 

67.6 

211.1 

-240.0 

15136.6 

67.9 

183.5 

-240.0 

13243.6 

68.2 

162.9 

240.0 

11709.8 

68.5 

177.2 

240.0 

12549.6 

68.7 

191.5 

240.0 

13399.8 

69.0 

205.7 

240.0 

14336.7 

69.3 

219.9 

240.0 

15274.1 

69.6 

233.9 

240.0 

16194.4 

69.9 

241.0 

-240.0 

16616.5 

70.1 

213.2 

-240.0 

14638.9 

70.4 

185.5 

-240.0 

12682.5 

70.7 

165.0 

240.0 

11220.4 

71.0 

179.3 

240.0 

12133.1 

71.3 

193.5 

240.0 

13030.0 

71.5 

207.7 

240.0 

13909.7 

71.8 

221.9 

240.0 

14779.9 

72.1 

235.9 

240.0 

15687.0 

72.4 

243.0 

-240.0 

16130.0 

72.7 

215.2 

-240.0 

14223.4 

73.0 

187.5 

-240.0 

12321.3 

73.2 

166.9 

240.0 

10813.7 
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73.5 

181.2 

240.0 

11659.6 

73. a 

195.5 

240.0 

12488.1 

74 . 1 

209.7 

240.0 

13325.2 

74.4 

223.8 

240.0 

14290.8 

74.6 

237.9 

240.0 

15265.0 

74.9 

244.9 

-240.0 

15774.6 

75.2 

217.2 

-240.0 

13759.3 

75.5 

189.5 

-240.0 

11706 8 

75. a 

168.9 

240.0 

10123.0 

76.1 

183.2 

240.0 

10932.5 

76.3 

197.5 

240.0 

11950.6 

76.6 

211.7 

240.0 

12976.6 

76.9 

225.8 

240.0 

14010.3 

77.2 

205.1 

-240.0 

12438 5 

77.5 

177.5 

-240.0 

10341.8 

77.7 

184.8 

240.0 

10596.5 

78 0 

199.1 

240.0 

11596.1 

78.3 

213.3 

240.0 

12659.1 

78.6 

227.5 

240.0 

13729.3 

78.9 

241.6 

240.0 

14741.3 

79.1 

220.8 

-240.0 

13070 5 

79.4 

193.1 

-240.0 

10926.0 

79.7 

165.6 

-240.0 

8778.8 

80.0 

173.0 

240.0 

9049.0 

80.3 

187.4 

240.0 

10157.1 

60.6 

201.7 

240.0 

11271.8 

80.8 

215.9 

240.0 

12392.8 

81.1 

230.1 

240.0 

13104.1 

81.4 

244.2 

240.0 

13716.0 

81.7 

223.5 

-240.0 

11973.1 

82.0 

195.8 

-240.0 

10210.1 

82.2 

168.3 

-240.0 

8557.5 

82.5 

175.8 

240.0 

9258.8 

82.8 

190.2 

240.0 

10336.2 

83.1 

204.5 

240.0 

10920.4 

83.4 

218.8 

240.0 

11486.3 

83.6 

233.0 

240.0 

12033.9 

83.9 

212.3 

-240.0 

10770.6 

84.2 

184.8 

-240.0 

9168.2 

84.5 

164.4 

240.0 

8052.6 

84.8 

178.8 

240.0 

9157.0 

85.1 

193.2 

240.0 

9689.1 

85.3 

207.6 

240.0 

10202.8 

85.6 

221.9 

240.0 

10698.1 

85.9 

236.2 

240.0 

11782.5 

86.2 

243.4 

-240.0 

12529.6 

86.5 

215.8 

-240.0 

10985.4 

86.7 

188.3 

-240.0 

9335.4 

87.0 

167.9 

240.0 

7512.7 

87.3 

182.4 

240.0 

7977.9 

87.6 

196.9 

240.0 

8424.7 

87.9 

211.3 

240.0 

8990.4 

88.2 

225.7 

240.0 

10237.1 

88.4 

240.0 

240.0 

11487.9 

88.7 

219.4 

-240.0 

10463.3 



G.3 TWO-PHASE SIMULATION RESULTS FOR 80% SPEED, 60 KW 

Source Program 
Plots 

Stator Voltage and Current 
Output Power 
Current Spectrum 
Power Spectrum 


Listing of Plot Data 
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FILE: OPTIM 


HATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PACE 001 


//OPTIM JOB (XXXXLM ) • 'GOOOARZX' » 

// PASSWORDS??*? 

/*JOBPARM T*4,L=6 
//•ROUTE PRINT UMVMA.EEPR1 
// EXEC MATFIV 
//GO.SYSIN 00 * 

♦JOB 

C********************************************************************* 


c # 

C SOX SPEED 60KM JPL PROJECT » 

C « 

C********* *****•*««««**««*««•«««*»««•«««««•*»••«•»««•*••••«••••»«««•«« 
DIMENSION AR1 ( 600 ) > AR21 800 ) ,AR3( 600 ) ,Y( 20 ) ,X( 2Q).6A( 20) >PE( 20 )• 

1 XX( 4,20 ) ,H( 9 ) .ZZH 20 ) >ZZ2( 20 ) ,ZZ3( 20) ,IFL( 7),AR4(6Q0 )» 

2 AR5( 600 ) , AR6( 600 > , AR7( 800 ) , AR6( 600 ) 

LOGICAL* 1 LR1( 600 ) , LR2( 600 ) , LR3( 600 ) ) LR41 600 ) > LR51 600 1 ,LR6< 800 1 . 

1 LR7( 600 ) )LR6( 600 ) 

COMPLEX Fl( 512 ) i F2( 512 ) 

REAL LOS.LOR.LM.IAS.IBS.ICS.IAR.IBR.XCR.IDS.IQS.IDR.XQR.LHDO.INS* 

1 POUT, IAREF.IBREF, LI. L2.M.IMAX, IMAGO, IMAGO, IBU, 

2 MAGI ( 256 ) ,MAG2( 256 ) 

IMAX=245 

01=15.0 
SLIP=0.02 
C WRITE(6»50) 

C 50 FORMAT! /i3X, 'TIME' i6X, 'VOS' ,6X, 'IDS* »8X, 'VOS' .BX.'IQS* ,6X» 

C 4 ' IOR ' , 9X , ' POUT ' ) 


NRUN=1 

IST0P=0 

KEEP=1 

IFL( 1 )=-l 

NINT=7 

IPNT1=0 

IPNT2=0 

IPNT3=0 

IPNT4=0 

IPNT5=0 

IPNT6=0 

IPNT7=0 

£=240. 

C PARAMETERS FOR TWOPHASE INDUCTION MOTOR 

C *****»«tHt**«* tt*********************************************^******* 
L0S=. 2940E-03 
LORS.3230E-03 
LM=26.600E-03 
R3=.0274 
RRs.0213 
PAX«3. 141593 
P=l. 

FREQ=160. . 

OMEGA=2 . 0*PAX*FREQ 
OMEGR=OMEGA*( l.-SLIP) 

L1=L03*LM 

L2=L0R*LM 
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FILE: OPTIM HATFXV A UNIVERSITY OF MISSOURI COMPUTER NETWORK RASE 00Z 


M*LH 

FACT»1./CLI«L2-M*M) 

FACT1*L2«FACT 

FACT2sL2«RS«FACT 

FACT22=L1*RR»FACT 

FACT3=M*M*0NEGR»FACT 

FACT4=M*L2*ONEGR»FACT 

F ACT4 ASM* L 1 *OMEGfi*F ACT 

FACT5sM*RR*FACT 

FACT5S=M*RS*FACT 

FACT6=M*FACT 

FACT7=L1*L2*0MEGR*FACT 

1 IFL(2)=-1 

2 CALL RKF0UR(NINT»Y»X»6A»PE,XX*H»ZZ1 »ZZ2>ZZ3>TINE >PNT3«KEEP»IFL) 
IF(IFL<4>)3,4,5 

3 INIT=I 
TIMEls-1 
NC=1 
NCC=1 
TIME=0. 

FINTIM=0.2500 

PTMAX=800. 

Y( 1 1=0. 

Y( 2 )=0. 

Y(3I=0. 

Y<4)=0. 

Y(51*0.0 
Y<6)=0.0 
C YC7)=0.0 

AREAINsQ.O 
ARAOUT=0.0 
H( I )= .050E-04 
C H(7)=-I 

VDS=0.0 
V0S=0.0 
SUM=0. 

INsO 

SUMIsO. 

NUMBERsO. 

SW*0.0 

DI1 S 0I 

SWITCHsO.O 

PSW=0.0 

TOTSNsO.O 

IAE»-1 

IBEs-I 

DO 61 1*1 ,256 
PKI)s(O.OtO.O) 

F2( I )*( 0.0 >0.0 ) 

61 CONTINUE 
C 

4 IF(TIME-TIME1 >7,8,7 
7 TIME1»TIME 

C 

OEG=FREO*TIME 



164 


FILE: OPTIC! 


MATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PASE 003 


I0EG=DE6 

DES=C0E6-IDE6)«360. 

C 

IAREF=IMAX*SINC0HE6A*TIME ) 

IBREF=IMAX»SINC 0MEGA»TIME-PAI/2 . I 
C 

IF(IAREF.GT.O.O) IAREF=IMAX 
IF(IAREF.LT.O.OO) IAREF=-IMAX 
IF( IBREF.GT.O.O) IBREF=IMAX 
IFCIBREF.LT. 0.00) IBREF=-IMAX 
C 

OIA=YC 1 J-IAREF 
DIB=YC 2 )-IBREF 
C 

VDS=E 

VQS=E 

IFCIAREF.LT. 0) VD5=-E 

IFC IBREF . LT.O ) VQS=-E 
C 

IFCTIME.LT. 0.02000) GOTO B 
IFCTIME.GT. 0.00000) GOTO 23 

ISW=0 / 

IFCOIA.LT. Dll )ISW=I 
IFCOIA.GT.DI ) DI1=DIA 

IFCDIA.6T.OI.AND.IAREF.6T.O.O.AND.ISH.EQ.O) GOTO 24 
IFCOIA.GT.DI) IAE=1 
GOTO 27 
24 IAE-1 

IFC TIME. LT.O. 24375) GOTO 27 
SW=SWM 

PSW=ABSC C 240. «YC 1 )*5.5E-06 )/C6»6.25E-03) )*PSH 
27 IFCOIA.LT. C-OI) ) IAE=0 

IFCVDS.6T.O.O.AND.IAE.EQ.1 ) GOTO 16 
IFCVOS.LT. 0.0. AND. IAE.EQ.0) GOTO 18 
GOTO 22 
18 VD5=-VDS 

IF CTIME.LT. 0.24375) GOTO 22 
SMITCH=SMITCH+1. 

22 IFC0IB.6T.DI ) IBE=1 

IFCDIB.LT. (-DID IBE=0 

IFC VQS.6T.0.0.AND.IBE.EQ. 1 ) GOTO 19 
IFC VQ5.LT. 0.0. AND. IBE.EQ.O ) GOTO 19 
GOTO 23 
19 VGS=-VQS 

23 IFC TIME. GT. 0.23125. AND. TIME. LT.O. 2312612 )AREAIN=Y<6) 

IFCTIME.GT. 0.23125. AND.TIME.LT. 0.231 26 12 )ARAOUT=Y( 5) 

IFC TIME.GT. 0.23125. AND.TIHE.LT. 0.231 2612 )ARMS-Y( 7) 

C 

IFCTIME. LT.O. 24375) GOTO 6 
T0TSM-T0TSW*1 . 

C 

8 CONTINUE 
C 

XC 1 )=FACT1«VDS-FACT2*YC 1 )+FACT3*YC 2 )«FACT4*YC4 )*FACT5«YC 3) 

XC 2 )=FACT1«VQS-FACT2*YC 2 )-FACT3*YC 1 )-FACT4*YC 3 >*FACT5»YC4 ) 
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FILE: OPTIM 


HATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PA6E 004 


XC 3 )=-FACT6*VDS+FACT55«Y ( 1 )-FACT44«Y< 2 )-FACT7»Y( 4 >-FACT22*Y( 3 ) 
X(4 )=-FACT6«VQS*FACT55«Y( 2 )*FACT44*Y< 1 )+FACT7*Y( 3 )-FACT22«Y( 4 ) 
X(5) = LM*(Y(2)«Y(3)-Y<1 )*Y<4) )«OMEGR/.01875 
XC6 )=< Y( 1 )«VDS+Y( 2 )*VQS )/. 01875 
X(7)=Y(1)«Y(1)/. 01875 
C 

GOTO 2 

5 IPR=0 
NC=NC-1 
IF(NC)9,9,10 

9 NC=NCC 
IPR=1 

10 IF(TIME.GE.FINTIM) IST0P=1 
IF( (IPR+ISTOPJ.EQ.O) GOTO 11 

C 

IFfTIME.LT. 0.23125) GOTO 11 
IDS=Y(1) 

I0S=Y<2) 

IOR=Y(3) 

IQR=Y(4) 

C 

IMAGO=IDS*IOR 

IMAGQ=IQS*IQR 

C 

TQ=LM«( IQS*IOR-IDS*IQR ) 

POUT =TQ*OMEGR 
C PIN=IOS*VDS*IOS«VQS 

C SUM1=SUM1+PIN 

C SUM=SUM+POUT 

C NUMBER=NUHBER*1. 

IF ( TIME. LT. 0.24375) GOTO 31 
IN=IN+1 

CC FMIN)=CMPLX<Y<1),0.0) 

CC F2(IN)=CMPLXCPOUT,0.0) 

31 CONTINUE 

WRITE(6,250)0EG, IDS. VOS, POUT 
C IF(TIME.LT. 0.19250) POUT=0.0 

CALL FPLOT(800 , IPNT1 . AR1 , LR1 .ISTOP, 1 . 1 .TIME ,105 ) 

CALL FPLOTC 800 , IPNT2.AR2 , LR2 , ISTOP .1.1 .TIME . VDS ) 

C CALL FPLOT( 800 .IPNT4.AR4.LR4. ISTOP. 1 .2. TIME .IQS) 

C CALL FPLOT( 800 , IPNT4 . AR4 . LR4 . ISTOP .2.2. TIME » VQS ) 

C CALL FPLOT(800,IPNT3,AR3,LR3, ISTOP. 1.2, TIME, IMAGO) 

C CALL FPLOT(800»IPNT3»AR3» LR3.ISTOP.2.2 .TIME >IDR ) 

CALL FPLOT( 800 , IPNT5 .AR5.LR5, ISTOP ,1,1, TIME , POUT ) 

C 

11 INIT=0 

IF( ISTOP) 1,1, 6 

6 CONTINUE 

250 FORMAT! ' * ,F6. 1 .3IF9.2.2X) ) 

C AVER=SUM/NUMBER 

C AVER1=SUM1/NUMBER 

C EFF=AVER/AVER1 

PIN=Y(6)-AREAIN 
P0UT=Y(5)-ARA0UT 
RMS=( Y( 7 ) -ARMS )*«0 . 5 
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FILE: OPTIH WATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK PAGE 00' 


C *»*»*»**»*#»»#***»»»»#***»#»»*»»**»#**#»**«*)*«»#«*#**»*****###**** 

C NOTE .7 IS THE VOLTAGE DROP ACROSS THE OIOOE 

C 1.85 IS THE CE VOLTAGE OF THE TRANSISTOR 

C SWITCH IS THE • OF INTERVAL OF THE SWITCHING ! DIODE CONDUCTION} 

C TOTSW IS TOTOAL * OF INTEGRATION INTERVALS 

C PDIODE IS THE LOSS OF DIOOE DURING CONDUCTION 

C PTPAN3 IS THE ON TIME LOSS OF THE TRANSISTOR 

C PSW IS THE SWITCHING LOSS OF EACH TRANSISTOR !E*I«T5W/T) 

C TOTSL IS THE TOTAL LOSS OF THE SWITCH 

C #***«**««*»»*«*#*«*«««*«*«*«##*«»**»**#**#***#****«# 

POIOOE=! SWITCH/TOTSW )«0 . 7»IMAX/PAI 
PTRANS= ( 1 -SWITCH/TOTSW . 85»IMAX/PAI 
TOTS L= PDIOO E ♦ PTRANS* PSW 
T0T4SL=4*T0TSL 
DTRIO=SWITCH/TOTSW 
C EFFIO=POUT/! T0T4SL+PIN ) 

EFFIO=POUT/PIN 

SW=SW*SW 

C 

WRITE!6,702)POUT,EFFIO,SLIP,IMAX,OTRIO,RMS 
WRITE ( 6 . 703 ) SW , PSW , PTRANS >P0I00E,T0T4SL>IN 

702 FORMAT!///, * POWER OUT=* ,F3. 1 ,2X, 'EFFIO=* .F5.3.2X. 'SLIP=* . 

1F7.5. 2X. ' IMAX= ' ,F6. 1 ,2X> 'DTRIO=* »F7. 1 .2X. ’ I!RM5) = ' »F7. 1 »//l 

703 FORMAT! / , ' #SW/CYC= * , F5 . 1 , 8X , • PSW= 1 , F6 . 1 , 3X , ' PTRANS= • , F5 . 1 , 

1 2X , 1 PDIOOE= 1 , F4. 1 , 2X, ' T0T4SL= ' »F6. 1 >5X. ' IN=* ,15 ./////// ) 

WRITE! 6, 705) 

705 FORMAT! IX, 'Z* . 1IX, *1 HARMONIC ' ,15X,*P HARMONIC*,//) 

C706 FORMAT! 1X,F5. 1 ,5X,F12.8,5X,F15.d) 

IST0P=0.0 

CC CALL FFT! FI ,8) 

CC CALL FFT! F2 ,8 ) 

CC 00 401 1-1,128 

CC MAGI! I )=CABS! Fit I ) )/256. 

CC MAG2II)=CABS!F2!I))/256. 

CC Z=FL0AT1I-1) 

CC IF II.EO.l) MAGI! 2)=275. 

CC AMP1=MAG1( I 1/MAG1I 2 ) 

CC AMP2=MAG2! I )/MAG2! 1 ) 

C WRITE! 6, 706) Z.AMP1 , AMP2 

CC IF(I.EQ.128) ISTOP=l 

C IF! AMP2.GT.300 ) AMP2=0.0 

C 

C CALL FPLOT!800»IPNT4»AR4,LR4,ISTOP> 1 , 1 ,Z, AMP1 ) 

C CALL FPLOT!800,IPNT3,AR3,LR3,ISTOP,1,1,Z,AMP2) 

401 CONTINUE 

STOP 
END 

SUBROUTINE FFTIX.M) 

COMPLEX XI 256 ) ,U,M,T 

N=2**M 

NV2=N/2 

NMl=N-l 

J=1 

DO 7 1=1, NMl 
IFII.GT.J) GOTO 5 
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FILE: OPTIM MATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PAGE 006 


T=X(J) 

X( J)=X(I) 

X(I)=T 

5 K=NV2 

6 IFIK.GE.J )S0 TO 7 
J=J-K 

K=K/2 
GO TO 6 

7 J=J+K 

PI=3. 1A159265358979 

DO 20 L=1,M 

LE=2»»L 

LEl=LE/2 

U=(I. 0,0.0) 

Z-COSt PI/FLOAT( LEI ) I 
Y=-SIN( PI/FLOATI LEI ) ) 

W-CMPLX(ZtY) 

DO 20 J=1 >LE1 
DO 10 I=JiN> LE 
IP=I+LEl 
T=X(IP)*U 
X( IP )=X( I )-T 
10 X(I)=X(I)*T 

20 U=U*W 

RETURN 
END 

GENTRY 

1 A CURRENT OF PHASE A OF STATOR 

1 V VOLTAGE OF PHASE A OF STATOR 

1 P POWER VS TIME 

1 I CURRENT HARMONIC MAGNITUDE VS HARMONIC ORDER 

1 P POWER HARMONIC MAGNITUDE VS HARMONIC ORDER 

♦STOP 
/* 

// 



-h -r rn au 33 c_ r> 
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VOLTAGE & CURRENT (80% S.GOKW.MOTOR OPERATION) 


2S0-J 



3mim 
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OUTPUT POWER (80% S.60KW, MOTOR OPERATION) 



OEGREE 



CURRCXT SPECTRUM (80~ SPEED, 60 K\V) 


'•°1 t 

i I 

3 { 

c- . 9- l 

1 I 

4 I 

\ I 


i ! 


1 



0 1 2 3 H 5 6 7 8 9 10 1 1 12 13 t M 15 


Hfiftf'ONlC ORCEfl 
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POWER SPECTRUM (80 5 : SPEED, 6GKW) 



HARM2N I C 



172 


FILE: MOTOR PLOT 


A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PAGE 001 


//MOTOR P JOB ( XXXXLH ) , ' GOOD ARZI ' » 

// PASSMORO =XXXX 
/•UOBPARM R=900,B=99000 
// EXEC PLOTSAS 
//SYSIH DD » 

GOPTIONS DEVICE=VARIANl 
DATA ONE! 

INPUT X YI-Y3! 

K»*«*««*«K«**ltM««*Wtt»«*«*«**»**M«»**«lt »«»*** 


* 


» 


* 8 OX SPEED. 60KW JPL PROJECT » 
« PLOT DATA * 
*111 * 

* DEGREE I CURRENT I VOLTAGE I POWER » 
###*###*# | ******** | ********* | ***************************************** 


CARDS! 



- 

0 9 

-233.96 

290.00 

73530 25 

0.9 

-227.92 

290.00 

72589.00 

1.5 

-221.93 

290 00 

71657.56 

2.1 

-215.98 

290.00 

70751.13 

2.9 

-207.63 

290.00 

69503.99 

9.1 

-195.96 

290 00 

67807 31 

5.7 

-179.75 

290.00 

65599.56 

8.0 

-157.91 

290.00 

62615.67 

11.2 

-126.97 

290.00 

59009 33 

15.8 

-86.28 

290.00 

59955.57 

22.9 

-33.59 

290.00 

51216.90 

31.6 

31.90 

290.00 

99532.19 

99.6 

109.19 

290.00 

52899.97 

57.6 

155.75 

290.00 

60361.97 

70.6 

188.23 

290.00 

68587.13 

83.7 

209.60 

290.00 

79058.00 

88.3 

207.21 

290.00 

79693.56 

89.9 

207.80 

290.00 

79708.81 

90.3 

207.92 

290.00 

79138.25 

90.7 

208.03 

290.00 

73965.38 

91.1 

208.13 

290.00 

72802.99 

91.5 

208.22 

290.00 

72199.25 

92.1 

208.33 

290.00 

71292.75 

92.9 

208.95 

290.00 

69995.00 

99.1 

208.57 

290.00 

68298.13 

95.7 

208.61 

290.00 

66039.99 

98.0 

208.99 

290.00 

63103.30 

101.2 

207.81 

290.00 

59992.22 

105.8 

206.29 

290.00 

55929.19 

112.9 

203.07 

290.00 

51671.75 

121.6 

197.75 

290.00 

99950.17 

130.8 

192.98 

290.00 

51719.76 

190.0 

190.93 

290.00 

55985.92 

199.2 

191.68 

290.00 

61636.06 

158.9 

198.18 

290.00 

67998.50 

167.6 

211.29 

290.00 

72177.50 

176.8 

231.96 

290.00 

79603.25 

179.1 

238.95 

290.00 

79719.19 

179.9 

290.87 

290.00 

79698.63 
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FILE: MOTOR PLOT A UNIVERSITY OF MISSOURI COMPUTER NETWORK PAGE 002 


ICO. 5 

236.11 

-240.00 

73889.38 

181.1 

230.08 

-240.00 

72941.25 

181.7 

224.10 

-240.00 

72012.88 

182.2 

218.16 

-240.00 

71104.63 

163.1 

209.82 

-240.00 

69854.19 

184. 2 

198.17 

-240.00 

68153.75 

185.8 

181.99 

-240.00 

65884.75 

188.1 

159.68 

-240.00 

62946.13 

191.4 

129.29 

-240.00 

59323.89 

196.0 

88.69 

-240.00 

55243.95 

202.5 

36.05 

-240.00 

51460.83 

211.7 

-28.73 

-240.00 

49697.88 

224.7 

-101.31 

-240.00 

52920.96 

237.8 

-152.69 

-240.00 

60219.04 

250.8 

-185.00 

-240.00 

68248.75 

263.6 

-201.25 

-240.00 

73508.50 

266.4 

-203.82 

-240.00 

74067.63 

269.6 

-204.24 

-240.00 

74071.44 

270.0 

-204.37 

-240.00 

74058.94 

270.4 

-204.49 

-240.00 

73484.50 

270.8 

-204.60 

-240.00 

72807.69 

271.2 

-204.70 

-240.00 

72140.88 

271.6 

-204.79 

-240.00 

71484.19 

272.2 

-204.90 

-240.00 

70572.56 

273.0 

-205.02 

-240.00 

69317.56 

274.1 

-205.13 

-240.00 

67611.13 

275.8 

-205.17 

-240.00 

65333.88 

278.1 

-205.00 

-240.00 

62384.63 

281.3 

-204.37 

-240.00 

58749.62 

285.9 

-202.80 

-240.00 

54655.75 

292.4 

-199.63 

-240.00 

50861.50 

301.6 

-194.33 

-240.00 

49102.36 

310.9 

-189.59 

-240.00 

50845.82 

320.1 

-187.09 

-240.00 

55110.38 

329.3 

-188.39 

-240.00 

60769.29 

338.5 

-194.96 

-240.00 

66604.94 

347.7 

-208.09 

-240.00 

71369.31 

356.9 

-228.90 

-240.00 

73842.38 

359.2 

-235.41 

-240.00 

73966.63 

359.8 

-237.12 

-240.00 

73962.13 


1 

PROC GPLOT i 

PLOT Y1«X Y2«X /OVERLAY I 
SYMBOL1 I-JOIN L=1 ! 

SYMBOLS I=JOIN L=13i 
LABEL X=0EGREE 

Y1=CURRENT VOLTAGE l 

TITLE VOLTAGE A CURRENT (SOX S.60KW, MOTOR OPERATION)! 
PROC GPLOT! 

PLOT Y3*X /VAXIS=0 TO B0000 BY 5000! 

SYMBOL! I=SPLINE! 

LABEL X=DEGREE 

Y3= POWER I 

TITLE OUTPUT POWER 180X S.60KW, MOTOR OPERATION)! 

/* 



G.4 TWO-PHASE SIMULATION RESULTS FOR 80% SPEED, 12 KW 


Source Program 
Plots 

Stator Voltage 
Stator Current 
Output Power 
Current Spectrum 
Power Spectrum 


Partial Listing of Plot Data 
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FILE: OPTIM HATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK PAGE 001 


//OPTIM JOB (XXXXLH ) > ' JPL PROJ', 

// PASSWORD =XXXX 
/*JOBPARM T=2 , L-4 

//*ROUTE PRINT UMVMA.EEPR3 
// EXEC WATFIV 
//GO.SYSIN DD * 

♦ JOB 

DIMENSION ARK 800 ) , AR2( 800 ) ■ AR3( 800 ) , Y( 20 ) ,X( 20 ) >GA( 20 ) ,PE( 20 ) , 

1 XX(4,20),H<9),ZZK20 ) ,ZZ2( 20 ) , ZZ3( 20 ) , IFL( 7 > , AR4C 800 ) , 

2 AR5( 600 ) , AR6( 800 ) ,AR7( 800 ) . AR8( 800 » 

LOGICAL*! LR1 ( 800 ) , LR2( 800 1 , LR3< 800 ) , LR4( 800 t , LR5( 800 ) , LR6( BOO > , 

1 LR7( 600 ) ,LR6( 800 ) 

COMPLEX Fl( 512 ) ,F2(512 I 

REAL LOS , LOP , LM , IAS , IBS , ICS , IAR , IBR , ICR , IDS . IQS , IDR , IQR , LMDQ , INS, 

1 POUT, IAREF.IBREF, LI. L2.M.IMAX, IMAGO, IMAGQ.IBU, 

2 MAGI ( 256 ),MAG2( 256) 

IMAX=65.0 

01=38.0 
SLIP=0. 00725 

C********************************************************************* 

c * 

C 80X SPEED, 12KW JPL PROJECT * 

C * 

C*«*»««**«K**««*«******««*«»*««««***«*«»««*«»K««W**««««»««»««»»«VK**»» 

C WRITE(6,50) 

C 50 F0RMAT(/,3X, ’TIME' ,8X, 'VOS' >8X, 'IDS' ,8X, 'VQS' ,8X, 'IQS' ,8X, 

C 4 'IDR' ,9X, 'POUT' ) 

NRUN=1 

ISTOP=0 

KEEP=1 

IFL(1)=-1 

NINT=7 

IPNT1=0 

IPNT2=0 

IPNT3=0 

IPNT4=0 

IPNT5=0 

IPNT6=0 

IPNT7=0 

E=240. 

C *******#*«****«***«**««#««»**#******#*«» ***,Hm**m*)»*«tHt*,HHHHt****** 

C PARAMETERS FOR TWOPHASE INDUCTION MOTOR 

C tf****************************************************************** 
L0S=.2940E-03 
L0R=.3230E-03 
LM=26.800E-03 
RS=.0274 
RR-.0213 
PAI=3. 141593 
P=l. 

FREQ-160, 

OMEGA=2 . 0*PAI*FREQ 
OMEGR=OMEGA*( 1 . -SLIP ) 

L1=L0S+LM 

L2=L0R*LM 
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1 

2 

3 


C 


c 


61 

C 

4 

7 

C 
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M=LM 

FACT=1./(L1#L2-M*M) 

FACT1=L2»FACT 

FACT2=L2*RS*FACT 

FACT22=L1*RR*FACT 

FACT3-M»M»OMEGR*FACT 

FACT4=M*L2*OMEGR*FACT 

FACT44=M*Ll*0NEGR»FACT 

FACT5=M*RR»FACT 

FACT55=M«RS«FACT 

FACT6=M«FACT 

FACT7=L1*L2*0MEGR*FACT 

IFL(2)=-l 

CALL RKF0UR(NINT,Y,X,GA,PE,XX,H,ZZ1,ZZ2,ZZ3,TIME,PNTS,KEEP,IFL) 

IF(IFU4))3.4,5 

INIT=1 

TIME1=-1 

NC=I 

NCC=1 

TIHE=0. 

FINTIM=0.2500 

PTHAX=800. 

Y(1)=0. 

Y<2)=0. 

Y(3)=0. 

Y(4)=0. 

Y(5)=0.0 

Y( 6 )=0.0 

Y(7)=0.0 

AREAIN=0.0 

ARA0UT=0.0 

H( 1 )=.050E-04 

H(7)=-l 

VOS=0.0 

VQS=0.0 

SUM=0. 

SUMI=0. 

IN-0 

sum=o. 

NUMBERsO. 

SW-0.0 

011=01 

SUITCH=0.0 

PSW=0.0 

TOTSW=0.0 

IAE=-1 

IBE=-1 

DO 61 1=1,256 
F1(I)=(0. 0,0.0) 

F2(I)=(0. 0,0.0) 

CONTINUE 

IFCTIME-TIME1 )7,fl,7 
TIME1=TIME 


\ 
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OEG=FREQ*TIME 

I0EG-DEG 

t)FG-(OFG-TOEG)*360. 

C 

IAREF=IMAX*SINCOMEGA«TIME) 

IBREF=IMAX»SINC OMEGA*TIME-PAI/2 . I 

IF( IAREF.GT.0.0 ) IAREF=IMAX 
IF( IAREF. LT.0.00 ) IAREF=-IMAX 
IF( IBREF.GT.0.0 ) IBREF=IMAX 
IFCIBREF. LT.0.00) IBREF=-IMAX 

DIA=YC 1 l-IAREF 
0IB=YC2)-IBREF 
C 

VOS=E 

VQS=E 

IF(IAREF.LT.O) VDS=-E 
IF( IBREF.LT. 0) VQS=-E 
C 

IF(TIME.LT. 0.02000) GOTO 6 
C IFCTIME.GT. 0.00000) GOTO 23 
ISW=0 

IFCOIA.LT. Dll )ISW=1 
IFCDIA.GT.DI) DI1=DIA 

IFCOIA.GT.DI. AND. IAREF. GT. 0.0. AND. ISW.EQ.O) GOTO 24 
IF( DIA.GT.DI ) IAE-1 
GOTO 27 
24 IAE=1 

IFtTIME.LT. 0.24375) GOTO 27 
SU=SW*1 

PSW=ABSC C 240. »YC 1 )»5.5E-06 >/<6»6.25E-03 > l+PSW 
27 IFCOIA.LT. C-DD) IAE=0 

IFC VDS.GT . 0 . 0 .AND . IAE . EQ 1) GOTO 18 
IFCVDS.LT. 0.0. AND. IAE.EQ.0) GOTO 18 
GOTO 22 
18 VOS=-VDS 

IFCTIME.LT. 0.24375) GOTO 22 
SWITCH-SWITCH*1 . 

22 IFCOIB GT.DI) IBE=1 

ifcoib.lt. c-Din ibe=o 

IFC VQS.GT.0.0. AND. IBE .EQ. 1 ) GOTO 19 
IFCVQS.LT. 0.0. AND. IBE.EQ.O) GOTO 19 
GOTO 23 
19 VQS=-VQS 

23 IF! TIME. GT. 0.231 25. AND. TIME. LT. 0.23166 12 )AREAIN=Y(6 ) 

IF CTIHE.GT. 0.231 25. AND.TIME.LT. 0.2 3166 12 )ARA0UT=YC5) 

IFC TIME. GT.0.2312S. AND. TIME. LT.O. 2316612 )ARMS=YC 7) 

C 

IFCTIME.LT. 0.24375) GOTO 8 
T0TSW=T0TSW*1. 

C 

8 CONTINUE 
C 

XC 1 )=FACT1«V0S-FACT2*YC 1 )+FACT3*YC 2 )*FACT4«YC4 )*FACT5»YC 3) 
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X( 2 )=FACTl*Vqs-FACT2»YC 2 )-FACT3*Yt 1 )-FACT4»YC 3 )*FACT5»YC 4 ) 

XC 3) a -FACT6*VDS»FACT55»Yt 1 )-FACT44*Y< 2 l-FACT7*YC4 )-FACT22*Y( 3 ) 
X(4) a -FACT6«VqS*FACT55*YC2)+FACT44*YCl)+FACT7*Y(3)-FACT22i*Y<4) 
X( 5 ) a LM*( YC 2 )*YC 31-YC 1 )*YC4) )*OMEGR/. 01875 
XC 6 ) = C YC 1 )*V0S*Y( 2 )»Vq3 )/. 01875 
X(7)sY(ll*Y(l)/. 01875 
C 

GOTO 2 

5 IPR=0 
NC-NC-1 
IFCNC >9,9,10 

9 NC=NCC 
IPR=1 

10 IF( TIME .GE . FINTIM ) IST0P=1 
IF t ( I PR* ISTOP ) . EQ. 0 ) GOTO 11 

C 

IFCTIME.LT. 0.24375) GOTO 11 
IOS=YU) 

IQS=Y(2) 

I0R=YC3) 

IQR=Y(4) 

C 

IMAGD-IDS+ICR 

IHAGQ=IQS*IQR 

C 

TQ=LM<*( iqS*IDR-IOS»iqR ) 

POUT=Tq»OMEGR 
PIN a I0S«VDS*iqS*VqS 
SUM1=SUM1*PIN 
SUM=SUM*POUT 
NUMBER=NUMBER*1. 

IF(TIME.LT. 0.24375) GOTO 31 
IF(DEG.GT.180. ) GOTO 41 

sumi=sumi*yc i ) 

IN=IN*1 

F1C INI a CMPLXC YC 1 ) iO.O ) 

F2( IN) a CMPLX( POUT >0.0 ) 

COUTINUE 

URITEC6,250)0EG, IDS, VOS, POUT 
IF CTIME.LT. 0.24400) POUT=0.0 

CALL FPLOTC 800 , IPNT1 ,AR1 , LR1 , ISTOP, 1 . 1 .TIME , IDS) 

CALL FPLOTC 800 , IPNT2 , AR2 , LR2 , ISTOP ,1,1 .TIME , VOS ) 

CALL FPLOTC 800 , IPNT4 , AR4 , LR4 , ISTOP , 1 . 2 .TIME , I0S ) 

CALL FPLOTC 800, IPNT4.AR4.LR4, ISTOP, 2, 2, TIME, VOS) 

CALL FPLOTC 800 , IPNT3.AR3 , LR3 , ISTOP ,1,2, TIME , IMAGO ) 

CALL FPLOTC 800, IPNT3.AR3, LP3, ISTOP, 2, 2, TIME, IDR) 

CALL FPLOTC 800 , IPNT5 , AR5 , LR5 , ISTOP, 1,1, TIME , POUT ) 

C 

11 INIT=0 

IFC ISTOP )1 ,1,6 

6 CONTINUE 

250 FORMATC 1 • ,F6. 1 ,3C F9.2 ,2X ) ) 

C AVER=SUM/NUMBER 

AVERI=2«SUMI/IN 
C EFF=AVER/AVER1 
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PIN=Y(6 J-AREAIN 
POUT=Y( 5)-ARA0UT 
RMS=(Y(7)-ARMSI«»0.5 

C *H**»)*«»*#)**»ti*»HM**«»«#**HHHM»**#**«IHHHHMH»***«l»»««***IHHMHH**«IHHI*** 

C NOTE .7 IS THE VOLTAGE DROP ACROSS THE DIOOE 

C 1.85 IS THE CE VOLTAGE OF THE TRANSISTOR 

C SWITCH IS THE • OF INTERVAL OF THE SWITCHING (DIOOE CONDUCTION! 

C TOTSW IS TOTOAL # OF INTEGRATION INTERVALS 

C POIOOE IS THE LOSS OF OIODE OURING CONDUCTION 

C PTRANS IS THE ON TIME LOSS OF THE TRANSISTOR 

C PSW IS THE SWITCHING LOSS OF EACH TRANSISTOR (E*I*TSW/T) 

C TOTSL IS THE TOTAL LOSS OF THE SWITCH 

PDIOOE=( SWITCH/TOTSW )»0.7*IMAX/PAt 
PTRANS=( 1 -SWITCH/TOTSW 1*1 .85*IMAX/PAI 
TOTSL=POIODE+PTRANS*PSW 
T0T4SL=4«T0TSL 
DTRIO=l-( SWITCH/TOTSW) 

C EFFIO=PQUT/( T0T4SL+PIN) 

EFFIO=POUT/PIN 

SW=SU+SW 

C 

WRITE ( 6,702 )POUT »EFFIO,SLIP> IMAX.OI >RMS, AVERI ,DTRIO 
WRITE ( 6 , 703 )SW , PSW , PTRANS , FOIODE . T0T4SL , IN 

702 FORMAT! ///• ' POKER OUT=* ,F6.0,4X» 'EFFIO=' ,F5.3,2X, 'SLIP 5 ' , 
IF7.5.2X, 'IMAX=' ,F4.0,2X, 'DI= * , F3. 0 ,4X , 

2 ' I! RH5 ) = ' , F5. 1 , 3X, ' I! AVE ) = ' >F5. 1 > 3X» 'DTRIO=' »F5. 3 »/////) 

703 FORMAT! /» 1 #SW/CYC= ‘ >F5. 1 ,8X, ' PSW= ' > F6 . 1 > 3X > ' PTRANS 3 ' ,F5.l, 

1 2X , ' PD IOOE= ' , F4 . 1 , 2X , ' T0T4SL= • , F6 . 1 , 5X , • IN= ' , 15 ,/////// ) 

CC WRITE! 6.705 ) 

C705 FORMAT! 1X,'Z'«11X,'I HARMONIC '.ISX.'P HARMONIC',//) 

C706 FORMAT! IX, F5. 1 ,5X ,F12 .3.5X, F15.8 ) 

CC ISTOP=0.0 

CC CALL FFT! FI ,8 ) 

CC CALL FFT(F2,3) 

CC 00 401 1=1,128 

CC MAGI! I )=CABS(F1( I))/256. 

CC MAG2! I )=CABS( F2( I ) 1/256 . 

CC Z=FL0AT(I-1) 

CC IF (I.EQ.l) MAGI ( 21=200. 

CC AMP1=MAG1( I )/MAGl( 2 ) 

CC AHP2=MAG2( I )/MAG2( 1 ) 

CC WRITE! 6 , 706 ) Z.AMP1.AMP2 

CC IF! I. EQ. 128) ISTOP=l 

C IF! AMP2 .GT. 300 ) AMP2=0.0 

C 

CC CALL FPLOT! 800 >IPNT4, AR4, LR4,IST0P, 1 , 1 ,Z, AMP1 ) 

CC CALL FPLOT! 800 , IPNT3, AR3, LR3, ISTOP, 1 , 1 , Z, AMP2 ) 

C401 CONTINUE 

STOP 
END 

SUBROUTINE FFT(X.M) 

COMPLEX X! 256 ) ,U,W,T 

N=2*«M 

NV2=N/2 
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NM1=N-1 

J=l 

DO 7 1=1, NMI 
IF(I.GT.J) GOTO 5 
T=X( J) 

X« J)=X(I1 

xiiir 

5 K=NV2 

6 IFCK.GE. J )GO TO 7 
J=J-K 

K=K/2 
GO TO & 

7 J=J+K 

PI=3. 14159265358979 

00 20 L=1,M 

LE=2*«L 

LEl=LE/2 

U=< 1.0,0.01 

Z=C0S( PI/FLOATI LEI ) ) 

T=-SIN( PI/FLOAT( LEI )) 

W=CMPLX( Z,Y) 

DO 20 J=I,LE1 
00 10 I=J,N,LE 
IP=I*LE1 
T=X( IP )*U 
X(IP»=X1IJ-T 

io xm=x(u*T 

20 U=U»W 

RETURN 
END 

SENTRY 

l A CURRENT OF PHASE A OF STATOR 

1 P POWER VS TIME 

1 I CURRENT HARMONIC MAGNITUDE VS HARMONIC ORDER 

1 P POWER HARMONIC MAGNITUDE VS HARMONIC OROER 

1 V VOLTAGE OF PHASE A OF STATOR 

♦STOP 
/* 

// 
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STATOR 


'OLTAGE 


(80- 


SPEED, 12 KW) 


25C-J 




Du ci uj: 
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c 


STATOR CURRENT (80S 


SPEED, 12 K1V) 



5C 


t /s 

iv r w v 


3CC 
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OUTPUT POWER (80T SPEED, 12K1Y) 





184 


CURRENT SPECTRUM (80% SPEED, 12 KW) 



0 3 6 9 12 15 19 21 24 27 30 33 35 33 42 45 ,3 51 

HfiP^ON [ C 
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POWER SPECTRUM (80°^ SPEED, 12KW) 



HPRMIJS ' 



FILE: PLOTP160 HATFIV 


A 


UNIVERSITY OF MISSOURI COMPUTER NETWORK 


//PLOTHS JOB (XXXXLHJi'JPL PROJS 
// PASSWORD? 

/•JOBPARM R=900,B=99000 
// EXEC PLOTSAS 
//SYSIN 00 * 

GOPTIONS 0EVICE=CA11077» 
DATA ONE I 
INPUT X Y1-Y3I 

CAROS I 


OATA FOR SOX SPEED 60 KU OPERATION 


I DEGREE I 
I I 
1.1 
1.6 
2.2 
2.8 

3.6 

4.7 

6.4 

8.7 

11.9 

16.5 

23.0 

23.4 

23.8 

24.2 

24.6 

25.2 

26.0 

27.1 

28.8 

29.2 

29.6 

30.0 

30.4 

30.9 

31.8 

32.9 

34.5 

36.8 

40.1 

40.5 

40.9 

41.3 

41.7 

42.3 

43.1 

44.2 

45.8 
48.1 


I 

CURRENT I 

I 

-22.81 

-20.59 

-13.39 

-16.23 

-13.22 

-9.07 

-3.40 

4.21 

14.19 

26.83 
41.97 

42.81 
43.65 
44.48 

45.29 

44.42 
47.99 

50.13 
53.03 

53.73 

54.42 

55.10 
55.77 

56.71 

58.00 
59.76 

62.14 

65.29 
69.33 

69.81 
70.26 

70.74 

71.20 

71.83 

72.71 
73.91 
75.54 

77.71 


I 

VOLTAGE I 

I 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 

240.00 


I 

POWER I 

I 

5920.99 

6231.94 

6548.00 

6363.89 

7330.49 

7997.52 

8965.61 

10373.63 

12412.63 
15306.36 
19177.87 

18808.73 
18330.34 

17862.74 

17405.93 
16778.39 

15927.89 

14799.00 
13350.31 

13531.00 

13812.66 

14090.93 

14365.63 

14747.93 

15275.69 
15995.06 

16954.89 

18187.43 

19653.47 

19466.14 

19217.69 

18979.47 

18750.66 

18443.44 

18042.04 

17539.42 

16957.45 

16387.70 
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48.5 

78.08 

240.00 

16533.59 

48.9 

78.45 

240.00 

16717.02 

49.3 

78.81 

240.00 

16893.27 

49.8 

79.17 

240.00 

17062.34 

50.3 

79.67 

240.00 

17288.99 

51.1 

80.36 

240.00 

17584.35 

52.3 

81.33 

240.00 

17950.66 

53.9 

82.65 

240.00 

18363.02 

56.2 

84.45 

240.00 

18727.66 

56.6 

84.76 

240.00 

18682.68 

57.0 

85.07 

240.00 

18627.90 

57.4 

85.38 

240.00 

18581.06 

57.8 

85.69 

240.00 

18542.06 

58.4 

86.12 

240.00 

18500.16 

59.2 

86.73 

240.00 

18467.24 

60.3 

87.58 

240.00 

18472.39 

62.0 

88.78 

240.00 

18530.52 

64.3 

90.48 

240.00 

18928.00 

64.7 

90.79 

240.00 

18958.41 

65.1 

91.09 

240.00 

18969 84 

65.5 

91.39 

240.00 

18971.66 

65.9 

91.70 

240.00 

18963.81 

66.5 

92.13 

240.00 

18936.08 

67.3 

92.75 

240.00 

18863.47 

68.4 

93.63 

240.00 

18693.46 

70.0 

94.91 

240.00 

16316.52 

70.4 

95.23 

240.00 

18347.82 

70.8 

95.56 

240.00 

13413.83 

71.2 

95.89 

240.00 

18485.64 

71.7 

96.22 

240.00 

18563.12 

72.2 

96.69 

240.00 

18682.33 

73.0 

97.38 

240.00 

18869.82 

74.2 

98.36 

240.00 

19171.70 

75.8 

99.80 

240.00 

19668.84 

76.4 

93.83 

-240.00 

18545.96 

76.9 

86.57 

-240.00 

17178.11 

77.5 

79.32 

-240.00 

15829.75 

78.1 

72.09 

-240.00 

14501.10 

78.9 

61.89 

-240.00 

12656.22 

80.1 

47.51 

-240.00 

10116.46 

81.7 

27.30 

-240.00 

6666.26 

82.2 

26.69 

240.00 

6761.41 

82.8 

27.39 

240.00 

7121.00 

83.4 

28.10 

240.00 

7488.83 

84.0 

28.83 

240.00 

7864.62 

84.8 

29.87 

240.00 

8409.14 

85.2 

30.41 

240.00 

8298.47 

85.6 

30.95 

240.00 

8100.25 

86.0 

31.49 

240.00 

7890.95 

86.4 

32.04 

240.00 

7670.58 

87.0 

32.84 

240.00 

7339.98 

87.8 

33.98 

240.00 

6834.56 

88.9 

35.65 

240.00 

6043.88 

90.5 

38.12 

240.00 

4773.80 

91.0 

38.76 

240.00 

4903.55 
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91.4 

39.40 

240.00 

5129.91 

91.8 

40.05 

240.00 

5358.91 

92.2 

40.71 

240.00 

5590.45 

’2 7 

41 66 

240.00 

5922.09 

93.6 

43.03 

240.00 

6398.95 

94.7 

45.03 

240.00 

7087.78 

96.3 

47.97 

240.00 

8087.09 

98.6 

52.39 

240.00 

9539.93 

101.9 

59.18 

240.00 

11643.59 

106 4 

69.95 

240.00 

14630.51 

112.9 

87.72 

240.00 

18636.06 

113.5 

82.94 

-240.00 

18127 67 

114.1 

76.89 

-240 00 

17470.94 

114.7 

70.88 

-240.00 

16835.94 

115.2 

64.89 

-240.00 

16222.70 

116.0 

56.49 

-240.00 

15392.52 

117.2 

44.71 

-240 00 

14292.68 

118.8 

28.28 

-240.00 

12835.57 

119.4 

29.04 

240.00 

13154.08 

120.0 

31.12 

240 00 

13562.98 

120.5 

33.23 

240.00 

13964 82 

121.1 

35.37 

240.00 

14359.29 

121.9 

38.44 

240.00 

14904.15 

123.1 

42.88 

240.00 

15647.33 

124.7 

49.35 

240.00 

16640.20 

127.0 

58.90 

240.00 

17917.97 

130.2 

73.20 

240.00 

19445. 14 

130.8 

69.33 

-240.00 

19197.31 

131.4 

64.19 

-240.00 

18870.05 

131.9 

59.09 

-240.00 

18562.40 

132.5 

54.02 

-240.00 

18273.94 

133.3 

46.92 

-240.00 

17898.85 

134.5 

37.00 

-240.00 

17433.36 

136.1 

23.22 

-240 00 

16903 72 

136.7 

24.93 

240.00 

17099.50 

137.2 

27.96 

240.00 

17353.89 

137.8 

31.02 

240.00 

17594.44 

138.4 

34.11 

240.00 

17321.02 

139.2 

38.54 

240.00 

18117.16 

140.3 

44.90 

240.00 

18486.34 

142.0 

54.12 

240.00 

18906.48 

144.3 

67.58 

240.00 

19289.04 

144.8 

64.52 

-240.00 

19189.06 

145.4 

60.20 

-240.00 

19071.40 

146.0 

55.92 

-240.00 

18969.96 

146.5 

51.68 

-240.00 

18334.70 

147.4 

45.74 

-240.00 

18791.24 

148.5 

37.46 

-240.00 

18712.46 

150.1 

26.02 

-240.00 

18705.63 

152.4 

10.35 

-240.00 

18893.12 

155.7 

-10.80 

-240.00 

19553.66 

156.2 

-7.91 

240.00 

19597.46 

156.8 

-3.69 

240.00 

19600.96 

157.4 

0.55 

240.00 

19584.65 

158.0 

4.83 

240.00 

19548.38 
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G.5 TWO-PHASE SIMULATION RESULTS FOR POWER REGENERATION AT 80% SPEED 


Plots 

Stator Voltage and Current 
Output Power 

Listing of Plot Data 


Comment: Source Program is the same as for G.3 except the slip frequency 


is negative. 
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//REGP JOB ( XXXXLH ) > ' GOOD ARZI ' > 

// PASSWORD=XXXX 
/•JOBPARM R=900,B=99000 
// EXEC PLOTSAS 
//SYSIN DD * 

GOPTIONS DEVICE=VARIAN; 

DATA ONE; 

INPUT X Y1-Y3; 

******* X*************** ****** (I**************************************** 
* 

« REGENERATION MODE 3 80Z SPEED JPL PROJECT 

* PLOT DATA 

•III 

• DEGREE I VOLTAGE I CURRENT I POWER 

********* | ******** | ********* | is*************************************** 

caros; 


0.5 

-234.49 

240.00 

-79927.88 

1.1 

-232.59 

240.00 

-79913.88 

1.6 

-230.73 

240.00 

-79883.56 

2.2 

-228.91 

240.00 

-79837.13 

3.0 

-226.41 

240.00 

-79744.81 

4.2 

-223.00 

240.00 

-79562.31 

5.8 

-218.45 

240.00 

-79205.38 

8.1 

-212.52 

240.00 

-78515.63 

11.3 

-205.15 

240.00 

-77208.75 

16.0 

-196.63 

240.00 

-74807.44 

22.5 

-188.09 

240.00 

-70622.75 

31.7 

-182.09 

240.00 

-64028.68 

40.9 

-181.74 

240 00 

-57917.55 

50.1 

-185.33 

240.00 

-53617.25 

59.3 

-191.03 

240.00 

-52301.88 

68.5 

-197.03 

240.00 

-54937.41 

77.7 

-201.52 

240.00 

-62237.88 

86.9 

-202.77 

240.00 

-74629.75 

89.2 

-202.38 

240.00 

-78543.88 

89.8 

-202.24 

240.00 

-79573.25 

90.2 

-202.12 

240.00 

-79703.56 

90.6 

-201.99 

240.00 

-79704.81 

91.0 

-201.85 

240.00 

-79698.00 

91.4 

-201.71 

240.00 

-79682.94 

92.0 

-201.48 

240.00 

-79647.94 

92.8 

-201.11 

240.00 

-79571.56 

94.0 

-200.52 

240.00 

-79411.38 

95.6 

-199.52 

240.00 

-79085.44 

97.9 

-197.77 

240.00 

-78438.06 

101.2 

-194.58 

240.00 

-77188.06 

105.8 

-188.51 

240.00 

-74862.19 

112.3 

-176.55 

240.00 

-70770.81 

121.5 

-152.17 

240.00 

-64283.78 

134.5 

•101.16 

240.00 

-56195.74 

147.5 

-29.33 

240.00 

-52653.23 

160.6 

63.87 

240.00 

-56617.87 

173.6 

177.76 

240.00 

-70080.00 

178.2 

222.68 

240.00 

-77304.44 

179.8 

239.11 

240.00 

-80170.69 
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180.4 

238.47 

-240.00 

-80355.13 

181.0 

236.56 

-240.00 

-80352.50 

181.5 

234.69 

-240.00 

-80333.38 

182.1 

232 87 

-240.00 

-80298.25 

182.9 

230.35 

-240.00 

-80221.63 

184.1 

226.93 

-240.00 

-80061.06 

185.7 

222.35 

-240.00 

-79734.38 

188.0 

216.40 

-240.00 

-79086.13 

191.3 

208.99 

-240.00 

-77834.75 

195.9 

200.42 

-240.00 

-75505.56 

202.4 

191.81 

-240.00 

-71410.00 

211.6 

185.73 

-240.00 

-64915.04 

220.8 

185.32 

-240.00 

-56869.37 

230.0 

188.85 

-240.00 

-54598.85 

239.2 

194.52 

-240.00 

-53279.17 

248.4 

200.51 

-240.00 

-55878.30 

257.6 

205.00 

-240.00 

-63112.98 

266.9 

206.26 

-240.00 

-75411.81 

269.2 

205.88 

-240.00 

-79299.13 

270.0 

205.67 

-240.00 

-80751 .25 

270.4 

205.55 

-240.00 

-80879.50 

270.8 

205.42 

-240.00 

-80878.31 

271.2 

205.27 

-240.00 

-80868.75 

271.6 

205.12 

-240.00 

-80851.06 

272.2 

204.88 

-240.00 

-80812.31 

273.0 

204.51 

-240.00 

-80730.44 

274.1 

203.90 

-240.00 

-80562.56 

275.8 

202.88 

-240.00 

-80225.38 

278.1 

201.10 

-240.00 

-79561.69 

281.3 

197.86 

-240.00 

-78288.13 

285.9 

191.72 

-240.00 

-75926.63 

292.4 

179.65 

-240.00 

-71782.19 

301.6 

155.10 

-240.00 

-65216.10 

314.7 

103.83 

-240.00 

-57010.73 

327.7 

31.73 

-240.00 

-53344.49 

340.7 

-61.74 

-240.00 

-57181.73 

353.7 

-175.89 

-240.00 

-70511.56 

358.3 

-220.89 

-240.00 

-77688.31 

360.0 

-237.35 

-240.00 

-80537.50 


PROC GPLOTJ 

PLOT Y1*X Y2#X /OVERLAY; 

SYMBOLI I=J0IN L= 1 ( 

SYMBOL2 I=J0IN L=13l 
LABEL X-DEGREE 

YI=CURRENT t VOLTAGE I 

TITLE VOLTAGE t CURRENT <80Z SPEED [REGENERATIVE ) ! 
PROC GPLOT; 

PLOT Y3»X /VAXIS=0 TO -85000 BY -50001 
SYMBOLI I=SPLINES 
LABEL X=OEGREE 

Y3= POWER | 

TITLE REGENERATIVE POWER (80X SPEED J» 



G.6 THREE-PHASE SIMULATION RESULTS 


Source Program 
Plots 

Stator Voltage and Current 
Output Power 
Listing of Plot Data 



195 


FILE: JPL330 WATFIV 


UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PA6E 001 


//B6330 JOS !XXXXEE), 'ABBAS GOOOARZI* >CLASSsA,MSGLEVEL=!0,0)> JPLOOOIO 

// PASSWORD =XXXX JPL00020 

/*J08PARM T=2 ,L=2 JPL00030 

/*ROUTE PRINT UMVMA.EEPR3 JPLOOOIO 

// EXEC WATFIV JPL00050 

//GO.SYSIH DO « JPL00060 

$ JOB JPL00070 

DIMENSION ARH800 ) , AR2! 800 ) , AR3! 800 ) > Y! 20 ) >X< 20 ) ,GA! 20 ) >PE! 20 ) > JPL00080 

1XX!4,20>,H<9),ZZ1!20>,ZZ2I20),ZZ3120),IFL(7),AR4<800> JPL00090 

LOGICAL* 1 LRl! 800 ) > LR2! 800 ) , LR3! 800 ) , LR4( 800 ) JPL00X00 

REAL LI 6,6 ),LR( 6,6), 01 6. 6). LI! 6,6), WKAREA( 54). LOS, LOR, LM, THETA, C.JPL00110 

lIAS,IBS,ICS>IAR,IBR,ICR>IER.RS,RR,DD,EE,FF,TQiAA,BB.CC>E,OMEGR, JPL00120 
2QL(6,6),IAREF,IBREF,ICREF,IMAX JPL00130 

C/ PASSWORO=XXXX,MSGCLASS=S JPL00140 

C*JOBPARM T=3 > L=3 JPL00150 

C*ROUTE PRINT UMMVSA.RO JPL00160 

C 50 FORMAT! />2X, 'TIME' i8X> 'VA' ,8X, 'IAS' >8X. 'IBS' ,8X, 'ICS 1 , 8X» JPL00170 

C 1 *IAR',8X,'IBR',8X, 'ICR', 5X, 'POWER', 7X, 'ANGLE* > JPL00180 

IMAX=100. JPL00190 

SLIP-0.0125 JPL00200 

DI=10. JPL00210 

E=240. JPL00220 

NRUN=1 JPL00230 

ISTOP=0 JPL00240 

KEEP=1 JPL00250 

IFL!1)=-1 JPL00260 

NINT=9 JPL00270 

IPNT1=0 JPL00280 

IPNT2=0 JPL00290 

IPNT3=0 JPL00300 

IPHTA=0 JPL00310 

C************************************************************ JPL00320 

C MOTOR PARAMETERS JPL00330 

C FOR BRIDGE INVERTER (JPL) DATA SIMPLIFIED MODEL JPL00340 

C*******************************************)***************** JPL00350 

L0S=75.36OE-06 JPL00360 

LOR=150.E-06 JPL00370 

LM=4200.E-06 JPL00380 

RS=0. 007345 JPL00390 

RR=0. 009813 JPL00400 

PAI=3. 14159 JPL00410 

PI2=2.*PAI/3. JPL00420 

THETA=0. JPL00430 

FREQ=330. JPL00440 

OMEGA=2*PAI*FREO JPL00450 

OMEGR=OMEGA*l 1 .-SLIP) JPL00460 

WRITE! 6 >64 ) JPL00470 

WRITE(6,66)L0S,L0R JPL00480 

WRITE(6>69)RS>RR JPL00490 

WRITE1 6 >68 )LM JPL00500 

WPITE(6>71 ISLIP JPL00510 

WRITE! 6 , 72 1FREQ JPL00520 

WRITE! 6 ,64 ) JPL00530 

C WRITE! 6>50 ) JPL00540 

64 FORMAT!' ' ,///,20X, '******MOTOR PARAMETERS******' ,///) JPL00550 
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66 FORMATS • ,20X, ’L0S=' ,F9.7,5X, ’LORs* ,F9.7> 

68 FORMATC • ,20X. , LM=' ,F9.7> 

69 FORMATC ' ,20X, 'RS S ' «F8.6,5X, *RR=' ,F8.6) 

71 FORMATC * >20X, 'SLIP 8 * «F6.4> 

72 FORMATC 1 ,20X, 'FREQUENCYs' ,F6. 1) 

WS=LM*L0S 

WR=LM+L0R 

T=SQRT(3. )/2.*OMEGR#LM 
Z=-0.5*LM 

U=SQRT( 3. )/3.*0ME6R*( 3./2.*LM*L0R ) 

DO 100 1=1,6 
DO 100 J=1 ,6 
L(I,J)=Z 
LR( I , J )-0. 

D(I,J)=0. 

100 CONTINUE 

DO 200 1=1,3 
L(I.I)=US 
LR( 1 ,1 )=RS 

200 CONTINUE 

00 201 1=4,6 
L( 1 , 1 J=WR 
LR( 1 ,1 1=RR 

201 CONTINUE 
L( 1 ,4 ) = LH 
L( 2,5 ) = LM 
LC 3,6 ) = LM 
L(4,l )=LM 
L(5,2)=LM 
L(6,3) = LM 
LR(4,2)=T 
LR(4,3)=-T 
LR ( 4 , 5 ) =U 
LR(4,6)=-U 
LR(5,1)=-T 
LR(5»3)=T 
LR(5»4)=-U 
LR(5,6)=U 
LR(6,1 )=T 
LR(6,2)=-T 
LR(6,4)=U 
LR(6,51=-U 
DO 400 1=1,6 
DO 400 J=1 ,6 

400 LI(I,J ) = L( I, J) 

CALL INVERS (LI,6) 

CALL MATMUL ( LI, LR, 6, 6, 6,0 ) 

1 IFL( 2 )=-l 

2 CALL RKF0UR(NINT,Y,X,6A,PE,XX,H,ZZ1 ,ZZ2 ,ZZ3, TIME, PNTS, KEEP, IFL) 
IF(IFL(4) 13,4,5 

3 INIT=1 
TIME1=-1 
NC=l 
NCC=3 
TIME=0. 


JPL00560 
JPL00S70 
JPL00580 
JPL00590 
JPL00600 
JPL00610 
JPL00620 
JPL00630 
JPL00640 
JPL00650 
JPL00660 
JPL00670 
JPL00680 
JPL00690 
JPL00700 
JPL00710 
JPL00720 
JPL00730 
JPL00740 
JPL00750 
JPL00760 
JPL00770 
JPL00780 
JPL00790 
JPL00800 
JPL00810 
JPL00820 
JPL00830 
JPL00840 
JPL00850 
JPL00860 
JPL00870 
JPL00880 
JPL00890 
JPL00900 
JPL00910 
JPL00920 
JPL00930 
JPL00940 
JPL00950 
JPL00960 
JPL00970 
JPL00980 
JPL00990 
JPL01000 
JPL01010 
JPL01020 
JPL01030 
JPL01040 
JPL01050 
JPL01060 
JPL01070 
JPL01080 
JPL01090 
JPL01 1 00 


i 
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PAGE 003 



FINTIH=0. 060505 


JPL01110 


Y(1)=0. 


JPL01120 


Y( 2 ) = 0 . 


JPL01130 


Y(31=0. 


JPL01190 


Y(4I=0. 


JPL01I50 


Y(5)=0. 


JPL01I60 


Y(6)=0. 


JPL01170 


Y(7) = 0.0 


JPL01180 


Y<8)=0 0 


JPL01190 


Y(9)=0.0 


JPL01200 


H(l)=0.lE-4 


JPL0I210 


E2-E/2. 


JPL01220 


IAE=-1 


JPL01230 


IBE=-1 


JPL01240 


ICE=-1 


JPL01250 


UPIN=0.0 


JPL01260 


UPOUT-O . 0 


JPL01270 


UPLOSS=0.0 


JPL01280 


POWER=0.0 


JPL0I290 

c 

SUM=0.0 


JPL0I300 

c 

PIN-0.0 


JPL01310 

c 

RMSI=0.0 


JPL01320 

c 

RMS2=0.0 


JPL01330 

c 

COUUT^O . 0 


JPL01340 


SW=0.0 


JPL01350 


PSW=0.0 


JPL01360 


PERIOD=l ./FREQ 


JPL0I370 


TF=FINTIM-PERIOD 


JPL01380 


4 IF(TIME-TIME1 )7,8,7 


JPL01390 


7 TIME1=TIME 


JPL01400 

C******************************* ********************************* ****** 

JPL01910 

C****** ************ *************** ####**####### ******* XXX ******* XXX XXX* 

JPL01420 

c 



JPL01430 

c 

’ 


JPL01AA0 

c 

REFFERENCE SIGNAL FOR GENERATING 

INPUT VOLTAGE 

JPL01450 

c 

AND THE CURENT CONTROLER REFRENCE 

SIGNAL 

JPL01960 

c 



JPL01470 

C*X*** *********************** ************************* ******** ********* 

JPL01480 

c************* *************************************** ****************** 

JPL01A90 


DEG=FREQ»TIME 


JPL01500 


I0EG=DES 


JPL01510 


DEG=( DEG-IDEG )#360 . 


JPL01520 

c 



JPL01530 


VA= 1 60 . #SIN( OMEGA«TIME ) 


JPL01580 


VB=160.*SIN(OMEGA*TIME-PI2I 


JPL01590 


VC=160.*SIN(OMEGA*TIME*PI2> 


JPL01600 


IF( TIME. LT. 2. 8890909) GOTO 88 


JPL01650 


VA0=E2 


JPL01540 


VB0=E2 


JPL01550 


VC0=E2 


JPL01560 

c 



JPL01570 


IAREF=IMAX«SIN< OMEGA«TIME ) 


JPL015B0 


IBREF=IMAX*SIN< OMEGA«TIME-PI2 1 


JPL01590 


ICREF=IMAX»SIN<OMEGA*TIME+PI2 1 


JPL01600 

c 



JPL01610 
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FILE: JPL330 MATFIV 


UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PA6E 004 


C 

12 


13 


C 

87 


88 


IFdAREF.LT. 0.0) VA0=-E2 
IFCIBREF.LT. 0.0) VB0=-E2 

IFCICREF.LT. 0.0) VC0=-E2 

IFCTIME.LT. 2. 8890909) GOTO 87 

DIA=TC 1 )-IAREF 
DIB=YC 2 )-IBREF 
OIC=Y(3)-ICREF 

IFCDIA.GT.DI) IAE=1 

IF ( TIME . LT.TF . AND.DIA.GT.DI ) GOTO 13 

SW=SW+1 . 

PSU=ABSC C 240. *YC 1 )«S.5E-06 )/C6*PERIOO I HPSW 

IFCOIA.LT. C-DD) IAE=0 

IFC VAO.GT. 0.0. AMD.IAE .EQ. 1 ) VA0=-VA0 

IFC VAO. LT. 0.0. AMD . IAE . EQ. 0 ) VAO=-VAO 


IFCDIB.GT.DI) IBE=1 

IFCDIB.LT. (-DI)l IBE=0 
IFC VBO .GT.O.O. AND . IBE .EQ. 1 ) 
IFCVBO.LT. 0.0. AND IBE.EQ.O) 

IFCDIC.GT.OI) ICE=1 

IFCDIC.LT. (-DID ICE=0 
IFC VCO. GT.O.O. AND. ICE. EQ.l) 
IFC VCO . LT.0 . 0 .AND . ICE . EQ. 0 ) 


VBO=-VBO 

VB0=-VB0 


VCO=-VCO 

vco=-vco 


UPIN=YC7) 

UP0UT=Y(8) 

UPL0SS-YC9) 


VN=( VAO+VBO+VCO )/3. 

VA=VAO-VN 

VB=VBO-VN 

VC-VCO-VN 

CONTINUE 

IFC TIME. GT. 0.03318. AND. TIME. LT. 0.03323) 
IF ( TIME .GT. 0 .03318. AND. TIME . LT.0 03323) 
IFC TIME. GT. 0.03318. AND. TIME. LT. 0.03323) 

THETA=OMEGR*TIME 
DO =SINC THETA) 

EE=SINCTHETA+2.*PAI/3. ) 

FF=SIN( THETA-2. *PAI/3. ) 


IAR=2./3.*( ( COSC THETA ) 40.5 )»YC 4 )*CCOSC THETA-2 .*PAI/3. )*0.5 )*YC 5) 
. l+CCOSCTHETA+2.«PAI/3. H0.5)»Y<6)> 


8 CONTINUE 

XC 1 )=-0C 1 » 1 )»Yl 1 )-Dl 1 ,2 )#YC 2 )-0C 
1 -OC 1 ,6 )*YC6)+LI( 1 , 1 )#VA+LI( 

XC 2 )=— DC 2,1 )*Yl 1 )-DC 2,2 )»YC 2 )-DC 
1 -DC 2,6 )*YC6 )*LI( 2, 1 )*VA+LIC 
XC 3)=-DC 3,1 )*YC 1 )-D( 3,2 )*Y( 2 )-DC 
1 -DC 3,6 )»YC6 )*LI( 3, 1 )*VA*LIC 

XC 4)=-0C4, 1 )*YC 1 )-D( 4,2 )#YC 2 )-0( 
1 -DC4,6)*Y(6) + LIC4,1 )»VA*LIC 

XC5)=-0C5, 1 )»YC 1 )-0C5,2 )#YC 2 )-DC 


1 >3)»YC 3)-Dl 1 >4)»Y(4 )-D( 1 ,5 )#Y(5 ) 

1 . 2 )*VB*LI( 1 , 3)»VC 

2,3)«YC 3 )-D( 2 ,4 )»YC4 )-DC 2 >5 )»YC 5 ) 

2.2 ) »VB*LI( 2,3)»VC 

3.3) »YC 3 )-D( 3,4 )*YC4 )-Dt 3,5 )#YC5) 

3.2 )*VB+LI( 3,3)»VC 
4,3>*Y(3)-0C4,4)*YC4)-DC4,5)»YC5) 

4.2) «VB*LK4,3)«VC 

5. 3 ) #YC 3)-D( 5,4)»YC4 )-DC5,5)*YC5 ) 


JPL01620 

JPL01630 

JPL01640 

JPL01650 

JPL01660 

JPL01670 

JPL01680 

JPL01690 

JPL01700 

JPL01710 

JPL01720 

JPL01730 

JPL01740 

JPL01750 

JPL01760 

JPL01770 

JPL01780 

JPL01790 

JPL01800 

JPL01810 

JPL01820 

JPL01830 

JPL01840 

JPL01850 

JPL01660 

JPL01870 

JPL01880 

JPL0I890 

JPL01900 

JPL01910 

JPL01920 

JPL01930 

JPL01940 

JPL01950 

JPL01960 

JPL01970 

JPL01980 

JPL01990 

JPL02000 

JPL02010 

JPL02020 

JPL02030 

JPL02040 

JPL02050 

JPL02060 

JPL02070 

JPL02080 

JPL02090 

JPL02100 

JPL02110 

JPL02120 

JPL02130 

JPL02140 

JPL02150 

JPL02160 
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MATFIV A UNIVERSITY OF MISSOURI COMPUTER NETWORK 


PAGE 005 


it 

6 
250 
C 12 


1 -D(5.6)»Y(6)UI(5,1)»VA4LI(5.2)«VB«LI(5.3)»VC 

X(6)=-0(6,l )«YC n-D(6,2)»Y(2)-0(6.3)»Y(3)-0(6,4>»Y(4>-D(6,5>«Y(5) 
1 -0(6,6)»Y<6>4LIC6,1)«VA4LI(6,2)*VB4LK6,3)«VC 

X(7)*< VA*Y< 1 )+VB*Y( 2 )4VC»YI 3 ) )/. 02632 
X(8) = POUER/.02632 

X(9) = 3.*(Y(1)»Y( 1 )*RS*IAR*IAR»RR)/. 02632 

GOTO 2 
5 I PR-0 
NC-NC-1 
IF( NC )9,9, 10 
9 NC=NCC 
IPR=1 

10 IF( TIME .GE . FINTIM )IST0P= 1 
IFUIPR4ISTOPI.EQ.O) GOTO 11 


IAR=2./3.*( (C0S( THETA 140.5 )«YC 9 )4(C0S( THETA-2. »PAI/3. )40.5)«Y(5) 
l4(C0S(THETA42.«PAI/3. )40.5 )*Y<6>) 

IER=2./3.»( (C0S(THETA42.*PAI/3. U0.5 )«Y(<* |4( COSt THETA (40.5 )«Y(5) 
14( COS( THETA-2. »FAI/3. )40.5)*Y<6>) 

ICR=2./3.»t ( COS( THETA-2. 4PAI/3. )40 .5 )*Y<9 )4(C0S( THETA42.»PAI/3. ) 
14 0.5 )*Y( 5 )4(C0S( THETA )40.5)*Y(6)) 

TQ=< 00*Y< 1 )«IAR4FF»Y( 2 )*IAR4EE»Y( 3 )*IAR4EE*Y( 1 )*IBR 
1 4D0»Y( 2 )*IBR4FF*Y( 3 )#IBR+FF»Y( 1 )*ICR4EE»Y( 2 )»ICR 

1 4DD#Y(3)*ICR)«(-LM) 

POWER=TQ*OMEGR 

IFUIME.LT. 0.05756) GOTO 11 

IAS=Y(1) 

IBS=Y( 2 ) 

ICS=Y<3) 

IAR=Y<4) 

IBR=Y(5) 

ICR=Y(6) 

WRITE! 6 >250 )TIME >VA, Y( 1) > Y( 2 ) ,Y( 3 ) ,Y<4) , Y( 5 ) >Y(6 ) > POWER >0EG 

PIN=(VA»IAS4VB»IBS4VC*ICS)4PIN 

SUM-SUM4 POWER 

Rt1Sl=RMSl4(IAS«IAS) 

RHS2-RMS24( IAR*IAR ) 

COUNT^COUNT 4 1 . 

IF C TIME. LT. 0.057597) POWERPC. 0 

CALL FPLOTl 800 > IPNT1 , API , LR1 > ISTOP, 1 ,2 ,DEG, IAS ) 

CALL FPLOT(800>IPNT1 > AR 1 > LR1 > ISTOP>2 >2 >DEG> VA ) 

CALL FPLOT(800,IPNT2,AR2,LR2, ISTOP, 1,2, DEG, IBS) 

CALL FPLOT(800»IPNT2,AR2 , LR2 , ISTOP.2 ,2 ,DEG, VB ) 

CALL FPLOT(800,IPNT3,AR3,LR3, ISTOP, 1,2, DEG, ICS) 

CALL FPLOT( 80 0,IPNT3,AR3,LR3, ISTOP.2 ,2 >0EG , VC ) 

CALL FPLOT( 800, IPNT4, AR4 , LR4. ISTOP, 1,1, DEG, POWER) 

INIT=0 

IF( ISTOPll , 1,6 
CONTINUE 

FORMATC F8. 6. 8( F8. 0. 3X ), F5. 1 ) 

FORMAT! • ' ,A20,/,6(F12.6,3X)) 

POWERI=Y(7)-UPIN 
POWERO= Y( 8 ) -UPOUT 


JPL02170 
JPL02 1Q0 
JPL02190 
JPL02200 
JPL02210 
JPL02220 
JPL02230 
JPL02290 
JPL02250 
JPL02260 
JPL02270 
JPL02280 
JPL02290 
JPL02300 
JPL02310 
JPL02320 
JPL02330 
JPL02390 
JPL02350 
JPL02360 
JPL02370 
JPL02380 
JPL02390 
JPL02400 
JPL02910 
JPL02420 
JPL02430 
JPL02990 
JPL02950 
JPL02460 
JPL02970 
JPL02A80 
JPL02490 
JPL02500 
JPL02510 
JPL02520 
JPL02530 
JPL02540 
JPL02550 
JPL02560 
JPL02570 
JPL02580 
JPL02590 
JPL02600 
JPL026IO 
JPL02620 
JPL02630 
JPL02690 
JPL02650 
JPL02660 
JPL02670 
JPL02680 
JPL02690 
JPL02700 
JPL02710 
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pioss=rm-uPLoss 

EFFIO=POMERO/POHERI 
C AVER=SUM/COUNT 

C PIN=PIN/COUNT 

C RMS1 =! RMS 1/COUNT )»». 5 

C RMS2=!RMS2/COUNT)*».5 

C PLOST= ( RMSl #RMS 1 1»RS*(RMS2«RMS2 l»RR 

C EFF=AVER/PIN 

EFFL=!POWERI-PLOSS)/POWERI 

EFFPSW=I POWERI-PLOSS )/C P0WERI*PSW»31 

WRITE! 6, 266 lPOWERI.POWERO.EFFIO, SLIP, IMAX.EFFL 

266 FORMAT!' ',///, ’«POWER IN=' ,F8.1 , AX, 'POWER OUT=' ,F8. I, AX, 

1 ' EFF! 101= ' »F5. 3>AX» 'SLIP= ' ,F7.5,AX, ' IMAX= ' 

1 ,F5. 1 ,3X,F5. 3, ' 1 ' ) 

WRITE(6,267ISW,PSW,EFFPSW 

267 FORMAT!' ’ , 'SW= ' ,F6. 1 , ' PSW/LEG= ' ,F8. 1 , 'EFFPSW= ' ,F8.A,////1 
C 

STOP 

END 

C 

SUBROUTINE INVERS! ZZBUS.N ) 

DOUBLE PRECISION ZBUSI6.6) 

DIMENSION ZZBUS! 6 ,6 ) 

DO 1001 1=1 ,N 
00 1001 J=1 ,N 
ZBUS! I, JI=ZZBUS( I, J I 

1001 CONTINUE 

DO 1070 1=1, N 
TRY=ZBUS< 1,1) 

IF! ABSlTRYniOOO, 1090, 1000 
1000 00 10A0 J=1 ,N 

IF! I-J 11010 , 10AO ,1010 
1010 DO 1030 K=1,N 

IF(I-K) 1020,1030.1020 

1020 ZBUS! J,K )=ZBUSI J,K l-ZBUS! J,I )»ZBUS( I ,K l/ZBUS! 1,11 
1030 CONTINUE 
10A0 CONTINUE 

DO 1060 K=1,N 
IFIK-I110S0, 1060, 1050 
1050 ZBUS! K,1 1=ZBUS(K,I l/ZBUS! 1,1 1 
ZBUS! I , K 1=ZBUSI I . K l/ZBUS! 1 , 1 1 
1060 CONTINUE 

ZBUS! 1,1 )=-l .O/ZBUS! I«I 1 
1070 CONTINUE 

DO 1080 1=1, N 
DO 1080 J=1,N 
ZBUS! I, J I=-ZBUS( X, J I 
1080 CONTINUE 

DO 1002 1=1, N 
DO 1002 J=1,N 
ZZBUS! I , J )=ZBUS! I , J 1 

1002 CONTINUE 

DPT1DM 

1090 WRITEI6, 11001 

1100 FORMAT! ‘ERROR IN MATRIX INVERSE' 1 


JPL02720 
JPL02730 
JPL027AO 
JPL02750 
JPL02760 
JPL02770 
JPL02780 
JPL02790 
JPL02800 
JPL0281 0 
JPL02820 
JPL02830 
JPL028A0 
JPL02850 
JPL02860 
JPL02870 
JPL02880 
JPL02890 
JPL02900 
JPL02910 
JPL02920 
JPL02930 
JPL029A0 
JPL02950 
JPL02960 
JPL02970 
JPL02980 
JPL02990 
JPL03000 
JPL03010 
JPL03020 
JPL03030 
JPL030AO 
JPL03050 
JPL03060 
JPL03070 
JPL03080 
JPL03Q90 
JPL03100 
JPL031 10 
JPL03120 
JPL03130 
JPL031A0 
JPL03150 
JPL03160 
JPL03170 
JPL03180 
JPL03190 
JPL03200 
JPL03210 
JPL03220 
JPL03230 
JPL032A0 
JPL03250 
JPL03260 
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RETURN 

END 

C 

SUBROUTINE MATMUL ( A.B.M.l.N.C) 
DIMENSION A(M>L) iBt L.NI tC(M,N) 
DO 100 1=1, M 

DO 100 J=l,N 

C(I,J)=0. 

DO 100 K=1,L 

100 C(I,J)=C(I,JHA(I,K1»B(K,J) 
RETURN 
END 

tENTRY 

1 AV IAS t VA VS. TIME 

1 P POWER OF THE MOTOR 

iSTOP 
/» 

// 


JPL03270 

JPL03280 

JPL03290 

JPL03300 

JPL03310 

JPL03320 

JPL03330 

JPL03340 

JPL03350 

JPL03360 

JPL03370 

JPL033S0 

JPL03390 

JPL03400 

JPL03410 

JPL03920 

JPL03430 

JPL03940 



Hzm33cn moD-ino< 


VOLTAGE & CURRENT (80% SPEED, MOTOR OPERATION) 


ucc- 





3H1XOTI 
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OUTPUT POWER (80% SPEED MOTOR OPE 



— 1 — n 1 1 

0 25 50 75 IOC 125 ISO I ’S 2CG JJ: _;0 :l 
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//THREEP JOB ( XXXXLH ) > 'GOOOARZI ' > 

// PASSWORO=XXXX 
/•JOBPARM R=900,B=99000 
// EXEC PLOTSAS 
//SYSIN DO * 

GOPTIONS DEVICE 1 VAR I AN; 

DATA ONES 
INPUT X Y1-Y3; 

********** «*««»«*** ****** ******************************** ***** «««***** 


* 


« 


THREE PHASE PLOT OATA JPL PROJECT 


********************************************************************** 


CARDS; 


2 4 

-238.1 

80 0 

62992 7 

4 8 

-215.4 

80 0 

61387 7 

7 2 

-193 3 

80.0 

59997 3 

9.5 

-172 0 

80 0 

58817.8 

11.9 

-151.4 

80 0 

57843 1 

14.3 

-131.4 

80 0 

57066.1 

16.6 

-112.2 

80.0 

56478.4 

19.0 

-93.8 

80.0 

56071 7 

21.4 

-76.1 

80.0 

55836 0 

23.8 

-59.1 

80 0 

55761.7 

26.1 

-42.9 

80.0 

55837.9 

28.5 

-27.5 

80.0 

56052.8 

30.9 

-12 8 

80.0 

56394.6 

33.2 

1.1 

80 0 

56850.5 

35.6 

14.2 

80.0 

57407.2 

38.0 

26.5 

80.0 

58051.1 

40.4 

38.1 

80.0 

58768.1 

42.7 

48.9 

60.0 

59543.6 

45.1 

59.0 

80.0 

60362.7 

47.5 

68.3 

80.0 

61210.5 

49.8 

76.8 

80.0 

62070.9 

52.2 

84.6 

80.0 

62929.0 

54.6 

91.6 

80.0 

63768.9 

56.9 

97.9 

80 0 

64574.7 

59.3 

103.5 

80 0 

65330.1 

61.7 

120.9 

160.0 

63515.1 

64.1 

135.4 

160.0 

62315.5 

66.4 

151.3 

160.0 

60915.1 

68.8 

166.6 

160.0 

59724.7 

71.2 

181.2 

160.0 

58741.3 

73.5 

195.1 

160.0 

57956.0 

75.9 

208.5 

160.0 

57359.9 

78.3 

221.2 

160.0 

56947.8 

80.7 

233.4 

160.0 

56707.7 

83.0 

244.9 

160.0 

56629 4 

85.4 

256.0 

160.0 

56705.8 

87.8 

266.5 

160.0 

56921.1 

90.1 

276.6 

160.0 

57265.6 

92.5 

286.2 

160.0 

57726.8 

94.9 

295.3 

160.0 

58291.0 
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97.3 

304.1 

160.0 

58944.8 

99.6 

312.4 

160.0 

59674.0 

102.0 

320.5 

160.0 

60464.4 

104.9 

328.1 

160.0 

61300.7 

106.7 

335.5 

160.0 

62168.1 

109.1 

342.7 

160.0 

63051.2 

111.5 

349.6 

160.0 

63934.3 

113. 8 

356.2 

160.0 

64802 . 2 

116.2 

362.7 

160.0 

65638.1 

118.6 

369.1 

160.0 

66426.8 

121.0 

375.3 

80.0 

67152.6 

123.3 

367.0 

80.0 

65000.2 

125.7 

360.7 

80.0 

63487.0 

128.1 

354.4 

80.0 

62183.1 

130.4 

348.1 

80.0 

61083.0 

132.8 

341.9 

80.0 

60180.0 

135.2 

335.7 

80.0 

59467.0 

137.5 

329.7 

80.0 

58935.5 

139.9 

323.9 

80.0 

58577.0 

142.3 

318.2 

80.0 

58381.4 

144.7 

312.7 

80.0 

58338.2 

147.0 

307.5 

80 0 

58437.0 

149.4 

302.6 

80.0 

58665.8 

151.8 

297.9 

80.0 

59012.7 

154.1 

293.6 

80.0 

59464.2 

156.5 

289.6 

80 0 

60007.5 

158.9 

286.1 

80.0 

60627.9 

161 3 

282.9 

80.0 

61312.5 

163.6 

280.2 

80.0 

620 -* 6.2 

166.0 

277.9 

80.0 

62814 4 

168.4 

276.1 

80 0 

63602.2 

170.7 

274.8 

80.0 

64393.3 

173 1 

274.0 

80.0 

65172.8 

175.5 

273.8 

80.0 

65923.1 

177.9 

274.2 

80.0 

66630.7 

180.2 

275.1 

- 80.0 

67280.4 

182.6 

247.6 

- 80.0 

65003.4 

185.0 

225.0 

- 80.0 

63374.9 

187.3 

203.0 

- 80.0 

61957.1 

189.7 

181.7 

- 80.0 

60747.6 

192.1 

161.1 

- 80.0 

59738.1 

194.4 

141.2 

- 80.0 

58922.0 

196.8 

122.0 

- 80.0 

58291.2 

199.2 

103.6 

- 80.0 

57837.4 

201.6 

85.9 

- 80.0 

57551.2 

203.9 

69.0 

- 80.0 

57422.9 

206.3 

52.9 

- 80.0 

57441.2 

208.7 

37.5 

- 80.0 

57594.9 

211.0 

22.9 

- 80.0 

57871.9 

213.4 

9.0 

- 80.0 

58259.6 

215.8 

- 4.1 

- 80.0 

58744.8 

218.2 

- 16.4 

- 80.0 

59314.1 

220.5 

- 27.9 

- 80.0 

59953.4 

222.9 

- 38.7 

- 80.0 

60648.3 

225.3 

- 48.7 

- 80.0 

61384.1 
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227.6 

- 57.9 

- 60.0 

62165.7 

230.0 

- 66.6 

- 80.0 

62918.1 

232.6 

- 76.1 

- 80.0 

63685.0 

236.8 

- 81.2 

- 80.0 

66631.7 

237.1 

- 87.6 

- 80.0 

65162.2 

239.5 

- 93.0 

- 80.0 

65801.0 

261.9 

- 110.3 

-160 0 

63900.6 

266.2 

- 126.8 

-160 0 

62610.6 

266.6 

- 160.7 

-160 0 

61121.6 

269.0 

- 155.9 

-160 0 

59362 6 

251.3 

- 170.5 

-160 0 

58768.7 

253.7 

- 186.6 

- 160.0 

57891.7 

256.1 

- 197.7 

- 160.0 

57205 8 

258.5 

- 210.6 

- 160.0 

56703 2 

260.8 

- 222.6 

- 160.0 

56370.8 

263.2 

- 236.1 

- 160.0 

56202.0 

265.6 

- 265.1 

-160 0 

56186 1 

267.9 

- 255.7 

- 160.0 

56310.3 

270.3 

- 265.7 

- 160.0 

56563.8 

272.7 

- 275.3 

-160 0 

56936 6 

275.1 

- 296.6 

-160 0 

57609.3 

277.6 

- 293.2 

- 160.0 

57973.1 

279.8 

- 301.5 

- 160.0 

58615.2 

282.2 

- 309.5 

- 160.0 

59318 3 

286.5 

- 317.2 

- 160.0 

60068 8 

286.9 

- 326.6 

- 160.0 

60851.5 

289.3 

- 331.7 

- 160.0 

61651 6 

291.6 

- 338.6 

- 160.0 

62652 3 

296.0 

- 365.3 

- 160.0 

63260 0 

296.6 

- 351.8 

-160 0 

63997 6 

298.8 

- 358.2 

- 160.0 

66710.0 

301.1 

- 366.6 

- 80.0 

65360.8 

303.5 

- 356.1 

- 80.0 

63169 2 

305.9 

- 369.8 

- 80.0 

61577.6 

308.2 

- 363.5 

- 80.0 

60218 0 

310.6 

- 337.2 

-80 0 

59065.0 

313.0 

- 331.0 

-80 0 

58112 6 

315.6 

- 326.9 

- 80.0 

57352.7 

317.7 

- 318.9 

- 80.0 

56777.5 

320.1 

- 313.1 

- 80.0 

56377.9 

322.5 

- 307.6 

-80 0 

56165.2 

326.8 

- 302.0 

- 80.0 

56068.6 

327.2 

- 296.8 

- 80.0 

56135.7 

329.6 

- 291.9 

- 80.0 

56336.6 

331.9 

- 287.3 

- 80.0 

56657.9 

336.3 

- 283.0 

- 80.0 

57088.6 

336.7 

- 279.1 

- 80.0 

57613 6 

339.1 

- 275.5 

- 80.0 

58219.1 

361.6 

- 272.6 

- 80.0 

58892.3 

363.8 

- 269.7 

- 80.0 

59617.1 

366.2 

- 267.5 

- 80.0 

60378.2 

368.5 

- 265.7 

- 80.0 

61162.9 

350.9 

- 266.5 

- 80.0 

61956 6 

353.3 

- 263.8 

- 80.0 

62739.2 

355.7 

- 263.6 

- 80.0 

63699.7 
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358.0 -264.0 -80.0 64218.1 

! 

PROC GPLOT; 

PLOT Y1«X Y2*X /OVERLAY! 

SYMBOL1 1= JOIN L=li 
SYMBOL2 I=JOIN L=13l 
LABEL X=0EGREE 

Y1=V0 LTAGE t CURRENT i 

TITLE VOLTAGE t CURRENT (80/ SPEED, MOTOR OPERATION); 
PROC GPLOT! 

PLOT Y3*X/VREF=0; 

SYMBOL1 I=SPLINE! 

LABEL X=0EGREE 

Y3= POWER ! 

TITLE OUTPUT POWER (80/ SPEED MOTOR OPERATION), 

/* 

// 
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1.0 OBJECT 
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For the majority of practical ac machines, including ideal induction 
machines, a great simplification is obtained by expressing the equations 
in a new reference frame and introducing certain currents and voltages which 
are different from, but are related to, the actual ones. The currents can 
have a physical meaning in that they can be considered to flow in windings 
acting along two axes at right angles, called the airect and quadrature axes. 
In this paper a two-axis d-q theory of the induction motor in a stationary 
reference frame is developed. 

2.0 INDUCTION MOTOR EQUATIONS IN THREE PHASE VARIABLES 
It is assumed that the induction motor is ideal and symmetrical, which 
means : 

1 . Urn form air gap 

2. Linear magnetic circuit 

3. Identical stator windings which produce a sinusoidal MMF wave 
in space 

4. No changes in resistances due to temperature and frequency 
changes 

The stator circuits are oriented along the a s> b $ and c $ axes which are 

120° apart. The rotor circuits are along the a . b and c axes which are 
r 3 r r r 

also 120° apart. They rotate relative to the stationary stator axes with 

dQ 

an angular velocity ,jj r = ^ > which is the motor shaft velocity, and the angl 
9 is the angle between the stator and rotor. 

The basic KVL equations for the three phase induction motor are - for 


the stator 
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di dX 


1 V 

I v as 1 


'bs 


R,.. i + L, 


R l . iu. + L, 


V. ' = R__ t + L. 


and for the rotor 


V '■ R_ + Li 


I V t _ = R l .. i. + Li 


^cr ”i R cr ^cr + '"l cr dt + dt 


dt 

f dt 

di bs 

, dX bs 

dt 

dt 

di 

dX„ 

cs 

, cs 

dt 

dt 

di 

dX,„ 

ar 

i ar 

dt 

dt 

di br 

d: 'br 

dt 

dt 

di 

d> 

cr 

+ , cr 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


(2.4) 

(2.5) 

( 2 . 6 ) 


where 


V„ A instantaneous line to neutral voltage of stator Phase a. 

Phase b 

" 11 " Phase c 


as 
V bs ^ 


V cs^ 


II II II 


V ar^ 
^br — 
V cr^ 


II II II 


It tl II 


II II II 


" " rotor Phase a_ 

Phase j) 
Phase c 


II II II 


R a _ 4 line to neutral resistance of stator Phase a 

a b — 


R bs A " " 

R cs 4 " " 


Phase b 


Phase c 
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R A line to neutral resistance of rotor Phase a 
ar = — 

R. A " " " Phase b 

br = — 

R„ A " " " Phase c 

cr = - 


^as 

A 

instantaneous 

phase or 

line 

curren 

it in stator 

Phase £ 


ibs 

A 

II 


II li 

li 

ii 

II II 

Phase £ 


i 


II 


ii li 

II 

ii 

li il 

Phase £ 


cs 










i 

A 

ii 


•i ii 

li 

ii 

" rotor 

Phase £ 


ar 










1 br 

A 

n 


ii li 

it 

n 

n ii 

Phase £ 


i 

A 

M 


li li 

n 

H 

it ii 

Phase £ 


cr 










A 

as 

A 

instaneous 

flux linkage 

of 

stator 

Phase a_ with 

i air gap 

magnetic field 

X. 

A 

li 

II 

II 

ii 

II 

Phase £ " 

II 11 

II ll 

bs 










X cs 

A 

II 

H 

li 

ii 

II 

Phase £ " 

II li 

il H 



" rotor Phase £ " 

" " Phase b " 

" 11 Phase c " 


A leakage inductance of stator Phase £ 

L-. A " " " " Phase b 

lbs = — 

L, A " " " " Phase c 

lcs = — 


L 

L 

L 


lar 

Ibr 

ler 


A 

A 

A 


Phase £ 
Phase b_ 
Phase c 


ii 


n 


Fig. 2.1 shows the three phase axes of the stator and rotor 
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Stator 


|h as 


"1 cs 


c s f: 


as 




x 


lbs 


c b s 




i i 


f 


‘lar 


R cr '/ Rar 

'1 cr ^ I 1/ 

> ox — — '" f ' J ' \ 




Rotor 


X L lbr 

*\ 


b„ 


Fia. 2.1 


i 


a: 


r 


de 

dt 


Because of the symmetrical and balanced windings: 

R as = R bs = R cs = R s (stator resistance/phase) 

R ar = R br = R cr = R r (rotor resistance/phase) 

L las = hbs = hcs = L ls ( stator leakage inductance/phase) 
L lar = hbr = her = hr ( rotor leakage inductance/phase) 
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The basic equations can be shown in matrix form as 
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where the variables written in vector form are 


as 
v bs ‘ 


cs 


(2.9a) 


iv. 


ar 


'br 


cr 


(2.3a) 


as 

, V 


i i 


cs 


(2.9b) 


ar 

br 

cr 


(2.8b) 


X 

i as 

|X = X, 

I s bs 


cs 


(2.9c) 


ar 


br 


cr 


(2.8c) 


Next, the equation for the total flux linkage for each stator and rotor 
phase and the air gap magnetic field is written as follows 


X = Li 


( 2 . 10 ) 


where 

x as = 'as + L ms 'bs cos (2 ' /3 > * L m s 'cs cos <- 2,/3 > 

+ L mr 'ar cos 9 * L mr 'br cos (9+27/3 > + >- mr i cr cos (6-2r/3) 

( 2 . 11 ) 

X. = L mc i, r cos (- 2 tt/ 3 ) + L p i Kc + L . i_ . cos (2ir/3) 
bs ms as ' ' ms bs ms cs 

+ L m „ i cos (9-2tt/3) + L i. cos 9 + L _ i* cos (9+2n/3) 
mr ar mr br mr cr 


( 2 . 12 ) 
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and 


X cs - Has 'as cos < 2 "' 3 > + L ms 1 bs cos <- 2t/3 > + L ms 'cs 

L mr 'ar cos (5t2,r/3) * L „r 'br cos < 9 - 2 ’ /3 > * L mr 'cr cos 9 

(2.13) 

V * L m 'as cos 9 * L ms 'bs ^ (9 - 2r/3> * L ms 'cs cos (9t2 ' /3 > 
L „r 'ar + L mr 'br cos <- 2 ^ 3 > + L „r 'cr cos <- 2 ’' 3 > 

(2.14) 

A br * Sns 'as (:t27/3 ) * l ms 'bs cos 3 + L n,s 'cs cos (S - 2 ' /3) 


L i cos ( 2^/3 ) + L i. + L i cos (-2V3) 
nr ar mr br mr cr 


(2.15) 


X = l„ c i, c cos (^- 2^/3 ) + L__ i Kc COS ( 9 + 2 TT /3 ) + L i‘ cos 5 
cr ns as v ' ' ms DS IDS CS 

( 2 . 16 ) 


cr ns as 
^nr ^r 


L «- cos + L mr ’ br cos (2t/3 > + L mr V 


L mS 4 magnetizing inductance/stator phase of air gap magnetic field 
L mr 4 magnetizing inductance/rotor phase of air gap magnetic field 


Equations (2.11)-(2.16) can be written in the matrix equation form of (2.17). 
For simplicity, (2.17) can be divided into four 3X3 as follows 



L ms 

-1/2L 

ms 

- 1/2L ™sl 


h - 

- 1/2L n s 

*"ms 

-'' 2L J 

(2.18a) 


-T^Lms 

-1/2L 

ms 

*"ms j 




cos9 

cos( 9 - 2 h 73 ) 

cos( 9 + 2 ti/ 3 ) 


"mr 


mr 


mr 


cos (9 +2tt/3 ) 
cose 

cos(9-2ir/3) 


L ms cos( 9 + 2 r/ 3 ) 

L mr cos( 9 + 2 Tr/ 3 ) 

L_ cos9 
mr 


(2.18b) 



L c C0S9 
ms 

L ms C0S ^ 8+27r/3 ^ 

L ms cos( 9 - 27 r/ 3 ) 


L ms c °s (9-271/3) 

L m .cos9 

ms 

L ms cos(9+27i/3) 


L ms cos ( 0+27T / 3 ) 

L ms cos ( 6 - 271 / 3 ) 

L m .cos9 

ms 


( 2 . 18 c) 



as 


'bs 


cs 


ar 


br 


cr 


L mr 

ms 

-1/2L 

ms 


-1/2L 

ms 

L ms 


-1/2 L 

ms 

— 1/2 L 

ms 

L ms 

L mc cos0 

ms 

L ms cos ^" 27r/3 ^ 

L ms cos (°+2 ti / 3 ) 

L rns cos (6+2Tr/3) 

L mr .cos0 

ms 

^ms COS ( n ” 27T ^ 3 ) 

LmsCOsCQ-Zir/S) 

L ms cos(0+2Ti/3) 

L m ^cosO 

ms 


L costl 
mr 

L mr cos ( 0+ 2i r / 3 ) 

L mr cos ^ 0-2ll/3 ^ 

L| nr cos (0- 2 it / 3 ) 

LcosO 

mr 

L mr c °s(0+ 2 n/3) 

L mr c °s (n+2n/ 3 ) 

L mr cos(0-27i/3) 

L cos0 
mr 

^mr 


-1/2L 

mr 

'^mr 

^mr 

-1/2L 

mr 

'^mr 


^mr 


(2.17) 
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L mr 

mr 

-l/2L mr 

mr 

-1/2L 

-l/2L mr 

mr 

L 

mr 

-1/2L 

-1/2L 

mr 

-1/2L 

mr 

L m 

mr 


(2.13d) 


Then (2.17) can be written 


y . ! L i y v 

A r, 4 L 4, V| 

Since this induction motor has no neutral connection 


(2.19) 


Vs + Vs + Vs “ 0 
Vr + Vr + Vr * 0 


Therefore 


Vs ~ "Vs ’ Vs 
Vr " _1 br " Vr 


(2.20a) 

(2.20b) 

(2.21a) 

(2.21b) 


If these equations are substituted into (2.17), matrixes Lj and become 


3/2 L 0 

ms 


( 2 . 22 ) 


0 3/2L„ 


3/2L mr 0 


0 3 / 2 ^ 

0 0 


mr 

0 3/2L 


(2.23) 
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It is more convenient to write the resistance matrix and 


inductance matrix of (2.7) 

as 

follows 


! R 

0 

1 


L, 

i r i * ; s 


| 

(2.24) 

,i ; Is 
L il - ■ 

■ 0 

R , 

r ' 


|0 

where 






R 

0 

0 


L. 

s 




Is 

; ' R s ' = '° 

R 

s 

0 

(2.24a) 

i L ls'- 0 

: 0 

0 

R s 


’ 0 

R 

r 

0 

0 


L lr 


R r « 0 

R 

r 

0 

(2.24b) 

Kir '’. 0 

1 

i° 

0 

R r 


1° 


lr 


0 ( 

■IS C 
0 l 


0 c 


■lr 


- 


By substituting (2.19), (2.24) and (2.25) into equation (2.7), 
motor equations can be written 

[V s ] = [R $ ] [T s ] + [L 1 ] d/dt [T s ] + d/dt [L 2 ] [T r ] + 
tf r ] = [R r ] (t r ] + d/dt (L 3 1 [t $ ] + [L 4 ] d/dt [T r ] + 


Equations (2.26) and (2.27) can be written in matrix form 
!! I[R _1 T. |[L,J. ] T. lILj] £L 2 3 ; 

[L 3 ] [L 4 ] r r 



[R s ] T s 

“•is 1 


-t 


. + 


P 

i + p 


[ R r i 

[hr> 


Vi 


the leakage 


(2.25) 


(2.25a) 

Is 


(2.25b) 


the induction 

[L ls ] d/dt [t s ] 
(2.26) 
[L lr ] d/dt [T r ] 

(2.27) 


(2.28) 


where 


P = d/dt 
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3.0 TRANSFORMATION TO STATIONARY d-q REFERENCE FRAME 
In this report, the following symbols are used: 

d-q => stationary reference frame (fixed with respect to 
the stator) 

D-Q => rotating reference frame (rotating at rotor speed - 
fixed with respect to the rotor) 

3. 1 Definition of d-o Stationary Reference Frame 

For the d-q stationary reference frame, the d-axis is m line with the 
a-axis of the stator as shown in Fig. 3.1. Then, 


v d = v as + v bs cos (120°) + v c cos (-120°) 
Vq = sin (120°) + v c sin (-120°) 


v d 


1 

-1/2 

-1/2 

V q 


0 

/3/2 

-/3/2 


(3.1) 


3.2 Definition of D-Q Rotating Reference Frame 

The D-Q reference frame rotates relative to the d-q stationary reference 
frame, where 0 is the angle between the d-axis and the D-axis and 

e ' /<- r Ur = f 


where ay is the rotor speed in radians/second. This is shown in Fig. 3.2. 


cos 0 -sin 0 


sin 0 


cos 8 


(3.2a) 
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Fig. 3.1 

d 

* O N ^ 



Fig. 3.2 



223 


V D 


cos 9 sin 9 


v d 

v n 


-sin 9 cos 9 


v 

Q 




q! 


(3.2b) 


The v is rotating clockwise at angular velocity uy. 

3.3 Induction Motor d-q-Q Variables 

For simplicity in the definitions of the d-q and D-Q reference frames, 
only two axes were considered. In general, it is desirable to consider three 
coordinate axes for two reasons. Matrix transformation methods can be used 
to transform variables from one reference frame to another, but only if the 
transformation matrices are square and nonsingular matrices. Thus, for 
three phase situations with three phase variables, it is necessary to con- 
sider transformations to coordinate axes where three variables exist. For 
balanced and symmetrical cases, it is often sufficient to consider only d-q 
variables. However, the three d-q-0 variables are necessary in the most 
general situations. 

In the remainder of this report we are concerned with the d-q-0 reference 
frame. Again, this is a reference frame fixed with respect to the stator. 

The O-variables will not only make it possible to use matrix transformation 
techniques, but they also will have important physical significance in many 
cases. 

For this report, the d-q-0 transformation from stator three phase voltages 
or currents is defined as follows 


v ds 


V 


qs 

= 3 

CO 

o 

> 



1 

0 

1/2 


cos 2tt/3 
sin 2tt/3 
1/2 


cos (-2 tt/3) 


v as 

sin (-2 tt/3) 


< 

cr 

to 

1/2 


v cs 


(3.3a) 
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^s 


i 

cos 2 tt/3 

cos(-2tt/3) 


^as 

V 


0 

sin 2-/3 

sin(-2m/3) 


V 

1 

. 

^s 


1/2 

1/2 

1/2 

1 1 cs ! 


Equations (3.3a) and (3.3b) transfern the three phase stator line-to- 
neutral voltages or line currents to d-q-0 coordinate voltages or currents. 
It should be noted that 


'os = 1/3 + * 


as 


1 bs 




(3.4a) 


and thus ig s is zero for all three wire machines where no neutral is present. 
This is a consequence of Kirchoff’s current law. However, 


v 0s " 


(v 


as 


+ v 


bs 




(3.4b) 


is zero when balanced and symmetrical three phase 1 ine-to-neutral voltages 
are involved, but it is not zero in general. 

Equations (3.3a) or (3.3b) apply to the transformation of stator three 
phase to d-q-0 quantities for all circuit variables including flux linkages 
as well as the voltages and currents. 

The d-q-0 transformation for rotor quantities is defined as follows 


V dr 


V 

A 2 
3 

< 

O 

“I 



cos 6 cos(8 + 2 tt/ 3) 
sin 9 sin(9 + 2m/3) 
1/2 1/2 


cos(0 - 2m-/ 3 ) 
sin(6 - 2tt/3) 
1/2 


ar 


'br 


cr 


(3.5) 


Equation (3.5) transforms the rotor quantities, which are given in a 
three phase coordinate frame rotating with angular velocity to the d-q-0 
reference frame fixed with respect to the stator. The same transformation 
as in (3.5) is used to transform currents and flux linkages. Fig. 3.3 illus 
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Fig. 3.3 
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trates this relationship. The 0-axis is orthogonal to the d and q-axes, and 
it may be considered directed into the paper. 

A number of additional transformation relations are given as follows 


for stator variables 


where 


.-V 

' V S 

= ;a , " 1 

i H S 1 

' V dqo' 

(3.6a) 

: V 

= !a !_1 

|A s' 

'i S 1 
1 'dqO 1 

(3.6b) 

;*s 

= !a 1-1 

l i 

'X s ' 
1 A dqO ' 

(3.6c) 

and 

for stator 

variables 


i 

lV r 

= lApj* 1 

I v dq0' 

(3.7a) 

i’r 

= ! A I' 1 

1 r 1 

r 

' ^dqO' 

(3.7b) 


= 1 A 

1 n p 1 

r 

> X dqO l 

(3.7c) 



1 

-1/2 

-1/2 

A I = — 
H S‘ 3 

0 

/J/2 

-/J/2 


1/2 

1/2 

1/2 


1 

0 

1 

i-l 

1 

-1/2 

/J/2 

1 


(3.8) 


(3.9) 


-1/2 -/J/2 
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i\i • f 


cos 9 cos(9 + 2 t t/3 ) cos(9 - 2ir/3) 


sin 9 sin(9 + 2tt/ 3 ) sin(S - 2-/3) 


1/2 


1/2 


1/2 


(3.10) 


cos 


Sin r 


1! 


I A r r 1 * j cos ( 9 + 2r/3) sin(5 + 2-/3) 


(3.11) 


cos(9 + -2-/3) sin(r - 2tt/3 ) 1 ! 


From (2.7), the motor voltage equations are 


l V cl 


IV 


R c nj + P Lie IT,; + P,XJ 


S ' S ' 


-Is >'s 


R r “VI + P L lr IV + ^r' 


(3.12) 

(3.13) 


Substituting (3.5) and (3.7) into (3.12) and (3.13) respectively gives 


IV I v dq0 I = P c i A e! I ^ HnO 1 + L le P ' A c ' I 1 rl n fl I + 


s 1 s 


dqO 1 u ls 


dqO 


p lAj' 1 u s 


dqO 1 


I A r I 1 I v dq0 I = R r l A rl 1 I i dqol + L lr P i A r' 1 I’dqol + 


P | A. 


■1 IX r 1 
r 1 ' aqO 1 


(3.14) 


(3.15) 


Equation (3.15) can be expanded as 


I A r I 1 I v dq0 I = R r l A rl 1 I i dq0 I + L lr l PA r' l*dqol + 


h r 1‘rl' 1 p tl^nnl> + < p IM' 1 ' 


dqO 


r 
dqO 1 


+ lArl" 1 tP hdqoll 


Equation (3.14) need not be expanded because the matrix | A g | “ 1 is constant 


(3.16) 
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It can be shown that [1] 





(3.17) 


Substituting (3.17) into (3.16) gives 


' v dqir * R r X ' 1 * hr ^X ' 1 tP i A r i } I i 


+ l a ”1 !p !i r i' + * a ■ "1 ;o i a •' ’a 1 ‘ 1 i . r 

L lr A r |P 1 1 dqO 1 1 1 |A r' |A r‘' 'V ' 'dqO 


+ !A. 


1-1 r D ’-i r i’. 


r 1 


jp v 
lP ■' dqO 1 


(3.18) 


Multiplying both sides of (3.14) and (3.18) by j A J and |A r | respectively 


gives 


iv 


s 

dqO 1 

" R s i 

1 1 dqO' + L ls P ! 1 dqO' + P jX dqO 

r 

d qO 1 

= R r : 

’dqO! + L lr P ^dqo' ' L lr * P 


+ P 

1 / dqO^ ” *- D ! A r^ 1 A r 1 |X drO 


(3.19) 

i r ’ 

dqO 1 

(3.20) 


The expressions !P J A p j } jA^!"* in (3.20) can be evaluated using matrix 
algebra, which is shown in Appendix I. The result is 

0 -1 oi 

1 

-1 


(p |A r ;i iA r , 


1 0 0 
0 0 0 


d0 

dt 


(3.21) 


d9 , 

dt s 


where 
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and 


I H | 


0 1 0 

-10 0 
0 0 0 


Substituting (3.21) into (3.20) gives 

! v dq0‘ = R r •*dq0 i + L lr P ^dq0' + L lr w r ;H| 1 ’ dqO • 

+ P 1 ' - 'dqO' + °r i H ' 1 X dqO 1 

Equations (3.19) and (3.23) are the motor voltage equations in the d- 
frame. 

Now, the flux linkage equations are expressed in terms of d-q-0 
as follows. Matrix (2.19) can be written 


i>v = ! i»i i ;i s i + |l 2 ! ji f j 
I* 


i _ 1 1 i i ,• ! + 1 1 ' I ,• i 

1 ” 1 ^3 • 1 1 s + 1 : 1 ' 


- 4 ; iyi 


Substituting (3.6c) and (3.6b) into (3.24) and (3.7c) and (3.7b) into 
gives 


i\! i*s q0 ! * mi iaj ' 1 mlj * mi mr 1 i> r 


1 • 1 s 1 


-1 


dqO 1 
■1 i,s 


■2' 1 r 1 


dqO 1 


i\i iA d r q0 i • il 3 i iA s r‘ n5 q0 i ♦ il 4 : iA r r 1 i. dq0 , 


Multiplying both sides of (3.26) by A $ and (3.27) by A p gives 


s 1 
dqO' 

■ |A S ! 

141 

M' 1 

I'dqol * 1 A s 1 

141 iA r r' 

1 - r 
1 1 dqO 

r I 
dqO' 

* |A r i 

i L 3 1 

i\r‘ 

Hdqol * ! fl rl 

|L 4 I 1a/* 

I'd* 


(3.22) 

(3.23) 
q-0 

current 

(3.24) 

(3.25) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 
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where 

|a s I |l 2 | 1a,.!' 1 


3/2 L 


mr 


0 0 
1 0 
0 0 


(3.30a) 


i 

i 


A r> 




n 


3 / 2 L ns 1° 


0 

1 

0 


°i 

o! 


I L-j | ; 


V 1 ■ 3/2 L «,s 


I 1 

i° 

!o 


0 

1 

0 


0, 

I 

o! 




3/2 L 


mr 


;i 

!° 

lo 


0 

1 

0 


o! 

i 



(3.30b) 


(3.31a) 


(3.31b) 


as shown in Appendix II. Substituting (3.30) and (3.31) into (3.28) and 
(3.29), respectively, gives 


s 1 
dqO 1 

■ 3/2 L ms l £ l 

1 .s I 
■’dqO 1 

+ 3/2 L 

mr 

IE! 

'i r 1 
1 dqO' 

(3.32) 

r i 

dqO 1 

■ 3 ' 2 L ,ns 

1 1 dqol 

+ 3/' 2 E 

mr 

|E| 

li r 1 
1 dqO 1 

(3.33) 


where 


1 0 0 


E 


0 


0 


1 0 
0 0 


(3.34) 


Equations (3.32) and (3.33) can be written in matrix form using (3.34) 
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Equations 

yielding 


where the 
matrix in 
form 


'V 


L 

0 

0 

L 

0 

0 


1 , ! 

! 

ds 


ms 



mr 



! 

ds< 


X 


0 

L „ 

0 

0 

L 

01 

1 

1 1 

i 


qsj 



ms 



mr 

1 

qs i 

! 

A n 

Os 

| 

= 3/2 

0 

0 

0 

0 

• 

0 

°! 

1 

1 

^Os ! 

(3.35) 

' 

'dr 

1 

1 

E ms 

0 

0 

E mr 

0 

0| 
j 1 

1dr 

> ! 

1 

t 

0 

L „ 

0 

0 

L 

l 1 

0; ■ 

1 

1 ' 


qP ! 

1 

1 

ms 



mr 

1 1 

1 1 

qr 1 

i 

1 

Or 

i 

1 

1 0 

0 

0 

0 

0 

O' i 
1 1 

v; 

(3.32) and (3.33) 

can 

be substituted 

into 

(3.19) and (3.23) 


'dqo ■ R s ■ 1 dqo + L i s p :4o + 3/2 p L - : e : i 1 


.s 


"ms 


dqO' 


+ 3/2 P L mr ^ 1 1 dqO 


(3.36) 


< r 
|V dqG 


R r 1 1 d q 0 1 + L lr P 1 ’dqO + L lr > ' H 1 1 dqO : 


+ 3/2PL ms ' E ' >dq 0 : 3/2 P E mr X ^qo' 

+ 3/2 > M L ms = E i • "*dqO * + 3/2 "r M L mr 


! c i I ,• r I 

dqO ' 


'Pi T -• 

l E . I 1 


(3.37) 


motor equations are written in terms of d-q-0 currents. The final 
d-q-0 coordinates is found by writing (3.36) and (3.37) in matrix 


v ds 


R s +PL s 

0 

0 

-PL 
2 L mr 

0 

0 


y 

v qs 


0 

R s +PL s 

0 

0 

fPl_ r 

2 mr 

0 


V 

v 0s 


0 

0 

R s +PL ls 

0 

0 



'os 

v dr 

= 

yPL 
2 ms 

f JJ r L Lms 

0 

R +PL- 
r r 

“r L r 

0 


^dr 

v 0 „ 

qr 


3 1 
iPVSns 

f PL ms 

0 

-u L 
rr 

V’ L r 

0 


'qr 

v 0r 


0 

0 

0 

0 

0 

R r +PL lr 


V 


(3.38) 
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where 

L s ’ hs + 3 ' 2 L ms 
L r ’ hr + 3 ' 2 L mr 

For the symmetrical three phase induction machine with wye-connected windings 
having no neutral connections, the O-components disappear so that d-q-0 
coordinates become d-q coordinates, -educing (3.351 to 


1 ‘ 
, ds 

l • 

;v pL s 

0 

-PL 
2 mr 

0 i 

„ 1 

' ds 

lv 1 
i qs | 

< 0 
_ 1 

R S+ PL S 

0 

|PL i 
2 mr 

!’«• 

'V. 1 
! dr 

'hi 

, 2 ns 

2 |_ 
2"r ms 

R r +PL s 

'j L i 

r r 

i 

1 i , 

, dr 

1 V ' 

, R r 

i-L L 
, 2 r ms 

-PL 

2 KL ms 

"*V L r 

R +PL i 
r s . 

! 1 qr 


4.0 RESULTS 

Matrix methods have been used to derive the d-q equations for the induc- 
tion machine. The resulting matrix equation (3.39) is much simpler than the 
conventional three phase representation. Only four voltage and current vari- 
ables are involved instead of six. Also, the representation of (3.39) is very 
useful in control analysis using state variable methods and optimization 
techniques . 
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APPENDIX I 
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Equation (3.21) is derived as follows. 

From (3.10), 

-sin 3 -si n(5 + 2-/3) 

P jA r ; = 2/3 cos 9 cos(0 + 2-/3) 

0 0 

Therefore, ."P ,A r ' ; } , A^ ' * can be expressed 

j - s i n 3 -sin(5 + 2-/3) -sin(5 - 2-/3) 

:p |A r : :A r ; _1 = 2/3 |cos 5 cos(e + 2-/3) co S (3 - 2-/3) ; ^ x 

I 0 0 o' 

Jcos 3 sin 3 1 I 

|cos(9 + 2-/3) sin(3 + 2-/3) l| 

| cos (5 - 2r/3) s i n ( 3 - 2^/3) l| 

Multiplying the first row of P jA r | by the first column of j A r ; 1 yields 

C 

-sin 9 - cos 9 - sin (S + 2-/3) cos(3 + 2tt/3) - sin(3 - 2r/3) cos(9 - 2-/3) 

= - |l/2 sin 29 + 1/2 sin 2 (9 - 2 tt/ 3) + 1/2 sin 2 (9 + 2-/3 ) J 

= - 1/2 j si n 29 + sin (29 + 2-/3) + sin (25 - 2 tt/3)', 

= - 1/2 | s i n 29+2 sin 29 (cos 2 t/3 ) ! 

= - 1/2 J si n 29 - sin 29j = 0 

Next, multiplying the first row of P ] A p j by the second column of |A r ’ * gives 

-sinQ 2 - sin 2 (9 + 2-/3) - sin 2 (6 - 2 tt/3 ) 

= - |l/2 - 1/2 cos 29 + 1/2 - 1/2 cos 2 (9 + 2n/3) + 1/2 - 1/2 cos2(9-2tt/3) 

= - 3/2 + 1/2 | cos 29 + cos (29 + 2 tt/ 3) + cos (20 - 2 tt/3) 

= - 3/2 + 1/2 | cos 29 + 2 (cos 29) (cos 2ir/3)| 

= - 3/2 + 1/2 | cos 29 - cos 29 1 = - 3/2 


- s i n ( 9 - 2-/3 ) | 

cos(5 - 2-/3) * d5_ 

! dt 

A * 
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Finally, the first row of P ; A p [ is multiplied by the third column of 
| A r | “ 1 which yields 

- sin e - sin (3 + 2-/3) - sin (5 - 2-/3) = 0 

Thus, the first row of (3.21) has been derived. The second and third rows 
of (3.21) are derived in a similar manner. 
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Equation (3.30) is derived as follows, using JA j from (3.8), |L_] 

5 C. 

from (2.18b) and A r from (3.10). 

|A S ! [L 2 : |A r ; = 



i 1 

-1/2 -1/2 j ; cos e 

cos(9 

+2-/3) 

cos(- 

-2-/3) 

= 2/3 

'o 

i 

*7/2 -*'3/2 J L Tr | cos ( •?- 2-/3 ) 

cos e 

cos(e 

+2-'3) 


;i/2 

1/2 1/2 1 , cos ( r+2-/3 ) 

cos (r 

-2-/3) 

CO 

s - 


Jcos e 

sine 1 ] 





X 

| COS (e 

+2-/3) sin(9+2-/3) lj 






i COS (r 

-2-/3) sin( c -2~/3) 1 \ 





rst, ;a s 

1 ' L ' 
1 * L 2 ' 

is determined. Multiplying the 

fi rst 

row of \ 

A s ; by 

the 

second and thi 

rd columns of 





cose 

- 1/2 

cos • 9-2 tt/ 3) - 1/2 cos (e+2ir/3) 





cose 

- 1/4 

cos 9 - v7/4sin9 + l/4cose + /J/4s 

;in9 




(1 + 

1/4 + 

1/4) cos 6 = 3/2 cos 6 




(1) 

cos (9+2-i/3) 

- 1/2 cos9 - 1/2 cos(6-2t/3) 





-1/2 

cose - 

3/2 sin9 - 1/2 cos0 + 1/4 cose 

- /7/4 

sine 



(-1/2 

-1/2 

+l/4)cos6 - 3 (v7/4) sine = - 3/4cos6 - 

3(v7/4) 

sin 9 

(2) 

cos( 

-2 /3) 

- 1/2 cos( +2 /3) - l/2cos 





COS (9 

- 2tt/3) 

- 1/2 cos (0+2tt/3) - l/2cose 





-l/2cos9 + 

✓7/2 sine + 1/4 cose + /J/4 sine 

- 1/2 

cose 




(-1/2 - 1/2 + 1/4) cos9 + 3(vT/4) sine = -3/4 cos9 + 3(/J/4) sine (3) 

Multiplying the second row of | A g | by the first, second and third columns of 
1 4 1 gives 

3/2 cos (0-2tt/3) - ✓3/2 cos8(tt+2tt/3) * 

- /3/4 cose + 3/4sin 6 + /7/4 cose + 3/4 sine 

= 3/2 sin 9 (4) 
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/J/2cose - /J/2e(-l/2 cos 9 + /J/2 sin9) 

/372cos0 + /J/4 cos0 - 3/4 sine 

= 3(/3/4) cos0 - 3/4 sin9 ( 5 ) 

/J/2 cos (0+2 t:/ 3) - /J/2cos9 
- 3/4 cos 9 - 3/4 sine - /T/2ccse 

= -3(v"J/4) cosS - 3/4 sine (6) 

Multiplying the third row of }A S ! by the first, second and third column of 
|L 2 | gives 

1/2 cose + cos (e-2r/3) + cos (e+2-/3) 

1/2 cose - l/2cose + vT/2 sine - l/2cose - /T/2 sine 


1 - 1/2 - 1/2 = 0 


(7) 

V2 I | = 

0 


(8) 

1/2 ! I = 

0 


(9) 



-1 


The resulting A, L 

J s t. 

, is now multiplied by A r 

' : 


3/2cos 

-3/4cos - 3( 3/4)sin 

-3/4cos + 

3( 3/4)sin 

2/3 L mr 3/2sin 

3( 3/4)cos - 3/4sin 

-3( 3/4)cos 

; - 3/4sin 

0 

0 

0 


cos 9 

sine 

1 

| 

X cos{9+2tr/3) 

sin(9+2^/3) 

1 

(10) 

cos(0-2tt/3) 

sin(9-2Tt/3) 

1 


The first row of |A | 

|L 2 | is now multiplied by 

the first, second and third 

columns of |A 
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3/2cos 2 0 + (-3/4cos0 - 3{/T/4)sin9) (-1/2cos9- v7/2sin9) 

+ (-3/4COS3 + 3(/T/4)sin6) (-1/2cos3 + /J/2sin9) 

3/2cos 2 8 + /I/8cos 2 9 + 3(/T/8)sin9cos9 + 3(v7/8)sin0cos3 + 9/8 sin 2 3 

+ v7/8cos 2 3 - 3(»7/8)sin6cos9 -3( v7/8)sin0cos9 - 9/8sin 2 9 

3/2cos 2 a + 3/8cos 2 9 + 3(>'7/8)sin?cos9 + 3(*7/8)sin9co:9 + 9/8sin 2 9 

+ 3/8cos 2 9 - 3(/T/8)sin9cos9 + 3(/T/8)sin5cos5 + 9/8sin 2 9 

3/2 + 3/8 + 3/8 = 9/4cos 2 9 + 9/4sin 2 3 = 9/4 {}]) 

3/2sin9cosr + (-3/4cos9 - 3(»7/4)sin9)(-l/2sin9 + /J/2cos9) 

+ (-3/4cos9 + 3(v7/4)sin9) (-l/2sin3 - >7/2cos9) 

3/2sin9cos9 + 3/8cos3sin9 - 3(v7/8)cos 2 9 + 3(v7/8)sin 2 9 
-9/8sin9cos9 + 3/8cos9sin9 + 3(v7/8)cos 2 9 - 3(v7/8)sin 2 3 


-9/8sin3cosr = 0 

02) 

= 0 

(13) 


The second row of |A S [ • L 2 • is now multiplied by the first, second and third 
columns of |A r }~^ 

3/2sin9cos9 + ((3v7/4)cos3 - 3/4sin3)(-l/2cos9 - 3/2sin0) 

‘ + (-3(v7/4)cos9 - 3/4sin9)(-l/2cos0 + /3/2sin0) 

3/2sin6cos9 + -3(v7/8)cos 2 9- 9/8sin9cosS + 3/8sin0cos3 
+ 3(/T/8)si n 2 e + 3(/3/8)cos 2 9 - 9/8 sin3cos9 + 3/8sin9cos9 - 
3(v / 7/8)sin 2 9 

2 

3/2sin0cos9 - 3(/3V8)cos 9 - 9/8sin0cos0 + 3/8sin0cos9 
+ 3(/I/8)sin 2 3 + 3(/3/8)cos 2 9 - 9/8sinecos0 + 3/8sin0cos0-3(/3/8)sin 2 3 
= 0 04 ) 

3/2sin 2 0 + (3(/I/4)cos3-3/4sin0)(-l/2sin9+ /I/2cos8) 

+ (-3(>7/4)cos0 - 3/4sin0)(-l/2sin0 - /7/2cos0) 
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? 2 2 
3/2sin 0 + (3(/J/8)sin9cos9 + 9/8cos 9 + 3/8sin 9 - 

3(/J/8)sin0cos9 + 3(/I/8)sin9cos9 + 9/8cos 2 9 + 3/8sin 2 0+ 3(v / J/8)sinecose 

3/2sin 2 9 - 3(/I/8)sin0cos9 - 9/8cos 2 9 + 3/8si n 2 0 - 3(/3/8)sin0cosS 
+ 3(/3/8)sin0cos9 + 9/8cos 2 9 + 3/8sin 2 0 + 3(/J/8)sin9cos9 
(9/4si n 2 9 + 9/4cos 2 0) = 9/4 

The third row of |A $ j [Lgj, multiplied by the first, second and third columns 
of |A f 1 , gives zero in each case. 

Therefore, 

I 

3 / 2L mr 0 0 

0 3/2L mr 0 

0 0 0 1 

Equation (3.31) can be derived in the same manner. 
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